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PILE NEUTRON CAPTURE CROSS SECTIONS OF Co*™* 
(10.5 MIN.) AND Co® (5.25 YR.)! 


By F. Brown, J. L. WoLFson, 
AND L. YAFFE? 


ABSTRACT 

Activity of approximately 105 min. half-life has been detected in samples of 

Co** irradiated at high neutron fluxes in the NRX reactor. By a study of the 

yield as a function of neutron flux and time of irradiation this activity was 

demonstrated to be due to Co* (105 min.) formed by neutron capture in Co®, 

The pile neutron capture cross sections were found to be 95 + 45 barns and 6.0 

+ 1.4 barns for the Co®* (10.5 min.) and Co® (5.25 yr.) isomers respectively. 

I. INTRODUCTION 

In determinations of the thermal neutron capture cross section of Co®® (1, 
3, 11, 14) it has been assumed that the thermal neutron capture cross sections 
of the Co®* (10.5 min.) (7) and Co (5.25 yr.) (12) isomers are negligibly 
small. Large values of the capture cross sections of these isomers would have 


a pronounced effect on the value for the Co®® capture cross section determined 
by either the activation method or the pile oscillator method. In the former 
case the value obtained would be spuriously low, and in the latter case spuri- 
ously high. The pile oscillator method might be expected to give a more 
reliable value since the neutron fluxes employed are usually low. The fact 
that the results obtained by both methods are in agreement provides some 
evidence that the Co®* (10.5 min.) and Co® (5.25 yr.) capture cross sections 


are both small. 

This conclusion is supported by the results of the present experiment, in 
which measurements have been made of the pile neutron capture cross sections 
of the two isomers for formation of Co® (105 min.). The values obtained were 
95 + 45 barns and 6.0 + 1.4 barns for the short and long lived Co® isomers 
respectively. 

These cross sections were determined from a study of the production of 
Co in Co®* samples irradiated in the NRX reactor, and information as to 
the disintegration of Co®™ was therefore required. The half-life of Co™ has 
been measured by several observers and values reported ranging from 99 min. 
to 108 min. (6, 9, 10, 13). Parmley et a/. (9) found the 8-ray end point to be 
' 1.3 + 0.1 Mev. but detected no y-rays. Smith et al. (13) found two groups of 
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B-rays and in addition detected a y-ray of approximately 0.5 Mev. They 
conclude that Co decays by negative B-ray emission, 55% of energy 1.42 
+ 0.02 Mev. to the ground state of Ni®, and 45% of energy 1.00 + 0.02 Mev. 
to an excited state in Ni® at 0.42 Mev., followed by emission of a y-ray of 
the same energy. 

In the present work the half-life has been taken to be 105 min. since the 
data obtained agree best with this value. Nevertheless the data do not provide 
an accurate measure of the half-life and would not be inconsistent with a 
somewhat shorter or longer value. The errors assigned to the cross sections 
take into account the uncertainty in the lifetime. The disintegration scheme 
proposed by Smith et al. (13) has been adopted. Choice of a simple beta 
spectrum, as indicated by the results of Parmley et al. (9), would have the 
effect of reducing the measured values of the cross sections by about 10%. 


II. APPARATUS AND METHOD 

Samples of spectroscopically pure Co** sponge or wire of known mass were 
irradiated in the “rabbit’’ of the Chalk River reactor for periods of time 
ranging from 1.5 to 12 hr. After irradiation, a sample was dissolved in nitric 
acid, made up to a known volume, and an accurately measured aliquot de- 
posited on thin aluminum foil. The nitric acid was evaporated under heat to 
leave a thin layer of the nitrate about 1 cm. in diameter and about 0.3 mgm./ 
cm.” thick. The sample was then covered with a solution of polystyrene in 
benzene, which after drying left a thin covering layer about 0.1 mgm./cm.? 
thick. 

The activity was counted with an anthracene scintillation spectrometer, 
consisting of a piece of anthracene approximately 1.3 cm. X 1.3cm. X 0.3 cm. 
thick applied to the face of a 5819 photomultiplier with DC 200 silicone fluid. 
The anthracene was covered with an aluminum reflector 0.15 mgm./cm.? 
thick. The assembly was housed in a light-tight container with an aluminum 
window about 0.5 mgm./cm.? thick. The sources were placed 4 cm. above the 
window, resulting in an effective solid angle of about 0.9%. 

The spectrometer was calibrated with a source of Cs!7 (37 yr.) of known 
activity, using a single channel pulse height analyzer. The activity of the Cs!87 
source was determined by comparison with a standardized source which had 
been measured in a 4m proportional counter. The position of the 624 kev. 
line due to the internally converted 661.6 kev. y-ray (8) provided the energy 
calibration. The intensity of the line, i.e. the integrated number of electrons 
per second counted in the line, along with the conversion coefficients (4) and 
8-branching ratio for Cs'*’ (5) gave the solid angle. Width of the conversion 
line was 15% at half maximum intensity. 

A typical experiment involved calibration of the spectrometer, followed by 
replacement of the Cs'*’ source with the cobalt source. The lower discriminator 
of the pulse height analyzer was then set to a pulse height corresponding to 
410 kev. and all pulses of heights larger than this were taken from the “in- 
tegral’’ output of the analyzer and fed into a scaler. The output of the scaler 
was arranged to record pips every 1,000 or 10,000 counts, depending on 
counting rate, on the chart of an Esterline Angus recorder. Counts were 
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recorded in this manner for approximately 20 hr., after which the spectrometer 
was calibrated again to determine whether drifts in the electronic apparatus 
had altered the effective bias setting. Corrections were applied to the data 
where changes in the position of the Cs'’ conversion line indicated bias 
changes larger than 1%. Fig. 1 illustrates the type of decay curve obtained. 

The high background is due to the 1.17 and 1.33 Mev. y-rays from Co®. 
In addition there is the 1.56 Mev. 6-ray transition of intensity 0.28% (2) 
from the 10.5 min. isomeric state in Co, which is detected in the initial 
period of counting. However the 320 kev. 6-rays from the 5.25 yr. Co® isomer 
are not counted since the discriminator was set to reject pulses of height 
lower than that corresponding to 410 kev. After subtraction of the background 
there remains a component with a half-life of roughly 105 min. In every case 
this component has been fitted with a line of slope corresponding to a half-life 
of exactly 105 min. This procedure assumes that the activity is due to Co®. 
Justification for this assumption is given below. The fraction of the 6-spectrum 
counted assuming the disintegration scheme of Smith et al. (13) is 538%. 

In all, eight trials were made: one at an irradiation time of 1.54 hr., three 
at an irradiation time of 3 hr., one of which was done at reduced neutron flux, 
one at an irradiation time of 6 hr., and three at an irradiation time of 12 hr. 

The amount of Co* activity formed by successive neutron capture in Co®® 
is given in the following expression, for Ast > 1: 


Aaa _ | aig t (gx sats) a. a 
(1] oN; = o34| (e 1) ro “i % Se, ee aa 


° s r eeet 
+ (12 + 913) | + ox i (l-—e™') +o ; a | 


where the subscripts 1, 2, 3, 4 refer respectively to Co**, Co®* (10.5 min.), 
Co® (5.25 yr.), and Co™ (105 min.), and the other symbols have the following 
meaning: 
A—disintegration constant (sec.~!) 
t—time of irradiation (sec.) 
N—number of atoms present at time ¢ 
o@—thermal neutron flux (cm.~ sec.~!) 

o12. and o;;—refer to the effective neutron capture cross sections of Co*9 
(at the position of the ‘‘rabbit’’ in the NRX reactor) to form 
the isomers Co®* (10.5 min.) and Co® (5.25 yr.) respectively. 
These are derived from the corresponding thermal neutron 
cross sections, 032’ and o,;’, by the addition of terms to each 
which take into account the contribution to the activities by 
the capture of fast neutrons. 

From other data it is estimated that approximately 2.5% of 
the total Co® (5.25 yr.) activity produced (after decay of the 
10.5 min. activity) is due to capture of epicadmium neutrons. 
This leads to a value of 35.1 barns for (012 + 013) using the 
value 34.2 barns (14) for the corresponding thermal neutron 
capture cross sections (o12’ + 013’). It is however not known 
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how the difference of 0.9 barns is to be divided between o12’ and 
013. The procedure followed here has been to divide the 0.9 barns 
between a1.’ and o;;’ in the ratio of their magnitudes. Since 
the measured value of o12’ is 18.3 + 1.7 barns (7), any error 
introduced on account of incorrect assignment of these terms 
will be less than 3% in each of o12 and 043. The value of o;3’ is 
obtained from (o12’ + 013’) and o12’ by subtraction. 
o24 and o34—refer respectively to the pile neutron capture cross sections of 
Co®* (10.5 min.) and Co® (5.25 yr.) to form Co. (105 min.). 
They are the effective neutron capture cross sections at the 
position of the “‘rabbit’’ in the NRX reactor. At this position 
the ratio of fast to thermal neutrons is 0.589. 
Approximations made in the derivation of equation [1] are: 
(a) the number of Co®* atoms (NV) does not change during the irradiation, 
(b) 100% of the Co®* (10.5 min.) disintegrations lead to the Co® (5.25 yr.) 
state, 
(c) the rate of loss of both Co® isomers by neutron capture is small com- 
pared with the rate of loss of Co®* (10.5 min.) by decay, 
(d) the rate of loss of Co® (5.25 yr.) by decay is small compared with the 
rate of loss of Co®* (10.5 min.) by decay, 
(e) the rate of loss of Co™ (105 min.) by neutron capture is small compared 
with the rate of loss of Co® by decay. 
The following values have been employed for the constants in equation [1]: 
Ae = (1.10 + 0.02) & 1073 sec, 
hy = (1.10 + 0.066) X 10-4 sec—!, 
o12 = 18.8 + 1.8 barns, 
16.3 + 2.3 barns, 
35.1 + 1.4 barns. 


ll 


013 


F12 + O13 


The thermal neutron flux was determined by comparison of the sample 
with a calibrated Co® (5.25 yr.) source after the experiment. 

From equation [1] it can be seen that the yield, \4N4, differs, as a function 
of neutron flux and time of irradiation, from the yield which would be ex- 
pected from single neutron capture in an impurity. Firstly, for sufficiently 
large o34, the activity will not approach saturation with increasing time of 
irradiation, but will increase almost linearly (for \2¢ >> 1) owing to the term 
(o12 + o43)t. Secondly, the activity is proportional not to the first power but 
rather to the square of the neutron flux. Both of these characteristics have 
been utilized to demonstrate that the activity detected is due to Co®™ formed 
by successive capture of two neutrons in Co®’. Finally, the activities measured 
at the various irradiation times were inserted in equation [1] and values 


obtained for o24 and o34. 


III. RESULTS 
Two samples were irradiated, each for 3 hr., one at a neutron flux approxi- 
mately one-half that employed for the other. The results of these trial runs 
are shown in Figs. 1 and 2, and in Tables I and IT. In the tables, the estimated 
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Fic. 1. Decay curve obtained from Co®* sample irradiated for three hours at thermal 
neutron flux of 5.38 X 10 cm.~? sec.“ The background is due to Co® y-rays. The compo- 
nent of half-life-~~ 11 min. is due to B-rays from the Co®* (10.5 min.) isomer. Discriminator 
bias has been set to reject pulses from the Co® (5.25 yr.) B-rays. Where not denoted by vertical 
bars, the standard deviations are less than the radii of the circles. 
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Fic. 2. Decay curve obtained from Co®® sample irradiated for three hours at thermal 
neutron flux of 2.8 X 10'3 cm.~? sec.~' Comparison with Fig. 1 indicates that the yield of 
the 105 min. activity is proportional to the square of the neutron flux. Where not denoted by 
vertical bars, the standard deviations are less than the radii of the circles. 
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TABLE I 
VARIATION IN AMOUNT OF es MIN. ACTIVITY FORME D WITH CH ANGE 











| t | M Az N3 | NaN | [Aa s) /( (As Ns) }e| 2) (" =), 
Trial | (hr.) (mgm. )} (104 dis. /sec.) (dis. /sec.) awo/ (A3NV. (| N M 
2 | 3.00/ 0.215 | I8 46 ; 40.544.1 | 1.75 £0.29 1.90 + 0.07 
3 3.00 | 0.182 | 8.21 10.34-1.3 | 
TABLE II 


AMOUNT OF 105 MIN. ACTIVITY FORMED AS A FUNCTION OF TIME OF 
IRRADIATION. THE SYMBOLS ARE DEFINED IN THE TEXT. 


| (Ag.N4)/(G? N a 








| | M ANs Fs | 
Trial | thr.) (mgm.)| (10* dis. /sec.) )} (I 0v/cm. 2/sec.) | G0-* cm. ) | | Ww eighted mean 
1 | 1.53 | 0.318 | nm | 4.32 | 440.3 1440.3 
2 | 3.00] 0.215 18.46 | 5.30 | 1.70.3) | 
3 | 3.00 | 0.182 8.21 | 2.79 | 1940.37 | 1.7+0.2 
4 | 3.00] 0.156 | 12.16 | 4.82 | 115 £0.28 | 
| | | 
5 | 6.00} 0.285} 44.4 | 4.81 | 4841.1 | 4841.1 
| | | 
6 | 12.00 | 0.340 | 100.7 4.59 } 8.0+1.9) | 
7 | 12.00| 0.243 | 67.6 4.30 | 944lde | 8.5408 
8 | 12.00 | 0.332 93.0 4.33 | 7.8241.2) | 


error in the mass, M/, of Cc®® is 2%; in A3N3, 2%; and in the neutron flux 
¢, 6.5%. As shown in the last column of Table I, the ratio of the flux during 
trial 2 to that during trial 3 was 1.90 + 0.07. This was determined after the 
experiment by comparison of the two samples. The subscripts 2 and 3 refer 
to the trial number. If the activity 444 were due to absorption of neutrons 
in an impurity, or to any process depending on the first power of the flux, 
then the ratio (AsN4/N1)2/(AsN4/M1)3 should also be the ratio of the fluxes. 
In this event the ratio shown in column 6 of Table I would be unity. On the 
other hand, if the activity \4NV4 is due to Co, then (AqN4/N)2/(AaN4/M1)3 
will be the ratio of the squares of the fluxes, so that the quantity shown in 
column 6 of Table I should be the ratio of the fluxes, i.e. 1.90. 

The results for all trials are set forth in Table II, and a plot of the weighted 
mean values of (AyN4)/(¢?N1) X 10” a time of irradiation is shown in 
Fig. 3. The mean values of (A4N4)/(¢?N1), when inserted in equation [1], yield 
four equations in the two unknowns o24 ee o34. Solution of these equations, 
weighting each according to the error assigned to the values of (AyN4)/(¢?™j), 
yields the values: 

oo, = 95 + 45 barns, 
6.0 + 1.4 barns. 


O34 


These values were then inserted into equation [l] and the function 
(AsN4)/(¢?N1) X 10” vs. ¢ plotted as the solid line in Fig. 3. Although not a 
straight line, the curve agrees very closely with the straight line derived by 
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Fic. 3. (AgN4)/(¢2.Mi) X 10% as a function of time of irradiation. The experimental points 
are shown as circles. The vertical bars represent the probable errors. The solid curve is the 
computed function using the values o24 = 95 barns, o34 = 6.0 barns. The dashed curve illus- 
trates the growth of activity of 105 min. half-life to be expected from single neutron absorption 
or any process depending on the first power of the neutron flux. The curve has been normalized 
to the solid curve at t = 3 hr. 


a least squares fitting of the experimental points. The dashed curve illustrates, 
for comparison, the growth of activity of 105 min. half-life to be expected for 
single neutron absorption or any process depending on the first power of the 
neutron flux. The curve has been normalized to the solid curve at ¢ = 3 hr. 


IV. DISCUSSION 


Within the estimated errors, the results shown in Table I indicate that the 
activity detected is due to absorption of two successive neutrons in Co** to 
form Co*. The same conclusion is reached from a consideration of the results 
of Table II as plotted in Fig. 3. The results of Table I do not however preclude 
the possibility that some small part of the activity detected is due to activation 
of impurities with lifetimes comparable with that of Co*. Manufacturer's 
estimates of the sample impurities indicate that perhaps a few per cent of 
the activity after 3 hr. irradiation, and considerably less after 12 hr. irradi- 
ation, may be due to Mn**(2.6 hr.). Mn however appears to be the only im- 
purity with a lifetime comparable to Co which can conceivably cause any 
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difficulty. Cu also occurs as an impurity (0.0001%-0.003%), but the half-life 
of Cu® is 12.9 hr. and that of Cu® is 5 min., both markedly different from the 
105 min. half-life of the Co®™ activity. Also, the thermal neutron capture 
cross sections of both Cu® and Cu® are quite low (4.3 barns for the former and 
2.1 barns for the latter), so that the presence of Cu in the estimated amounts 
does not introduce appreciable error in the measurement of the intensity of 
the 105 min. activity. 

It will be noted that the fractional error in o24 is much larger than that in 
a34. This is due to the fact that on the whole the measured values depend on 
experiments done at relatively long irradiation times. From equation [1] it 
can be seen that for lengthy irradiations, the yield is sensitive to the value of 
o34 but insensitive to the value of oo4. Thus for an irradiation time of 12 hr. 
over 80% of the activity is due to capture in the 5.25 yr. isomer. For irradi- 
ations of several days capture in the 10.5 min. isomer may be neglected 
entirely with small error. 
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THE INTENSITIES OF ATOMIC AND MOLECULAR FEATURES 
IN THE AURORAL SPECTRUM! 


By W. PETRIE? AND R. SMALL? 


ABSTRACT 
A considerable number of auroral spectrograms (dispersion 85-28 A/mm.) 
have been analyzed with a recording microphotometer, and integrated intensities 
obtained for the main features. Relative population rates of the relevant levels 
are determined from the measured intensities. Excitation mechanisms are dis- 
cussed and the absence from the spectra of certain low level atomic lines does 
not appear to be anomalous. It seems that atomic lines are emitted as the results 
of several processes, while the ionic lines are produced during simultaneous 
ionization and excitation through electron exchange. Vibrational temperatures 
of First Negative, First Positive, and Second Positive nitrogen bands are 
computed. 
INTRODUCTION 
A considerable number of moderate dispersion spectra of the aurora have 
been secured at Saskatoon, and all features of appreciable intensity identified. 
The majority of the spectra were from arc and band displays, and the greater 
part of the radiation originated from near the 100 km. level. A report on these 
spectra has already been published (18). In this paper we compare measured 
and calculated intensities of the permitted atomic lines and a number of 
nitrogen bands. 
MEASURED INTENSITIES 
The spectra were secured with an instrument which employed a step-slit 
of 11 steps, ranging in width from 0.02 to 2.0 mm. This type of slit has certain 
advantages for several spectroscopic problems. Both weak and strong features 
characteristic curve’’ of the emulsion 


“a 


may be measured on the same plate. The 
and its variation with wave length may be determined by microphotometering 
selected lines in the spectrum. The wave length variation of the emulsion 
sensitivity, and radiation losses in the spectrograph due to absorption and 
scattering, may be evaluated by measuring the intensity distribution in the 
spectrum of a standard lamp, and comparing with the theoretical distribution 
as given by Planck’s law. 

It was not feasible to superimpose the standard lamp spectrum on the 
auroral spectra owing to the small size of the plates used by the Schmidt 
camera; hence the lamp spectrum was recorded on separate plates. This 
procedure is satisfactory when using a step-slit, since “characteristic curves”’ 
may be determined for each individual plate. The auroral and standard lamp 
spectra were analyzed with a Jarrell-Ash recording microphotometer. The 
line and band profiles were converted from a logarithmic to a direct intensity 
scale, the areas under the profiles measured with a planimeter, and corrections 
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made for the varying wave length sensitivity of the emulsion and losses in the 
spectrograph. Finally, a correction must be made to remove the effects of 
atmospheric extinction; this point has been discussed in an earlier paper (17). 
In view of the uncertainty of this latter correction, it is well to check on the 
final intensities. This may be realized as follows. The relative intensities of 
bands in a v”’ progression are independent of excitation conditions, and may 
be computed if the relevant transition probabilities are available. Hence, for 
such a progression computed and measured intensities should agree. Con- 
sidering the uncertainties in the A values from which the computed intensities 
are obtained, the tables given later in the paper show that the agreement is 
reasonably good. 


CALCULATED INTENSITIES OF ATOMIC FEATURES 

In order to compute the intensity of a spectral feature it is necessary to 
know the Einstein A value or some related quantity. Information on the A 
values of atomic oxygen and nitrogen transitions is now available, as the 
result of recent calculations. Since the expressions which must be evaluated 
include the radial wave-functions, it is necessary to have acceptable functions 
for all the states of the atoms in which one is interested. One approach is to 
adopt Slater wave-functions, and these have been used by Petrie (16) to 
calculate OI intensities. More accurate wave-functions have been derived by 
Bates and Damgaard (3), and these authors compute data from which A 
values may be obtained. Using their notation, A = 2.02 & 10!° s/god* sec.“!, 
where s is the “‘line strength’, go the statistical weight of the upper level of 
the transition, and \ the wave length of the spectral line. The “‘line strength” 
is the product of the factors aSo*s/Ls. The meaning and method of obtaining 
these factors have been discussed in an earlier paper (14). In this paper we 
use the o? values as obtained from the tables given by Bates and Damgaard (3). 
Although these are as accurate as can be obtained at present, it must be realized 
that there are uncertainties in the values, when the arguments which follow 
later in the discussion are being considered. The term values necessary for the 
use of the tables were taken from ‘‘Atomic Energy Levels’? by Moore (12). 
In the case of a few transitions the tables had to be extended by extrapolation 
in order to obtain a”. The accuracy of the resulting A values is very questionable, 
especially for the transitions to the ground configuration. In this paper we discuss 
multiplet rather than line intensities, hence we have summed over all lines of 
each multiplet in the case of both measured and calculated intensities. 

When dealing with a source in thermal equilibrium, the populations of the 
various levels are given by the Boltzmann law, and the intensity of a spectral 
line may be expressed as 
(1] T = N,Ann hay = No mt 78°"? Any hinn’ 

Wo 
where 
N, = the number of atoms in the upper level involved in the transition, 


No = the number of atoms in the ground level of the atom, 
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w, = the statistical weight of level n, 
wo = the statistical weight of the ground level, 
E,, = the excitation potential of level n, 


T = the temperature of the source. 


If the source is not in thermal equilibrium, we cannot say what are the relative 
populations of the various levels, but must instead consider the rate of exci- 
tation to these levels. One approach to the problem is to compute the relative 
intensities of lines assuming that all the relevant levels are populated at the 
same rate. Since a “‘cascade”’ process will follow the excitation of an atom to 
a certain level, the rate of population will be the sum of the 
rates due to the direct excitation and ‘‘cascade’’ processes. As 
an example, suppose that N atoms enter the levels A, B, C, D, 


D 
and £E per second, from the ground level. If transitions ED, 
C EC, DB, and CB are allowed, then the population rate of level 
B is 
B 
Acs ( Arc ) App ( Arp ) 
N a eee N . —— 7 “+ ences N anata eens N . 
LA Cx r LA EX LA DX . DA ex 


The intensity of the line arising from the transition BA is 


A A A A A 
2) I= |v Ace (yy — v) os w a v) | 24 hvpa. 
2] *¥AN "Lia ‘tae eae 


It is also possible to proceed in the reverse manner, that is to determine the 
population rates of the levels from which the observed lines arise, from the 
relative intensities of these lines. This procedure has been followed here, and 
each set of calculations will be discussed separately. 


(a) OI 

The highest level line found in the auroral spectrum (18) is from the transi- 
tion 2p°3p(°P) — 2p*6d(°D°). All levels between °D° and the ground level, as 
well as the permitted transitions, are shown in Fig. 1. This scheme assumes that 
all levels between the ground and *D° levels are populated by the excitation 
mechanism. The transitions giving lines in the spectral range mentioned, the 
computed A values to the nearest whole number, measured intensities, and 
calculated relative population rates are given in Table I. These latter data are 
determined from the equation which expresses the intensity of a line, that is, 


,_A 
[3] T=Nyyl 


where JN is the total population rate of the levels concerned, from direct 
excitation and cascade processes. In this and following tables only one line 
of each multiplet has been listed, whereas the intensities refer to the complete 
multiplets. The multiplets which we believe are present in the auroral spectrum 
are those assigned a measured intensity. 
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TABLE I 
OI INTENSITIES 
| Excitation | Relative 
Transition Wave length (A) | potential | A value | Measured | population 

| (e.v.) | intensity rate 
= Se ee b 
2p'3s(5S)-2p93p5P) | 7774 | 10.69 | 3 X 107 100 1.0 
2p°3p(3.S°)-2p°3p(3P) | 8446 | 10.94 3 X 107 150 1.6 
2p33s(5S°)-2p%4p(5P) | 3947 12.23 | 4x 105 
2p835(3S%)-2p4p(8P) | 4368 | 12.31 | 9X 108 25 0.23 
2p°3p(2P)-2p33s(3D°) | 7995 | 12.49 2 xX 107 
2p'3p(5P)-2p85s(8S) | 6456* | 12.61 | 1X 107 
2p'3p(3P)-2pr5s(2S%) | 7254* 12.64 | 6X 108 
2p'3p(3P)-2p4d(3D°) | 7002 | 12.70 | 3X 108 
2p33p(5P)—-2p34d(5D°) | 6156 | E278 | 8 >< 00° 12 0.18 
2p33s(85°)-2p35p(sP) — | 3692 | 12,82 4 X 10° 
2p°3p(5P)—-2p%6s(8S°) | 5436 12.96 | 4X 10° 10 0.07 
2p'3p(3P)-2p¥s(8S%) | 6046* | 12.98 | 3x 108 
2p33p(P)-2p35d(8D°) | 5958* 13.01 | 8 X 105 
2p*3p(5P)-2p%5d(§D°) | 5330 13.01 | 3 X 108 8 0.08 
2pi3p(5P)—2pi7s(8S%) | 5020 | 13.15 | 2X 108 
2p33p(8P)-2p7s(8S%) | 5555* | 13.16 2X 105 | 
2p'3p(5P)-2p%6d(5D°) 4968 13.18 | 2X 108 | 7 | 0.06 


* Obscured by other spectral features. 


It has long been considered anomalous that certain low level OI lines appear 
to be missing from the auroral spectrum. Of the multiplets listed in Table I 
which are not in the spectrum, only those with wave lengths 3947, 7995, 7002, 
3692, and 5020A are in regions unobscured by other features. Of these \\ 3947, 
3692 have relatively small A values. In addition, by considering the transitions 
illustrated in Fig. 1 and the A values associated with them, it can be seen 
that the levels from which the missing lines originate will be populated rela- 
tively slowly from the cascade process. Hence the absence of these lines is to 
be expected. It should be mentioned that Meinel (10) has observed part of 
the 7995A multiplet in his spectra. We should also note the anomalous in- 
tensities of the lines 7774 and 8446A. The intensity ratio of these lines varies 
considerably from one spectrum to another. The latter is always the stronger 
of the two, and the intensity ratio may be as great as 3:1. It would be in- 
structive to compare the intensity ratio of these lines and the intensity of one 
of the Balmer lines, but the authors have not yet been able to do this. 

We will now consider the manner in which the permitted oxygen levels can 
be excited. Oxygen atoms in the ground *P levels will experience collisions 
with incoming protons, and undergo transitions to higher triplet levels. It is 
important to realize that a change of electron spin is very unlikely during the 
collision with a proton (9). Hence, if the excitation resulted solely from this 
process, only the triplet levels would be excited, but such is not the case. A 
few oxygen atoms, relative to the total number, will at any given time be in 
the excited metastable levels 1D and !S, and may be raised to higher singlet 
levels by collisions with protons. However, this process will be of negligible 
importance. The hydrogen atoms which result from the capture of electrons 
by protons may by collision with oxygen atoms in the ground *P levels excite 


ee 


aera 
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Fic. 1. OI levels and transitions. 


both singlet, triplet, and quintet levels. We can only give a crude estimate as 
to how effective this process will be. Suppose at the 100 km. level we have an 
incoming hydrogen atom density of 0.1 per cm.*; a velocity of 10° cm./sec.; an 
oxygen atom density of 10"/cm.*; and a collisional cross section for the exci- 
tation of permitted oxygen levels of 107'§ cm.? Then the number of exci- 
tations/cm.®/sec. is ~ 10°. Oxygen atoms will also be excited by electron 
collisions. If electrons accompany the protons and have the same velocities, 


_ the former particles will be stopped in the atmosphere above the region from 


which the major part of the auroral light is emitted. It has been pointed out 
elsewhere (15) that if a proton is to penetrate to the 100 km. level, its energy 
and velocity of entry must be of the order of 3 Mev. and 2.4 X 10° cm./sec. 
An electron of the same velocity has energy of 1640 e.v. The range of such an 
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electron will be much less than that of the proton and the electron will be 
stopped near the 200 km. level (5). However, even if primary electrons cannot 
reach the main auroral region, fast protons produce electrons of a wide range 
of velocity, and these particles will be able to excite a number of oxygen levels. 
Bates (4) has shown that for any given spectral series the excitation cross 
section falls off rapidly with increase of the principal quantum number. For 
both electron and proton excitation the cross section is approximately pro- 
portional to o” (2). We would then expect the population rate of the higher 
levels to be small, and reference to Table I will indicate that this is indeed the 
case. A noticeable feature of the OI data is the presence of lines from several 
high d levels and the absence of lines from higher s levels. The identifications 
cannot of course be considered final until the multiplet structure is resolved, 
and present equipment is not able to do this. At any rate, the lower level lines 
which are certainly due to OI are likely excited by the several collisional 
processes discussed above. 


(b) OIT 

Data on the low level OII transitions are given in Fig. 2 and Table II. The 
line at 6641A, although near the (6,3) First Positive band, appears to be 
absent. The line at 3973A also appears to be absent. Superimposed on Hy are 
two features which are likely the lines AA 4347, 4351 of the (?D-?D°) multiplet. 
The line at 4072A, although weak, seems to add to the emission of the VK 
(2,13) band. The line at 4415A may be anomalously intense owing to an over- 
lapping (3,8) Second Positive band. 


34(*F) 
3p(?D°) 
3p (?F°) 
3s(?S) 
3p(? P*) 
2s2p*(?P) 
3p(*S°) 
3p(? D°) 
3p (* P°) 
3s(?D) 
3p (* D°) 
3p(?S°) 
2s2p*(*S) 
3s(? P) 
3s(°P) 
2s2p*(?D) 
~ 2s2p*(*P) 
2p3(?P*) 
2p5(?D°) 
2p3(4S°) 
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TABLE II 
OIL INTENSITIES 








| 
| 


| Excitation 





| 
| Relative 
Transition | Wave length (A)| potential | A value | Measured population 
(e.v.) | | intensity rate 
— cease sa societies co a | | | ad 
2p'3s(2P)-2p°3pPS) | 6641 | 25.18 | 3x 107 | 
2p3s(4P)-2p°3p(4D) | 4639 | 25.55 | 9x107|/ 200 | 1.0 
2p3s(*P)-2p3p4P%) | 4317 | 25.74 | 1108] 175 0.8 
2p°3s(2P)-2p°3p(2D") 4415 | 26.14 |} 1108 | 250 1.6 
2p93s(4P)-2p"3p(45%) | 3713 | 26.19 | 2x10°| 6 | O85 
2p'3s(2P)-2p°3p2P%) | 3973 | 26.45 | 1X 108 | | 
2p°3s(2D)-2p°3pCP) | 4591 | 98.24 |9x107| 80 | 04 
2p°3s(2D)-2p°3p(2D°) | 4347* 28.39 1X 108 | 
28.58 | 2X 108 | | 


2p°3p(4D°)-2p°3d(*F) 4072* 58 | 


* Obscured by other spectral features. 

In view of these uncertainties, little can be said about the interpretation of 
the OII intensities. However, it is apparent that the excitation process or 
processes will differ from those operating in the case of OI. There are relatively 
few O* particles, and hence the direct excitation of the ion will not be im- 
portant. An alternative process. is that oxygen atoms are simultaneously 
ionized and excited. This can only effectively occur through the process of 
the removal of one electron and the exchange of a second with the incident 
electron. By considering the possible spin changes during these processes, it 
can be seen that both doublet and quartet levels of OII will be excited. 
The probability of a proton capturing one electron from the atom and exciting 
a second is very small. The capture process is most effective at low proton 
energies (~ 2800 e.v.) and it requires much higher energy hydrogen atoms 
to be effective in inelastic collisions. In view of these remarks, it seems likely 
that OII permitted levels are excited as oxygen atoms are simultaneously 
ionized and excited through the process of electron exchange. 
(c) NI 

Data on the relevant NI transitions are given in Fig. 3 and Table III. 
There are few observed NI lines of appreciable intensity. 
TABLE III 


NI INTENSITIES 








| Excitation | Relative 
Transition Wave length (A)| potential | A value | Measured | population 
(e.v.) intensity | rate 
a eae ee eee ae en a i a eS ee 
2p°3s(4P)-2p°3p(tD) | 8680 11.70 | 3X10? | 100 =| 1.0 
2p?3s(4P)-2p’3p(4P°) 8186 | TE.7e | -35¢107 30 | 0.29 
2p?3s(*P)-2p°3p(4S°) | 7468* | 11.94 | 4X10? 
2p°3s(?P)-2p'3peP) | 8629 | 12.07 3 X 107 
2p3s(?2P)—-2p*4p(2S°) | 4915 | 13.14 4X 10° 25 0.22 
2p*3s(4P)-2p?4p(4D°) 4253* | 13.19 2 X 10° | 
2s2p'(4P)-2p24 p(*D) 5329 | 43.19 
2s2p*(*P)-2p*4 p(4P°) 5281 |. eae 
2p?3s(*P)-2p*4p(+P) | 4223 | oe 7X 108 | 
2p*3s(*P)-2p*4p(4S°) 4151 | 13.26 | 3X 10° | 


* Obscured by other spectral features. 
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4p(*0°) 
4p(*S°) 
3d (7D) 
3d (*D) 
3d (*P) 
3d (?F) 
3d (*F) 
3d (?P) 
4s(*P) 
4s (*P) 
3s(?D) 
3p(*P*) 
3p(70°) 
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3p(*S?) 
2s2p*(*P) 
3s(*P) 
3s(*P) 
2p*(?P*) 
2p?(? D°) 
2p%*S°) 


The excitation process of NI and OI levels may differ. Atomic nitrogen lines 
may be produced during collisions with the molecule, and since the concen- 
tration of NI in the atmosphere is unknown, we cannot say whether the major 
part of the atomic radiation is produced by collisions with nitrogen atoms or 
nitrogen molecules. Investigations in progress at the Queen’s University of 
Belfast (2) indicate that certain atomic nitrogen lines are readily produced 
during collisions with the molecules. The writers have carried out an investi- 
gation of the spectra of low pressure air discharges at the University of Saskat- 
chewan. The four low level infrared multiplets given in Table III appeared 
with considerable intensity in the spectra of uncondensed, condensed, and 
electrodeless discharges through air. The available evidence (8) suggests that 
a nitrogen molecule dissociates into a neutral and ionized nitrogen atom, 


rather than two neutral nitrogen atoms. The reaction is 


Ne +e—N + Nt + 2e. 


Hence, neutral and ionized atomic lines may be produced in a collision with 
the molecule. On the other hand, it is almost certain that nitrogen atoms are 
constituents of the normal night atmosphere, and atomic lines will be pro- 
duced during collisions with these particles. The arguments presented re- 
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garding the excitation of oxygen atoms will also apply to the nitrogen atoms. 
It appears then that the production of nitrogen lines is a complex problem. 


(d) NII 
Data on the relevant NII transitions are given in Fig. 4 and Table IV. The 
line at 3995A overlaps the (1,4) Second Positive band, but from the appearance 
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Fic. 4. NII levels and transitions. 





TABLE IV 
NII INTENSITIES 


_ | Excitation Relative 
Transition Wave length (A)| potential | A value | Measured | population 
(e.v.) intensity | rate 
2p3s('P°)-2p3p('P) 6482 20.32 3 X 107 100 1.0 
252p5(!D°)-2p3p(P) 1895 20.32 30 0.68 
2p3s(°P°)-2p3p(3D) 5667 20.58 5 X 10 180 1.3 
2p3s(2P°)-2p3p(3S) 5011 20.85 8 X 107 20 1.6 
2p3s(@P")-2p3pP) 4601 |} 21.07 1 xX 108 50 0.35 
2p3s(1P°)-2p3p('D) — | 3995* | 21.51 | 1X 108 
2p3s('P%)-2p3pUs) 3437 | 22.01 | 2X 108 
2p3p(@D)-2p3d(3 F°) 5001 23.04 1 X 108 150 0.87 


_ Obscured by other spectral features. 
of the blended feature the line must be weak. The line at 3437A is in the 
region of the (1,10) VK band, but may be seen although it is very weak. 

Arguments regarding the excitation of NII levels parallel those presented 
in the discussion of OII intensities, and if nitrogen atoms are involved in the 
excitation, the most likely process is the simultaneous ionization and excitation 
of nitrogen atoms through electron exchange. In addition, as we have pointed 
out, collisions with nitrogen molecules may produce excited N* ions. 
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As this discussion has indicated, an analysis of the intensities of oxygen 
and nitrogen lines in the auroral spectrum does not lead to an unambiguous 
explanation of the excitation processes operating. It is apparent, as has long 
been realized, that the excitation is complex in nature. It appears that the 
ionic lines are produced largely as the result of inelastic collisions between 
the neutral atoms and electrons whereas the atomic lines are emitted as the 
result of several processes. It is still desirable that spectra of sufficient dis- 
persion to resolve multiplet structures be secured, since only then can the 
identifications of some of the atomic features be considered certain. 


Unidentified Features 
There are a number of weak, but clearly seen sharp features, which appear 
to be atomic in origin. The wave lengths of these features are given in Table V. 
TABLE V 


UNIDENTIFIED FEATURES OF POSSIBLE ATOMIC ORIGIN 
Wave length (A) 





4042 | 4358 | 4950 6137 

4047 | 4684 | 4958 6141 

4180 | 4687 | 5916 | 6402 
| | 


6414 


A number of these features have been measured by Vegard (22) and attributed 
to NII and OI! multiplets and NO bands. Nicolet (13) has drawn attention 
to coincidences between the wave lengths of a number of auroral features and 
NO bands. On the whole, the identifications are not too convincing for a 
number of reasons, and further discussion at this time does not seem profitable. 
The features in the red region of the spectrum may be portions of OH airglow 
bands. 


Nitrogen Bands in the Auroral Spectrum 

(a) Negative Group Bands 

Relevant data on these bands are given in Table VI. The transitions listed 
here are those with the greatest theoretical strengths (21). The ones marked 
with an A are definitely present in the auroral spectrum, those denoted by P 
are probably present, and the remainder seem to be absent or obscured. The 
same notation is used with the tables which describe the other band systems. 
All the bands in Table VI have been observed in laboratory discharges (6). 
Notes regarding a number of these features and others which might be ex- 
pected to appear are given as footnotes to the table. As these notes indicate, 
there are relatively few of the stronger bands which are unobscured and suita- 
ble for intensity measurement. Those which have been measured and their 
relative intensities are given in Table VII. It is to be noted that these in- 
tensities fall to smaller values with increasing v’’ than do those of Vegard and 
Kvifte (22). The same result was found in an earlier discussion of the intensities 
of the Second Positive bands (17). This is at least partly due to the fact that 
we have measured spectra of higher resolution, and have given a different 
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TABLE VI 
NEGATIVE GROUP BANDS 


























~ v”’ 
os Ss 1 | 2 | 3 4 5 6 7 8 9 
0 | 3014 | 4278A | 4700A | 5228 | 
1 | 3582A"| 3884A | 4236A | 4652A | 5149A | 5754P2 
2 | 3308P?| 3564A | 3858A*/ 4199P* | 4600A | 5076P?| 5653 
3 3299P* | 3549A | | 4167 | 4554A8) 50127 | 55648 
4 | | 3203P* | 3538P> | 4141P3| 4516A | 4958A | 5486° 
5 | | | | | 3533P? | 4121 | 4486+) 4018% | 5421P" 
| | { 








1 The (83,12) VK band has the same wave length, but since it is weak, the N.G. band will not 
be appreciably affected. 

2 This band is in the region of a 1 P.G. band. 

3 This band is in the region of a 2 P.G. band. 

4 This band is likely blended with a P.G. (4,7). 

5 There is a weak feature in the spectrum centered at \8296 which is likely a blend of the 
(3,1) and (4,2) bands, although Barbier and Williams (J. Geophys. Research, 55: 401. 1950) 
attribute this to an airglow band. 

6 This band may be blended with some other feature since it seems abnormally strong. 

7 This band is in the region of an NII line. 

8 This band ts in the region of O.* (2,1). 

9 There is a feature near this wave length but it appears to be broad and diffuse. 

10 This band is in the region of an NI line. 

1 There is a feature near this wave length but it appears too broad. 














TABLE VII 
MEASURED AND CALCULATED INTENSITIES OF N.G. BANDS 
ae | Pee ie a 
vy” | 0 1 2 3 4 | 5 6 | 7 
aN — _ - _ nr aa loan —j———— —— 
z j | | | | | | | 
[Ml eC ruMEE u|c|M|c|M|c|m|c a Lee 
sonics | i seinem eit = | | | = anne 
0 | 100/100] 30 | 22} 6 | 3.3] 2 lo.9| 2 a a 
1 10} 62| 7} 37| 6 | 2.6) 2 |4.6/0.8/1.7| 
2 2/79 /0.4|10.5 = 10.2} 2.1| | 
3 11.0 | | 0.6| | 
4 | | | | | | 0.4} 0.2) 
5 | | | | | | 0.2) 
| 


Oho Ps ed 





evaluation of atmospheric extinction. The intensities of several of the weaker 
“ 


bands cannot be considered reliable. 

Bates (1) has discussed in detail the information which may be derived 
from band intensity data, and we will apply his methods to the intensities 
given here. The band intensities are computed from the theoretical py» 
values and the equation 


i) Ign = Sebeeroer 


where g,’ is the population rate of the vibrational level of the upper electronic 
state. These intensities are given in Table VII for the case g,, = 1, i.e. for 
equal population rates of the excited levels. The theoretical p,»» values of 
Turner and Nicholls (21) have been used. By comparing the observed and 
calculated intensities, it is possible to compute the relative population rates 
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of the levels of the upper electronic state. The rates for the levels v’ = 0, 1, 2, 
are 1.00, 0.15, 0.020. Bates (1) has determined theoretically the population 
rates which would result from several different excitation processes. By ex- 
amining his tables, it appears that if the band system is emitted as the result 
of simultaneous ionization and excitation of the neutral molecule by electrons 
the vibrational temperature must be near 1000° K. Although other views have 
been expressed (19, 20), the band intensities favor the combined ionization 
and excitation process, the high vibrational temperatures perhaps indicating 
that protons are involved. 


(b) Second Positive Group Bands 
Data for the stronger bands of this system, including the bands present 
and likely present in the auroral spectrum, are given in Table VIII. 


TABLE VIII 
SECOND PosITIVE GROUP BANDS 








wt 1 at es a, eS 
41 ed ew) ae) a ee aes 
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0 | 3371A | 3577A | 3805A | 4059A | 4344P"| | | 




















1 3159A | 3339A | 35387A | 3755A | 3998A | 4269? | 4574A | 4917P3| 
far | 3136A | 3309P? | 3500P* | 3710A | 3943A | 4201P?| 4490P?| 4815° 
oH | 3117A | 3285A | 3469P8 | 3672P4 | 3894P? | 4142P? | 4416P7 
eo | 3267P4 | —_ 3642 | 3857P?| 4094P*| 4356P8 
| | 

! This band is in the region of Hy and is not well resolved. 

2 This band is in the region of a N.G. band. 

8 This band may contribute to the intensity of the NI line at \4915. 

4 This band ts in the region of a VK band. 

5 This band has the wave length of a weak auroral feature. 

6 This band is in the region of the NI line at \3466. 

7 This band is in the region of an OI line. 

8 There is an auroral feature at the wave length of this band but it appears to be anomalously 

Strong. 


The measured intensities of the unobscured bands are given in Table IX. 
The theoretical intensities which result from equal population rates of the 
excited levels are calculated from [4] and are given in Table 1X. A comparison 


TABLE IX 
MEASURED AND CALCULATED INTENSITIES OF SECOND POSITIVE GROUP BANDS 


0 | 1 2 | 3 | 4 





5 | 6 
| l | | | al re ee 
M|/C|M|C/|M/C;}|M|C|M|C|M/C/|M]|C 
cco eae nem nore an mre preci eam 
0 | 100] 100} 59 | 40 | 21 | 14] 8 | 4.5] eS a 
1 | | 4] @ | 201 ai wet 8 1 8 | 3 
2 i } | 3]91] 3 | 
3 | 3 4 | | | 
| | | 





of measured and calculated intensities indicates that the relative population 
rates of the levels v’ = 0,1 are 1, 0.45. These values are much the same as 
those found from lower dispersion prism spectra (17) which extended further 
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to the ultraviolet. If we compare only those bands which have comparable 
wave lengths in order to minimize errors due to atmospheric extinction, 
somewhat higher relative population rates are found. Hunten (7) has measured 
the intensities of a number of the bands using a photomultiplier scanning 
spectrograph, and has obtained still higher population rates. In the earlier 
paper already referred to the conclusion was reached that the vibrational 
temperature is < 500° K. However, the relative population rates and hence 
the temperature are very sensitive to changes in the relative intensities of the 
bands, and the conclusion seems to be that the data are not sufficiently re- 
liable or extensive to yield accurate vibrational temperatures. Nevertheless, 
the intensities are of the order expected if the excitation is by collision from 
the ground level. It should be noted that excitation from the ground level 
can be achieved by electrons but not by protons and hence a lower vibra- 
tional temperature is not unexpected. 


(c) First Positive Group Bands 


Data for this system including the bands present and likely present in the 
auroral spectrum are given in Table X. 











TABLE X 
First PosiItIVE GROUP BANDS 
a | Loe ft ee -% ae 
Ar ee ee ee 
0 | 
1 | 8912A | 
2 | 7753A | 8722A | 
3 | 6875A | 7626A | 8542A | 
4 | 6185A' | 6789A | 7504A | 8369 | 
B'| | 6127A | 6705A | 7387A | 8204 
6 | | 6070A | 6624A | 7273A | 8047P? | 
oI | 6O14A | 6545A | 7165A | 7896P? 
8 | | 5959A | 6469A | 7059A 
9 | | 5906A | 6395A | 6968A 
10 | 5854A | 6323A 
11 | | 5804A | 6253A 
| | | { 





| 
1 ‘This band is in the region of 1 P.G. (12,9). 
2 This band is in the region of a Meinel N2* band. 

A considerable number of these bands are partially obscured by other 
features, including OH bands which appear on the longer exposure spectro- 
grams. The intensities of a number of bands are given in Table XI. Certain 
of these are partially obscured, but we have corrected for the blending effect 
of other features. Data on py» values for this system have been given by 
Bates (1) and Montgomery and Nicholls (11). Using Nicholl’s theoretical 
Py» values, intensities are calculated from [4] for equal population rates of 
the vibrational levels of the upper electronic state. The results are given in 
Table XI. The relative population rates of the levels v’ = 1, 2, 3, 4, 5, 6, 7, 8 
appear to be of the order of 1.0, 0.58, 0.50, 0.34, 0.30, 0.11, 0.13, 0.030. 
Because of the uncertainties in the p,-» values, no great weight can be given 
to these rates. As Bates (1) has pointed out, excitation of the nitrogen molecule 
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TABLE XI 


MEASURED AND CALCULATED INTENSITIES OF First PosITIVE GROUP BANDS 


| 
@ || 



































| oO  .. & Poe 4 5 6 7 8 9 

\ l Pa ce, = = l Pea 
v\|M;}C|M/C|M|C|M/C;/M/|C/M|C/M/C|M|C/M|CjMI[C 
0 

1 | 100] 100 | 

2| 35! 60 

3 | 3/ 6 | 

4 12 | 43] 29 | 73 

5 0.9| 2| 12 |67| 15 | 48 

6 11] 8| 9/93 | 

7 | 11.3117] 6 |102/1.8] 15 

8 | = 1.3| 38} 3 | 109 

9 Pa 1.1 1.8 1.3 
10 1.1 
ll 

















| 0.9 


| { 








from the ground level would populate the upper vibrational levels at a high 
rate. 
(d) Vegard-Kaplan Bands 


Data for this system including the bands present and likely present in the 
auroral spectrum are given in Table XII. 


TABLE XII 


VEGARD-KAPLAN BANDS 








11 12 13 14 15 17 
wae S Jan —— — am — - —_ 
| 3352A! 3602A 3887 P? 
| 3198A 3425A 3682A 3979A 
3503A 3767A 4072A 4837A 


ARPwWHOes | 


| 
| 3855P? | 4171A | 4535A 
4772A 


1 This band can be clearly seen near the violet edge of the (0,0) 2 P.G. band. 
2 This band is in the region of N.G. (0,0). 
3 This band likely contributes to the (2,2) N.G. band. 


(e) Other Band Systems 


Vegard and Kvifte (22) attribute a number of features to Goldstein-Kaplan 
Nz bands and to NO bands. We conclude from an examination of our spectra 
that few if any of these bands are present. There are several coincidences, but 
the agreement is not convincing. Weak features of uncertain origin which have 
the appearance of bands are listed in Table XIII. 





| 
| 








| 
i 
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TABLE XIII 
WEAK AURORAL BANDS 








Wave length (A) | Possible identification Notes 
3484 Appears to be degraded to red 
3614 N.G. (10,9) Appears to be degraded « red 
3632 Broad feature, appears to be degraded to red 
3659 Appears to be degraded to red 
4378 
4424 VK (2,14) 
4433 
4467 N.G. (6,8) 
4725 2 P.G. (3,9) 
4750 | Broad feature, appears to be degraded to red 
4805 | Broad feature, appears to be degraded to red 
4832 | Broad feature, appears to be degraded to red 
4882 N.G. (6,9) 
4977 2 P.G. (4,11) 
4991 
5087 | Appears to be degraded to red 
5163 | Broad feature 
5356 Broad feature 
5373 1 P.G. (12,7) Broad feature 
| N.G. (6,10) 
5405 1 P.G. (11,6) Broad feature 
5422 N.G. (5,9) | Broad feature 
5476 1 P.G. (9,4) | Broad feature 
5488 N.G. (4,8) | Broad feature 
5514 1 P.G. (8,3) Broad feature 
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THE NEUTRON CAPTURE 7-RAYS FROM POTASSIUM! 


By G. A. BARTHOLOMEW AND B. B. KINSEY 


‘ABSTRACT 


The capture y-rays from potassium have been re-examined with greater resolu- 
tion than was used in previous experiments. The upper end of the spectrum 
has been carefully studied both with a sample of natural potassium carbonate 
and with another in which the potassium was enriched in K*°. From a comparison 
of the spectra two y-rays with energies of 9.39 + 0.06 and 8.45 + 0.02 Mev. are 
assigned to capture by that isotope. The strong y-ray at 7.757 + 0.008 Mev. 
previously ascribed to the ground state transition in K*° is now found to represent 
a transition to a low-lying excited state in that nucleus. 


INTRODUCTION 


Measurements of the neutron capture y-rays from potassium obtained with 
the aid of a pair spectrometer have been described elsewhere (4). The results 
obtained in that investigation (line width 170 kev.) may be summarized as 
follows: 

(1) a strong y-ray, B, with an energy of 7.77 + 0.03 Mev. and an intensity 
of 11 photons per 100 captures, which was ascribed to the direct transition 
to the ground state in K*°; 

(2) a weak y-ray, C, which was ascribed to the ground state transition 
in K®; 

(3) a complex and partially resolved spectrum extending from 6 Mev. 
downwards; 

(4) several weak y-rays with energies above that of B which were accounted 
for by capture in the rare isotope K*. 

It was pointed out that the intensity of the y-ray B is surprisingly high if 
this y-ray is, in fact, emitted in the ground state transition in K*° for it must 
then be of the magnetic quadrupole (72) or electric octupole (£3) type 
according to whether the spin of the compound nucleus formed by capture 
is | or 2 units. Such radiations are expected to be highly forbidden. Also if 
the weak y-rays with energies above that of B are actually produced by 
capture in K* then that nucleus must have a high capture cross section. 
Since the publication of these results, Pomerance (5) has shown that the 
capture cross section of K*° is 66 barns and contributes about 0.5% to the 
total potassium cross section. 

NEW MEASUREMENTS 

The present report contains the results of further measurements of the 
potassium capture y-rays made with better resolution (line width 100 kev.). 
The coincidence spectrum obtained with a sample consisting of about 1 kgm. 
of anhydrous potassium carbonate is shown in Fig. 1. It will be seen that the 
main features of the spectrum previously obtained are reproduced. The ener- 


1 Manuscript received May 15, 1953. 
Contribution from Atomic Energy of Canada Limited, Chalk River, Ontario. Issued as 


A.E.C.L. No. 66. 
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TABLE I 
ENERGIES AND INTENSITIES OF THE POTASSIUM y-RAYS 
| 
| Intensity in photons 
y-Ray Energy in Mev. | per 100 capturesin | Origin 
| natural potassium | 
isi ie — | acai — — 
Af 9.39 +0.06 0.02 kK" 
A’ 8.45 +0.02 0.1 Kl 
B 7.757 + 0.008 3.5 K+o 
c 7.34 +0.02 | 0.1 K# 
D | 6.994 + 0.007 | 1.3 K+ 
E | 6.31 +0.06 | 0.3 | 
F | 5.740 + 0.012 | 6 K* 
F’ 5.66 +0.02 | 4 K+0 
F” 5.50 +0.02 | 2.5 K+ 
G 5.38 +0.03 | 6 K+ 
G’ 5.18 +0.02 | 2 | K+ 
H 5.06 +0.02 | 3 | K« 
I 4.39 +0.03 | 4 K" 








gies and intensities of the strong y-rays are listed in Table I. The intensities 
were calculated from the original measurements, in which the peak counting 
rate of the potassium peak F was compared with that of the 9.0 Mev. y-ray 
of nickel emitted from a sample consisting of an intimate mixture of NixO; and 
anhydrous K.CO;. In this calculation we used a new and more reliable 
determination (2) of the variation of the counting efficiency of the spectro- 
meter with energy, and we assumed that when the original nickel comparison 
was made, the two peaks F and F’ were exactly superposed. 

An inspection of Fig. 1 will show that the strong coincidence peak F, which 
originally was not resolved, now appears to be caused by two y-rays, F and F’, 
the energies and intensities of which may be separately estimated. In addition, 
another y-ray F’’ appears between F and G. Owing to a breakdown in the 
apparatus, the part of the spectrum on the upper edge of the peak G was not 
plotted and unfortunately not repeated, and thus the energy of the y-ray G 
cannot be determined with any better accuracy than that obtained in the 
previous work. No attempt was made to explore the spectrum at energies 
below that of the y-ray J. 

At the upper end of the spectrum, the y-ray C which was previously ascribed 
to the ground state transition in K® is again indicated though its existence 





COINCIDENCES PER MINUTE 


° 





8.0 ; 85 9.0 Mev. 9.5 


Fic. 2. Higher energy coincidence spectrum produced by"potassium. Line width: 100 kev. 
Full circles, from natural potassium; open circles, from potassium enriched in K*°, The peak 
coincidence counting rate for the y-ray B is the same for both samples. 
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is not established with complete certainty. Above 7.7 Mev. (see Fig. 2, full 
circles) the coincidence spectrum is better defined than in the previous work. 
The peak (A) previously found at 8.03 Mev. (not shown in Fig. 2) is now 
known to be spurious. It was a ghost peak derived from the y-ray B by scatter- 
ing of electrons from the frames supporting the slits in the spectrometer (2) 
and this fault has now been greatly reduced. The peak A’ is now well defined 
but the counting rate in the peak A”’ is very low and its energy and intensity, 
therefore, are difficult to measure. 
MEASUREMENTS WITH THE SAMPLE ENRICHED IN K* 

To verify that the peaks A’ and A” are produced by neutron capture in 
K*°, the coincidence spectrum in this region was examined using a sample 
of anhydrous potassium carbonate which had been enriched* in K**. The 
amount of the material used, about 1 kgm., was approximately equal to that 
used in the natural potassium samples. A mass spectrographic analysis** 
showed that the K*° content had been increased by a factor of 2.2 (to an 
accuracy of 10%) and a chemical analysis showed that chlorine was present 
to the extent of 0.12% by weight (as potassium chloride). 

The height of the peak B produced by this sample was equal to that 
produced in similar flux conditions by the natural potassium sample (about 
50 coincidences per minute) and the contours of the two peaks could be very 
exactly superposed. The coincidence spectrum produced by the enriched 
sample beyond the peak B is shown by the open circles in Fig. 2. It will be 
seen that the peak A’ in the enriched sample is enhanced by a factor of nearly 
two when allowance is made for the tail of peak B (shown by the broken 
line). It is therefore clear that A’ must be produced by neutron capture in 
K*°. The peak at 8.61 Mev. which appears just above A’ may be identified 
(4) with the ground state y-ray in Cl*, for its energy is in agreement (within 
the limits of the experimental errors) with that of the chlorine y-ray (A) and 
its intensity is in close agreement with that expected from the known chlorine 
content of the sample. Unfortunately, we were not able to verify that this 
y-ray is due to chlorine by measuring chlorine y-rays of lower energy, for 
one of them, at 7.77 Mev., lies under the peak B and is masked by it, and the 
next, at 7.42 Mev., was not covered by the range of measurements made 
with the enriched sample. There seems to be no doubt that the low intensity 
peak A’’ must also be caused by neutron capture in K*® for, although the 
enhancement of this peak in the enriched sample is not very clear, there are 
no strong y-rays emitted at this energy by any of the possible impurities which 
the material might contain. 

The present measurements with the natural potassium sample were 
extended to 10.2 Mev., which is close to the neutron binding energy of K“. 
No y-rays were detected between the energy of A” and this value. We should 


*This sample was made from a quantity of potassium chloride which had been exposed for a 
period of nearly a year to the neutron flux of the NRX pile. We are indebted to the Isotope Division, 
Atomic Energy of Canada Limited, for the use of this enriched material and for the conversion 
of the chloride to carbonate. 

**We are indebted to Mr. W. H. Walker of the Department of Chemistry, McMaster University, 
Hamilton, Ontario, for this measurement. 
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not expect to detect the ground state y-ray in K*!, for this y-ray, like that in 
K*, should be 1/2 and relatively weak. The two y-rays A’ and A” presumably 
produce excited states in K*! at 1.8 and at 0.85 Mev., although no excited 
states are known at these energies. The intensities are roughly 14 and 4 
photons respectively per 100 captures in K*°. 

DISCUSSION OF THE y-RAY B 

We have already pointed out that the intensity of the y-ray B is anoma- 
lously high if this y-ray is indeed that corresponding to the direct M2 transition 
to the ground state of K*°. Its energy, 7.757 + 0.008 Mev., is in agreement 
(within the stated errors) with the neutron binding energy of K*°® as deter- 
mined by Sailor (6), viz. 7.71 + 0.08 Mev., and within the limits of the errors, 
80 kev., there is no evidence in Fig. 1 for the existence of another y-ray with 
an energy slightly above that of B which instead could be the ground state 
transition. 

The recent discovery of a low-lying state in K*° near 30 kev. accounts for 
the difficulty, for the energy of B is in close agreement with that of the 
transition which produces this state directly.* In both K* and K® the odd 
proton and odd neutron are in d3,2 and f7,2 states respectively. In the ground 
state of K®, which has a spin of 2 units (1), it is clear that the total angular 
momenta of the odd neutron and proton are antiparallel. The spin of K* is 
4 units which indicates that for some reason, the antiparallel condition of 
the spins does not give the most stable configuration. We should expect, 
therefore, the excited state of K*° to be the analogue of the ground state 
of K*, i.e., that it should have a spin of 2 units and odd parity. The y-ray (B) 
emitted in the direct transition to the 30 kev. state should then be of the £1 
type regardless of the spin of the compound state produced by neutron 
capture. This interpretation is consistent with the high intensity of y-ray B 
for it has been shown elsewhere (3) that El y-rays, in this energy range, are 
strong. It follows that the transition from the 30 kev. state to the ground 
state should be of £2 type, and according to Weisskopf’s formula, it should 
have a lifetime (when account is taken of internal conversion) of about 107° 
sec. (7). No isomeric state formed by neutron capture has been found in 
potassium nor would it have been readily detected for such a lifetime is too 
short for direct activation measurements. 
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A QUANTUM-STATISTICAL THEORY OF LIQUID HELIUM! 


By K. BuCKTHOUGHT 


ABSTRACT 


A model for liquid He‘ is set up in which the atoms are regarded as an assembly 
of Bloch waves following Bose-Einstein statistics. The potential field inside the 
liquid is estimated. Numerical calculation of the thermodynamic functions 
shows fair agreement with experiment. The A-temperature and its variation with 
pressure are calculated, and are in good agreement with observed values. The 
model is also applied to liquid He’. The variation of entropy with temperature is 
calculated. 


1. INTRODUCTION 


London (16a) and Tisza (18) have attempted to explain the anomalous prop- 
erties of liquid He* by the postulate that an assembly of He* atoms obeys 
Bose-Einstein statistics. Their model for the liquid is an ideal Bose-Einstein gas, 
and they ignore the strong interatomic forces which exist in the liquid. 

Phenomenological models in which the He* atoms are regarded as Bloch 
waves (2) moving through a quasi-crystalline potential field have been pro- 
posed by London (16) and Goldstein (6). These models incorporate no infor- 
mation about the interatomic potential, although the latter is known with fair 
accuracy from the equation of state of He‘ gas (4, 11, 17). 

In the present theory an attempt is made to evaluate the mean potential 
field encountered by each He atom in the liquid. The wave function and 
potential field are regarded as functions solely of the coordinates of the single 
atom, in the sense of the Hartree self-consistent field approximation for atoms 
and molecules (7). The wave function for the helium atom is then of the Bloch 
wave type, i.e. a plane wave whose amplitude is modulated by the potential 
field due to all the other atoms. The assembly of Bloch waves, regarded as 
independent systems, is to obey Bose-Einstein statistics. 

2. SPACE-DISTRIBUTION OF ATOMS IN LIQUID HE‘ 

Owing to the disordering effect of the large zero-point energy, He* cannot 
form a space-ordered crystalline lattice, even at 0° K., except under high pres- 
sure. However, the potential field encountered by a single atom must have, on 
the average, some of the periodicity characteristic of a crystalline lattice field. 
London (16) has investigated the stability of certain ordered configurations 
in He‘, and has demonstrated that most such configurations correspond to too 
high an internal energy to be probable (see also Appendix A). 

We expect that most disordered configurations which depart widely from a 
distribution consistent with the most stable one will also be highly improbable. 
This most stable configuration is that of a face-centered lattice, each of whose 
sites is occupied 50% of the time. Such a structure has been proposed by 
Keesom and Taconis as a result of their X-ray investigations of He* (10). The 
coordination number of such a structure is six, and the interatomic distance, for 
the observed density of liquid He‘, is therefore 3.2A. 


1 Manuscript received June 10, 1958. 
Contribution from the University of Toronto, Toronto, Ont. 
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We therefore propose to replace the actual field inside liquid He‘ by that 
smoothed-out field which would result from the following model. The atoms are 
to be represented by Bloch waves of the form 


(1] v(t) = u(r)e™* 


where u(r) is a function having a three-dimensional periodicity in the space 
vector r, and k is a constant vector playing the role of a set of three quantum 
numbers, analogous to the wave vector for a plane wave. 

We assume that the maxima of u(r) occur at points of aspace lattice consist- 
ent with the X-ray experiments and the theoretical considerations raised by 
London. 

For self-consistency, we should require that the wave functions used in 
evaluating the potential be the same as the above type (equation 1). Such a 
program would require a prohibitive amount of calculation, since the assign- 
ment of atoms of the assembly to various wave functions, given by different 
values of k, depends on the temperature. Thus a separate calculation should be 
made for each temperature. However, we are interested in a small range of 
temperatures (0-2.5° K.) for which the thermal energy is about one tenth of the 
zero-point energy. For the purpose of calculating the potential field, we have 
replaced the actual set of wave functions, with finite but small k-values, by the 
wave functions of lowest energy, with k = 0. These wave functions have been 
calculated, using the Wigner-Seitz method (20). The potential field employed 
for this calculation was the spherically symmetrical one used by Lennard-Jones 
and Devonshire in their theory of the liquid state (15). The Wigner-Seitz wave 
function was calculated from the Schrédinger equation by applying the WBK 
method (3a, 12, 19). The wave function, in terms of distance from a lattice- 
point, is shown in Fig. 1. 
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Fic. 1. Wigner-Seitz wave function for liquid He*. | ; 
(Units of probability density: atoms per cubic Angstrom unit.) 
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3. THE POTENTIAL FIELD INSIDE LIQUID HE‘ 
The potential field was calculated from 
N 2 . 
[2] Vir) = > J he J eulvsltar, 


jFi 
1 





where V,(r,;) is the hypothetical fixed potential field encountered by the 7 
atom of N atoms in the liquid, v;;(7;;) is the interaction energy between the i” 
and j" atoms, ¥4(r;) is the Wigner-Seitz wave function for the j‘" atom, and the 
integration is taken over all the space available to the liquid. In practice, 
because of the crystal symmetry it is necessary to calculate this function only 
for a small volume surrounding one lattice-point. At each point of this volume, 
there are contributions from wave functions whose maxima of probability 
density lie at neighboring lattice-points. Because the pair potential falls off 
very rapidly with increasing distance, the contributions from only a few (about 
10) wave functions are significant at a given point. 

The pair-potential takes on large positive values for interatomic distances 
less than 2.6 A. Here the mutual wave function of two He atoms must go to 
zero, imposing an additional correlation effect outside the scope of the Hartree 
method. This has been taken into account crudely by excluding all contribu- 
tions to V; inside a sphere of radius 2.6 A, centered around the field point r;. 
Also we assume that the 7‘ lattice point is vacant. 

The pair-potential used was that of Buckingham, Massey, and Hamilton (4). 
There is some evidence that this function gives the best fit to the data for He! 
gas (11, 17). 

4. SOLUTION OF THE SCHRODINGER EQUATION 

Solutions of the Schrédinger equation in which the potential has a three- 
dimensional periodicity in space have been discussed by many authors (30, 5, 
21). We shall adopt the procedure and notation of Wilson (21). 

We expand the potential and one-atom wave function in a three-dimensional 
Fourier series based on the reciprocal lattice. Denoting the reciprocal lattice 
unit vectors by b;, be, b;, and the general reciprocal lattice vector by 
@, = gibi + gob. + g;b;, where the g; are integers, we may write the potential 
field V(r) and wave function y,(r) in the form 


3) 


[3] Vir) = >> Veexp (27ig,.r), 
919293=—m@ 

[4] vu(r) = >> Cpe exp (27ig,. 1). 
919293=—c 


Substituting these in the Schrédinger equation for one atom, we have an infinite 


set of linear equations 


[5] | = V o00 —- . (k + 2 xe)" | Ce = :s Ce_hVh. 


2m h #000) 


These equations are homogeneous in the Cg as unknowns, and have solutions 
only if the infinite determinant for the coefficients of the C, vanishes. This 


determines the eigenvalue E, as a function of k. 
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For a free atom we would have 


WP 
ono ee 
2m ' 


E= V 000 -- 


ik.r 
yk = const. X e 
Using these as a zero-order approximation in equation 5, we obtain for a first 
approximation, 


9 


m |Vgl" 


. + 77 ho Sa) gal as 
[6] E= Vooo + k wip» [? +ék: 


2m 2 2 7S 


This holds if all the | C,| are very much less than | Cooo |, i.e. if the wave 
function does not depart much from that for a plane wave. 

For the energy level density, given by g(£) dE, that is the number of 
quantum states in the energy range E to E + dE, we have 


E+dE 


(7] g(E) dE = as J | { dk, dko dk. 
7 


We have evaluated this by expanding the last series of terms of equation 6 ina 
power series in the k;, taking advantage of the fact that 
gibik; < w| By. kK |. 


For any term corresponding to a set of integers g1, go, gs, there is a corresponding 
term given by —gi, —g2, —gs; in the summation over all values of g, terms in 
odd powers of the k; consequently disappear. Similarly, summation over all 
permutations of the + and — signs in +g), +g2, +g3 cancels out all cross 
terms k;k;, where 747. We get finally 


—— hk” = 4m Jae I Re 

E = Voo + om ah 2 lo. >» ba J a) \ xG- 7 3 rG* 
[8] ee ae \ 
ee we 


where G = | ga | 
In this rapidly converging series the ratio of successive terms in k; is of the 
order of 10-! at most, Thus it will suffice to keep only those terms shown above. 
We have evaluated equation 7 approximately, using expression 8 for £, 


giving 


re 9.V/9 3/2 512m? Vel2(E — Ep)?” 

( t( fe = ah p ie) | yy. 4 a a, 2 ( as 40) 
[9] g(k ) h° a (E Eo) a Shia? z G ’ 
| h* Val" 

where a=1- Fk at 


, ‘ i’ ~ |Vel 
Eo = Vooo — - ry as . 


9 
mmr 
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The last expression is the net energy at 0° K. With the data of Appendix B, 
resulting from the numerical evaluation of the Fourier coefficients for the 
potential field, we find Ey = — 10.8 cal./mol., as compared with the accepted 
value of —14.2 cal./mol. (9). 
5. STATISTICAL MECHANICS OF LIQUID HE4 

The He‘ atoms are independent systems in this model, each moving in its 
own potential field. The distribution law for a Bose-Einstein assembly of 
independent systems is 


(EAE 


-(E=#)/RT : 
e —1 





[10] n(E) dE = 


where yu is the chemical potential (per mol.). E is now counted from the lowest 
quantum state as energy zero, and we will allow for the zero-point energy by 
including it as a constant U», in the internal energy of the assembly taken as a 
whole. 

To obtain yu as a function of temperature, we write the condition that the 
total number of atoms in volume V be a constant, i.e. 


2 eas cee a E) dE 
[11] N= : n(E) dE = J ———S ar 
” oe — JE 
The total internal energy of the assembly is then 
; ; ; @  Eg(E\)dE 
[12] U = Ut J a 
o¢ i 


From these two equations all the thermodynamic functions may be deduced. 
London (16a) has shown that below a certain temperature 7) given by 
V = . g(E)dE 


E/RT)” 
e 


— |] 


0 


a finite fraction of the systems ‘‘condense’’ into the lowest energy state. The 
number of systems in this state, for temperatures less than 7), is 
1 


No = yar" 

e = 

The thermodynamic functions below the \-point are determined by putting 
nu = 0 in equation 12, as shown by London. 

The thermodynamic functions were evaluated using equation 9 for the 
energy-level density. The entropy and specific heat are shown in Fig. 2, as 
functions of temperature. 

The internal energy was calculated as a function of volume, in the following 
way. In the analysis of the potential field the contributions from the attractive 
and repulsive parts of the potential were kept separate. For the purpose of 
calculating their variation with density, they were assumed to follow a Lennard- 


Jones potential (14), viz. 
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Fic. 2. Entropy and specific heat of liquid He?. 

-— Calculated values. -— Observed values. 
where &, 7 are positive constants, 7;; is the interatomic distance. This relation 
enables us to deduce the variation in the potential field resulting from changes 
in interatomic distance. 

Thus the compressibility K 7 and expansion coefficient 8 could be calculated 
from the thermodynamic formulae 
1/Kr = V@F/dV?, where F = U — TS, 
T(dp/dT)y — p = (dU/dV)z, 
B = K7(dp/dT)y. 


The results of these calculations are shown in Fig. 3. 

Similarly, the variation of 7, with density was calculated, and its variation 
with pressure deduced from the calculated compressibility (Fig. 4). At zero 
pressure we find 7, = 2.23° K., in fair agreement with the observed value of 
2.186° K. (9). 

The temperature variation of the compressibility of liquid He‘ is anomalous, 
in that there is a discontinuity at the A-point. It is therefore particularly satis- 


‘factory that this discontinuity, together with the rise in compressibility toward 


the A-point, is indicated in this theory. The negative expansion coefficient 
below 7) and the sign of d7\/dp are also accounted for. 

Exact agreement with experiment is not to be expected from a theory which 
takes into account the interatomic interactions in an oversimplified way. 
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Fic. 3. Compressibility and expansion coefficient of liquid Het. 
——————— Calculated values. --------—~—- Observed values. 
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Fic. 4. Variation of \-temperature with pressure. 
———————. Calculated values. @ Observed values. 
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6. LIQUID HE? 

The same type of model was applied to the problem of liquid He*. In the 
absence of direct experimental evidence, it was assumed that the space distribu- 
tion is the same as that in liquid He‘, except for the larger interatomic distances, 
corresponding to the larger molar volume of He’. In fact, a calculation of the 
internal energy for various possible configurations shows that the only stable 
one is the 50% occupied face-centered lattice used for He* (Appendix A). 

The Fermi-Dirac statistics was applied to this case. A system obeying 
Fermi-Dirac statistics is characterized by a degeneracy temperature (3d) 


r,=3 tor. 
°" 5 2mR\82xV/) * 


Here the constant a, characteristic of the present model, has been introduced to 
take into account that the energy level density of the Bloch wave energy 
spectrum differs from that of an ideal gas. The value of a was calculated by a 
method analogous to that used for He*, and was found to be 0.94. We find 


3 


t * 
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Oo I 
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ie 
Fic. 5. Variation of entropy of liquid He* with temperature. 
—— Calculated values. @ Observed values. 

Ty) = 1.8° K. This gives a Fermi-Dirac entropy as shown in Fig. 5, which 
departs from the linear relationship, where the temperature is not small 
compared to 7». 

The internal energy at 0° K. was calculated to be —1 cal./mol., compared 
to an experimental value of —4.4 (1). 

7. DISCUSSION 

The concept of liquid He‘ as a Bose-Einstein assembly seems to provide a 
satisfactory basis for the explanation of the thermodynamic properties of the 
liquid. The striking difference in behavior of the two isotopes of helium is 
immediately explained, as in the ideal gas model proposed by London. But the 
model considered here corrects some of the defects of the ideal gas theory. The 
general behavior of the compressibility and expansion coefficient is explained, 
and the sign of d7,/dp is predicted correctly. 

But the specific heat curve has a much sharper maximum at the A-point than 
predicted. This points to some ordering process not taken into account— 
perhaps outside the realm of the Bose-Einstein statistics, and thus probably 
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present in He’ as well. For the latter, the exclusion principle presumably pre- 
vents a condensation in momentum space. 

An additional weakness of the present theory lies in its prediction of a specific 
heat varying as JT?” near absolute zero. According to the theory of Landau (13) 
He‘ behaves like a Debye solid near 0° K., with a specific heat going to zero 
as T*. The experimental evidence favors Landau’s theory (8). 

Although Landau has had considerable success in predicting the behavior of 
liquid He* near absolute zero, his theory fails to give a satisfactory account 
of the \-transition. An important feature of the theory presented here is its 
ability to deal with the discontinuities at the A-point. 
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APPENDIX A 
Potential and Kinetic Zero-point Energies for Various Types of Lattice 
The potential energy for each structure has been determined by summing up 
the contributions to the potential energy at a lattice-point, of nearest neigh- 
bors and next nearest neighbors. This energy has been multiplied by Avo- 
gadro’s number divided by two, to get the energy per mole, U. 
The zero-point energy was calculated from the formula (16d) 
Z = h?/8Mr’, 


where M is the mass of the atom, 7 is the mean “‘radius of action”’ of the atom, 
conceived of as confined to a sphere. r was calculated as follows: (a) if the 
number of nearest neighbors is 8 or 12, 7 is to be the nearest neighbor distance 
less 2.4 A; (b) if the number of neighbors is 4 or 6, 7 is to be the mean of the 
nearest neighbor distance less 2.4 A. If there are three nearby shells of lattice 
sites, the radii of all three are averaged. 

These criteria were adopted to allow for the lower kinetic energies associated 
with the more open structures, having lower coordination numbers. 

In Table I all energies are in cal./mol. The Keesom-Taconis lattice is that 
proposed in Reference (10). The last two types listed are essentially the same, 
except that the “‘holes’’ of one are the lattice-sites of the other. 








Lattice type 


Cubic close-packed 
Body-centered cubic 
Simple cubic 
Diamond 

50% body-centered 
Keesom-Taconis 
50% face-centered 
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APPENDIX B 
Analysis of Potential Field Inside Liquid He* 

We may conveniently represent the lattice used as a basis for the calculation, 
by a set of simple cubic lattices having each the same lattice constant, but dis- 
placed with respect to one another. If the Fourier coefficient for the 50% 
occupied face-centered lattice is Vg, and that for any one of the simple cubic 
ones is ?g, we have 

Vg = g{1 + exp wi(g: + go) + exp wi(ge + gs) +... 
+ 2 cos $a(g; + go) + 2 cos 34r(go + gs) +... 
+ exp }m1(g; — go + 2g3) + exp 3mt(go — gs + 2g1) +...} 


giving 





V 000 = 16 vo00, | Vin | = 6 | Oy |, | Vor | = 16 | V229 


| V 00 | ™ V no es | 


The vg were determined from numerical integration of 


v% = : J J fre exp —27i(g,.1) dr 


over the volume of the liquid 7, giving the following values (cal./mol.) 
Vooo = —5.6, Vin = 4.9, Vizg = 26.4, Viss = 0. 
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AN ATOMIC MASS STUDY OF NUCLEAR SHELL STRUCTURE IN 
THE REGION 28 < n < 50 AND 28 <Z < 40' 


By BenyAMIN G. HoGc? AND HENRY E. DucKWoRTH 


ABSTRACT 


New mass spectrographic masses are reported for some of the isotopes of 
gallium, germanium, arsenic, and selenium. These results are combined with 
previously reported atomic mass data in order to study atomic mass trends in 
the region 28 <n < 50 and 28 < Z < 40. The only pronounced mass effects 
which occur in this region are those associated with 28 and 50 nucleons. We 
have found no extra stability which can be associated with the filling of the 
fsj2, Psj2, and p1/2 subshells. A list of mass values used in this study is included. 


I. INTRODUCTION 

Some time ago one of us gave a summary (10) of the atomic mass evidence 
for the existence of certain notably stable nuclear configurations. These stable 
configurations or ‘nuclear shells’ were attributed to the presence in the nucleus 
of 28, 50, or 82 protons, or 50, 82, or 126 neutrons, in agreement with the pre- 
dictions of the single particle nuclear model of Mayer (34, 35) and Haxel, Jensen, 
and Suess (24). Since that time the mass effect associated with the closure of the 
28- and 50-proton shells has been confirmed by Nier’s group at the University 
of Minnesota, who have also demonstrated (8) that there is atomic mass evidence 
for a 28-neutron shell. 

These pronounced mass effects are a// explained in terms of a spin-orbit 
coupling which results in large splitting of a nuclear energy level with orbital 
angular momentum / into two distinct levels } = /+ 3. Thus, the 28, 50, 82, 
and 126 shells mark the divisions between the f7/2—f5/2, 29/2-27/2) Ai1j2-A9/2, and 
113/2-111;2 groups, respectively. 

It is now of interest to determine the degree of stability which is associated 
with the filling of nuclear shells other than those resulting from strong spin-orbit 
coupling. Probably the best region in which to study this matter is that where 
28 <n < 50 and 28 < Z < 40, which extends from nickel to zirconium. Here 
is a region which is almost entirely free of major shell effects and, for this reason, 
is well suited to the study of effects which are presumably smaller than those 
which have been observed to date. At the same time, it is a region in which 
three subshells, namely, f5/2, 3/2, and p1,2, are filled. In this paper we report new 
atomic mass measurements in this region, and combine these with existing data 
in order to investigate the nuclear stability variations in which we are interested. 

II. EXPERIMENTAL 

The new atomic mass measurements upon which this paper is based were 
made with a Dempster-type double focusing mass spectrograph. Both the 
instrument (9) and the method (26) have been previously described in some 
detail. 

1 Manuscript received May 20, 1953. 
Contribution from Department of Physics, McMaster University, Hamilton, Ontario. This 
work has been supported by the National Research Council of Canada, the Research Council of 
Ontario, and by the Office of Scientific Research, Air Research and Development Command, 


United States Air Force. 
2 Holder of a National Research Council of Canada Fellowship. 
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The mass differences obtained in the present work are given in Table I. 


TABLE I 
List OF NEW ATOMIC MASS DIFFERENCES 











Mass difference | Previous 
Nuclides in m.m.u. | measurements 
Pb?" — 3Ga* | 196.6 + 6 
Ce — 2Ge7° | 56.8 + 6 | 59.524 4° 
Prt — 3T7i"” 50.9 + 3 5134+ 4° 
Nd — 3Ti8 66.7 + 3 
Nd'#4 — 2Ge”? 66.3 + 2 | 69.8 + 10° 
3Ti*s 0444 
3Tiv 85.2 + 5 | 77.4 +17° 
2As% 77.1+4 
3TiP? 8.1+6 
2Ge7® 76.2 + 2 
2Se76 81.4+6 | 
Se7® 26+3 
Gd"6 — 3Cri2 100.4 + 4 99.4 + 17° 
Gd®* — 2Se78 87.5 +2 
2Se78 — 3Cr*? 12.9+4 
4 Reference (12). 6 Reference (20). 


III. STANDARD MASSES 

In computing the masses of the nuclides with which we are concerned in this 
paper we have regarded the masses of certain of the lighter nuclides as sub- 
standards, in particular, H' and C”. This has necessitated choosing between the 
masses for these nuclides computed by Li ef a/. (29) on the basis of nuclear re- 
action data, and those obtained by mass spectroscopic methods by Nier and his 
collaborators (7) and by Ogata and Matsuda (37). Since we shall be making use 
of many other masses obtained by Nier and his group, we have chosen to base 
all of our present calculations on their masses for the substandards. These 
values are 

H!' = 1.008146 + 3 a.m.u., 
CY = 12.003842 + 4 a.m.u. 

The accuracy of the new mass comparisons described in this paper is such 
that the use of the reaction-derived values of H! = 1.008142 + 3 a.m.u. and 
C® = 12.003804 + 17 a.m.u. would generally make inconsequential changes in 
the atomic masses computed therefrom. These changes would never be outside 
the probable errors associated with the masses in question. 

IV. COMPUTATION OF ATOMIC MASSES 
(a) Argon, Potassium, Calcium, Scandium, Titanium, Vanadium, Chromium, 
Manganese, Iron, Nickel, Copper, and Zinc 

The masses of all the stable isotopes of these elements, with the exception of 
K* and Ca**, have been accurately measured by Collins, Nier, and Johnson (7, 8) 
and have been used by them to study the stability effects associated with the 
28-proton and 28-neutron shells. We had previously measured a number of 
these masses in our laboratory and find that the agreement is generally good 
between the new measurements and our older, less accurate ones. The two sets 
of measurements of the eight nuclides which have been studied in both labora- 
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tories are shown in Table II. These masses are based in each case on H! 
= 1.008146 + 3 a.m.u. and C® = 12.003842 + 4 a.m.u. 
TABLE II 


COMPARISON OF MEASURED MASSES FOR SOME NUCLIDES OF 
TITANIUM, CHROMIUM, IRON, AND NICKEL 


Collins, Nier, and Johnson? | Duckworth et al. 
Nuclide | (a.m.u.) | (a.m.u.) 


Tis | 47.96317 + 6 | 47.9630 + 3 ° 
Cro 49,96210 + 6 49.9602 + 2 ° 
Cr? 51.95707 + 9 | 651.9571 + 2 4 
Fe* =| 53.95704 + 5 | 53.9567 + 2 ¢ 
Fes | 55.95272 + 10 55.95279 + 17 ° 
Nis 57.95345 + 10 57.95362 +18 ¢ 
Ni®o 59.94901 + 29 | 63.94931 + 14 / 
Nis | 63.94755 + 8 | 63.94728 +16 


: Reference (7, 8). 

’ This is a new value resulting from a remeasurement of our old 
O' — Ti plates. ; 

° This is the mass we originally reported and is known to be in 
error. See Reference (11). 

“Reference (16). 

Reference (14). 

’ Reference (13). 

9 Reference (12). 

In the light of this agreement, we conclude that the Minnesota mass values 
and ours are equivalent for this part of the mass table, albeit the Minnesota 
values are more accurate. We are, therefore, adopting in this paper the mass 
values obtained by Nier and his collaborators for the elements argon to zinc, 
except for the nuclides Fe®®, Ni’, and Ni®. In the case of this trio, weighted 
means seem to be called for: these are Fe®® = 55.95274 + 9 a.m.u., Ni®3 = 
57.95349 + 9 a.m.u., and Ni® = 59.94925 + 13 a.m.u. 


(b) Cobalt 

Cobalt consists of a single stable isotope at mass number 59. Its mass can be 
deduced from Ni® = 59.94925 + 13 a.m.u. (see above), knowing that Co® — 
Co*® = 1.00070 + 5 a.m.u. (3, 23, 4) and Co® — Ni® = 0.00302 a.m.u. (42), 
to be Co®® = 58.95136 + 14 a.m.u. A second value for its mass can be deduced 
from Ni®® = 57.95349 + 9 a.m.u. by using Ni5® — Ni’ = 0.99930 + 1 a.m.u. 
(23, 28), and McCue’s figure (31) of 1.888 + 0.003 Mev. for the Co®(p,n) 
threshold. These data lead to Co®*® = 58.95164 + 9 a.m.u., which is in fair 
agreement with the Ni®-derived value. The weighted mean for Co®® = 58.95157 
+ 10 a.m.u. 


(c) Gallium 

Gallium exists in two stable isotopic forms at mass numbers 69 (60.2%) and 
71 (39.8%). We have determined the Pb?°? — 3 Ga® mass difference as noted in 
Table I. Using Pb?°’ = 207.0405 + 10 a.m.u. (40), we compute Ga®® = 68.9480 
+ 4 a.m.u. This value agrees well with the mass of Ga®® computed from the 
Minnesota masses for Zn®* and Zn7° and appropriate transmutation data, as 
will now be shown. First, one can combine Zn® = 67.94686 + 7 a.m.u. (8) with 
the 10.10 + 0.2 Mev. Ga® (vy, 2) threshold (39) and the Ga® — Zn*® decay 
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energy, known to be 2.90 Mev. (42), to obtain Ga®® = 68.94810 + 22 a.m.u. 
Then, also, using Zn7° = 69.94779 + 6 a.m.u. (8), the Zn?°(y,7) threshold of 
9.2 + 0.2 Mev. (22), and the mass difference Zn®’—Ga*® = 0.00092 (42), one 
calculates Ga®® = 68.94777 + 22 a.m.u. 

These three different values for the mass of Ga®® are in good agreement and 
lead to a weighted mean of Ga® = 68.94794 + 16 a.m.u. It is not possible at 
present to compute a mass for Ga”, the other stable isotope of gallium. 


(d) Germanium, Arsenic, and Selenium 

These three elements are being discussed together since they are linked by 
certain key transmutation data in addition to a mass spectrographic mass 
comparison. 

Germanium has five stable isotopes—70(20.55%), 72(27.37%), 73(7.61%), 
74(36.74%), and 76(7.67%). Arsenic has but one at mass number 75 
while selenium consists of six—74(0.87%), 76(9.02%), 77(7.58%), 78(23.52%), 
80 (49.82%), and 82(9.19%). 

The mass of As7° is found from the 2As7> — 3Ti*® mass difference listed in 
Table I and the mass Ti*® = 49.96077 + 4 a.m.u. (8) to be As?> = 74.94522 
+ 30 a.m.u. Kinsey and Bartholomew give 7.30 + 0.04 Mev. for the energy (2) 
of the most energetic gamma ray following slow neutron capture in As’. This 
will be assumed to represent the transition to the ground state of As’®. Combining 
these data with 3.12 Mev. (42) for the As7® — Se7® decay energy, one obtains 
Se7® = 75.94801 + 30a.m.u. A second value for this mass will now be computed. 

The mass of Se?’ is found from Table I and Cr® = 51.95707 + 9 a.m.u. (8) 
to be Se’® = 77.94210 + 22 a.m.u. From the Chalk River (,y) experiments 
the Se75(y,7) threshold has been found to be 10.483 + 0.014 Mev. (2), whence 
Se” = 76.94438 + 22 a.m.u. Further, Kinsey and Bartholomew give a probable 
value of 7.416 + 0.009 Mev. for the Se77(y,7) threshold (2). This leads to 
Se’® = 75.94337 + 22 a.m.u. 

The weighted mean of these two values is Se’® = 75.94324 + 18 a.m.u. 
Assuming that this represents our best value for Se’®, we now adjust the previous 
values of As”, Se77, and Se7® to 74.9454 + 3 a.m.u., 76.9442 + 2 a.m.u., and 
77.9420 + 2 a.m.u., respectively. 

The threshold for the As’*(y,7) reaction is 10.20 + 0.2 Mev. (38, 39) and 
the resulting As’ decays both by positron and negatron emission. There is a 
total decay energy of 2.54 Mev. in the positron branch and 1.36 Mev. in the 
negatron branch (42). These data lead to Se™ = 73.9460 + 4 a.m.u. and 
Ge™ = 73.9447 + 4 a.m.u. 

One may now use these masses for Se™4 and Se’* together with ratios of the 
isotopic masses of selenium found in the microwave absorption experiments of 
Geschwind and others (19) to compute Se*® = 79.9421 + 6 and Se® = 81.9431 
+ 6 a.m.u. 

Turning now to germanium, the mass of Ge™ has been computed above, and 
that of Ge” is found from Table I and Ti*® = 47.96317 + 6 a.m.u. (8) to be 
71.94493 + 20 a.m.u. These masses set the scale for the microwave-derived 
mass ratios for germanium (18) and lead through them to Ge? = 69.9464 
+ 7 a.m.u. and Ge? = 75.9456 + 7 a.m.u. These computed masses for Ge7° 
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and Ge’ will now be compared with ones that we have determined experi- 


mentally. 

Using Ti‘? = 46.96668 + 10 a.m.u. (8), 9.5 + 0.2 Mev. for the Pr''(y,7) 
threshold (22, 38), 3.22 Mev. for the total Pr'*#° — Ce™® decay energy (42), 
together with the Pr’! — 3Ti‘7 and Cel® — 2Ge7° mass differences given in 
Table I, one obtains Ge? = 69.9460 + 5 a.m.u. The weighted mean of the two 
values for Ge” is 69.9461 + 4 a.m.u. 

The mass of Ge7® may be obtained from Se7* using the mass spectrographic 
measurement of the Ge*® — Se7® mass difference listed in Table I. This leads to 
Ge7® = 75.9458 + 4 a.m.u. which, with the microwave-derived value, gives a 
weighted mean of Ge’® = 75.9457 + 4 a.m.u. 

These results indicate that the Ge7® — Se7® mass difference is 0.0025 + 
4a.m.u., or 2.3 + 0.4 Mev. This leaves 0.8 Mev. available for the decay of As7® 
by K-capture, if the total energy for negatron decay is correctly given by 3.1 
Mev. (42). Experimentally, the ratio of positrons to negatrons has been set at 
less than 0.007% (42). 

(e) Bromine 

There are two stable isotopes of bromine at mass numbers 79(50.52°%) and 
81(49.48%). We have not obtained mass spectrographic values for the masses 
of these nuclides. However, it is possible to compute the mass of Br?® from Se’ 
and that of Br*! from both Se*® and Se®. 

In the case of Br7’, one uses the Q = 4.55 + 0.2 Mev. for the Se7*(p,7) reaction 
(5) and 10.7 + 0.2 Mev. (32, 39) for the Br7°(y,”) threshold. In this way one 
computes Br’? = 78.9435 + 4 a.m.u. 

In obtaining Br*! from Se**, one makes use of Q = 2.57 + 0.2 Mev. for the 
Se*°(p,2) reaction (5) and 10.07 + 0.2 Mev. for the Br''(y, n) threshold (32, 39). 
Thus, Br*! = 80.9422 + 7 a.m.u. In working from Se*, one knows the Se*(-y,7) 
threshold to be 9.8 + 0.5 Mev. (1) and the total Se*! — Br*! decay energy to be 
1.5 Mev. (42). In this way, Br*! = 80.9430 + 8 a.m.u. The weighted mean of 
these two values is 80.9425 + 5 a.m.u. 

There is good agreement between these computed masses for Br7® and Br! 
and those obtained by Hays, Richards, and Goudsmit, who found, using their 
time-of-flight mass spectrometer (25), that Br7® = 78.944 + 1 a.m.u., and 
Br®! = 80.943 + 1 a.m.u. 


(f) Krypton 

There” are six stable isotopes of krypton at mass numbers Rnd 
80(2.27%), 82(11.56%), 83(11.5%), 84(56.90%), and 86(17.37%). Collins has 
obte and: masses (6) for all of these isotopes. The masses of Kr a and Kr*¢ 
had been previously measured in our laboratory (27) and, nichinenth there was 
good agreement between our value and that of Collins for Kr*, there were large 
discrepancies in the cases of Kr® and Kr**, We have re-examined the plates from 
which these measurements were made and have remeasured them. The results 
of the new measurements on these plates are compared in Table III with Collins’ 
values. 

The actual doublet measurements upon which these new masses are based are: 
2C3H; — Kr®, Am = 0.1653 +3 a.m.u.; 2C3Hg — Kr™, Am = 0.1826 + 3 
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TABLE III 


New measurements 
Nuclide Collins (a.m.u.) 


Kr*? 81.93961 +11 | 81.9392 + 4 
Krist 83.93836 + 9 | 83.9382 + 4 
Krsé | -85.93820+ 8 | 85.9384 + 4 


a.m.u.; ;Xe”* — Kr86, Am = 0.0255 + 3 a.m.u. In computing the mass of Kr*® 
we have used Xe”? = 128.94594 + 15 a.m.u., which is a weighted mean of 
Halsted’s 128.9460 + 2 a.m.u. (21) and our earlier value (27) of 128.94570 
+ 26 a.m.u. 

These new measurements remove two worrisome discrepancies and justify 
the assumption that Collins’ masses for krypton are compatible with masses 


got by us in this region. 


(g) Rubidium 

The only masses available for the two stable isotopes of rubidium, Rb** and 
Rb‘’, have been obtained by Hays, Richards, and Goudsmit (25). These do not 
appear to be consistent with the other masses which have been determined in 
this region and we are, consequently, not including them in the present set. 


(h) Strontium, Yttrium, and Zirconium 

Strontium exists in four stable isotopic forms at mass numbers 84(0.56%), 
86(9.86%), 87(7.02%), and 88(82.56%). Preliminary measurements of the 
masses of Sr8* and Sr*8 were made in 1937 by Mattauch (33). Recently, the mass 
of Sr®6 has been found, by a study of the C;H; — Sr**and C.0OH; — Sr** doublets, 
to be 85.93552 + 52 a.m.u., and that of Sr’* has been found, by a study of the 
CO. — Sr** doublet, to be 87.93364 + 35 a.m.u. (15). 

The Sr’? — Sr** mass difference can be computed using Harvey's Q of 6.29 
+ 0.2 Mev. (23) for the Sr**(d,p) reaction, and, also, using 8.40 + 0.20 Mev. (39) 
for the Sr*“(y,2) threshold, as measured by Sher ef al. These two values are 
0.99984 + 21 a.m.u. and 0.99996 + 21 a.m.u. respectively, giving a weighted 
mean of 0.99990 + 15 a.m.u. 

The Sr*8 — Sr‘? mass difference can be computed from the Sr**(y,7) threshold. 
This has been measured by Sher ef a/. to be 11.15 + 0.20 Mev. (39), corre- 
sponding to a mass difference of 0.99700 + 21 a.m.u. This datum, combined 
with the information in the previous paragraph, leads to a Sr** — Sr*® mass 
difference of 1.99690 + 29 a.m.u., which is in poor agreement with the mass 
spectrographic value of 1.99812 + 63 a.m.u. This difficulty will be dealt with 
later after a discussion of the mass of Zr®® and the Zr*® — Sr** mass difference, 
which information is useful in assessing the reliability of the masses of 
Sr®* and Sr®*. 

There are five stable isotopes of zirconium at mass numbers 90(51.46%), 
91(11.23%), 92(17.11%), 94(17.40%), and 96(2.80%). The mass of Zr has 
been determined mass spectrographically (13), by photographing the Si®° — Zr®° 
doublet, to be 89.93286 + 36 a.m.u. To ascertain the compatibility of this value 
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with the Sr** mass given above, one can compute the Zr®® — Sr** mass difference 
as follows: 

Sr8? — Sr88 = 1.00194 + 21 a.m.u. Reference (23) 

Sr’ — Y*? = 0.00161 Reference (42) 





Y — Sri = 1.00033 + 21 a.m.u. 
Zr” — Zr® = 0.99620 + 32 Reference (38) 
Zr®® — YY = 0.00303 Reference (42) 


Zr — Y® = 0.99923 + 32 





Zr® Zr**™ = 0.99578 + 16 References (22, 38) 
Zr8™ — Y89 = 0.00364 Reference (42) 


(c) Zr — Y% = 0.99942 + 16 
In this way Zr*® — Sr*3 = 1.99968 + 25 a.m.u.—(a) plus weighted means of 
(b) and (c). 

This mass difference is to be compared with the mass spectrographic value of 
1.9992 + 5 a.m.u. Since the agreement here is satisfactory, whereas it was not 
so in the case of the Sr** — Sr** difference, we are assuming that the Sr** mass 
is in error. We should have liked to remeasure the Sr** doublets in order to 
justify this assumption, but the plates upon which they are recorded have 
unfortunately been misplaced. We have, therefore, adjusted the values for Sr** 
and Zr®® in the light of the transmutation value for the Zr°° — Sr** mass dif- 
ference. The results of this adjustment are shown in Table IV. From these ad- 
justed values, using the transmutation results quoted above, one computes 
Sr®6 = 85.9366 + 6 a.m.u., Sr’? = 86.9365 + 4 a.m.u., and Y%* = 88.9336 
+ 4 a.m.u., the latter being the single stable isotope of yttrium. It is, unfortu- 
nately, not possible to compute a mass for the remaining stable isotope of 
strontium, namely Sr**, 

TABLE IV 
COMPARISON OF UNADJUSTED AND ADJUSTED MASSES 
FOR Sr°8 AND Zr*? 


Nuclides Tnadjusted values | Adjusted values 


Sr88 87.93364 + 35 87.9335 + 4 
Zt: Goas 1.99968 + 25 1.99958 + 20 
Zr 89.93286 + 36 89.9330 + 4 


Harvey has determined (23) the Q’s of the Zr®°(d,p) and Zr®'(d,p) reactions 
to be 4.93 + 0.05 Mev. and 6.50 + 0.10 Mev., respectively. These values lead 
to the masses Zr*! = 90.9343 + 4 a.m.u. and Zr® = 91.9339 + 4 a.m.u. The 
masses of Zr*4 and Zr** can be found from those of Mo“ and Mo*® (see next 
section) using Zr** — Mo = 0.00122 + 20 a.m.u. and Zr®* — Mo** = 0.00363 
+ 28 a.m.u. (17). The resulting values are Zr™ = 93.9365 + 5 a.m.u. and 
Zr = 95.9394 + 5 a.m.u. 
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(7) Nb, Mo”, Mo™, Mo, and Ru%® 

For the purpose of this paper it is also desirable to compute masses for Nb*, 
Mo”, Mo”, Mo®, and Ru®® in order to determine more accurately the stability 
associated with stable nuclides possessing 50 and 52 neutrons. 

The mass of Nb® can be deduced from that of Zr? using 8.7 + 0.2 Mev. for 
the Nb%(y,7z) threshold (39) and 1.8 + 0.2 Mev. for the Nb® — Zr® mass 
difference (17). The mass so deduced is Nb* = 92.9355 + 4 a.m.u. 

A recent measurement of the Mo®” — Zr mass difference has given 1.25 + 0.2 
Mev. (17), hence, Mo” = 91.9352 + 4 a.m.u. The mass of Mo™ can be calcu- 
lated from that of Nb®* knowing that Q = 5.03 + 0.10 Mev. for the Nb*(d,p) 
reaction (23) and that 1.26 Mev. represents the total Nb*t — Mo” decay 
energy (42). This leads to the value Mo” = 93.9353 + 4 a.m.u. The mass of 
Mo* has been found mass spectrographically to be Mo** = 95.9358 + 4 a.m.u. 
(18). 

Finally, the mass of Ru® is found, by combining the recent value Ru®* — Mo* 
2.8 + 0.2 Mev. (17) with the mass of Mo*, to be Ru®® = 95.9388 + 4 a.m.u. 


TABLE V 
ATOMIC MASSES ADOPTED IN THIS PAPER 


oS 


39.97513 + 3 Ni® 59.94925 +13 | Br*! 80.9425 
40.97490 + 4 Ni®! 60.94907 + 23 | Kr7® 77.94519 
39.97545 + 9 | Ni® 61.94681+ 9 | Kr*® 79.94246 
41.97216+ 4 | Ni® 63.94755 + | Kré? 81.93961 
42.97251 + 6 Cu® 62.94926 + 6 | Kr* 82.9419 
43.96924 + 6 | Cu® 64.948385 + 6 | Kr* 83.93836 
47.96778 + 1¢ Zn™ 63.99455 + | Kr8® 85.93820 
44.97010 5 Zn® 65.94722 + 6 | Rb® 

45.96697 Zn™ 66.94815 + 6 | Rb* 

46.96668 Zn® 67.94686 + | Se 

47.96317 Zn 69.94779 + 6 Sr** 85.9366 
48.96358 Ga®® 68.9478 | Sr’? 86.9365 
49.96077 Ga"! Sr*8 87.9335 
49.96330 Ge? 69.9461 Y* 88.9336 
50.96052 Ge? 71.9449 Zr® 89.9330 
49.96210 Ge Zr™ = 90.9343 
51.95707 Ge" 73.9447 Zr® 91.9339 
52.95772 Ge®® 75.9457 Zr™ =93.9365 
53.956¢ As® 74.9454 Zr =95.9394 
Se! 73.9460 Nb*® 92.9355 
Se7® 75.9432 Mo” 91.9352 
Se7? 76.9442 Mo” 93.9353 
Se*® 77.9420 Mo** 95.9358 
Se*® 79.9421 Ru 95.9388 
Se 81.9431 
Br®? 78.9435 
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V. DISCUSSION 
(a) Mass Values 
The atomic mass values adopted in this paper are shown in Table V. Each of 
these has been discussed in the previous text. 
(6) Representation of the Data 
Since Table V contains a fairly complete list of the stable isotopic masses in 
the region 28 < n < 50 and 28 < Z < 40, it should be possible to use this 
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information to study nuclear shell effects. In attempting to do this, we have 
represented the data in several different ways, the most revealing of which 
appears to be that shown in Figs. 1 and 2. In these figures we have plotted as 
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Fic. 1. <A plot of neutron number versus experimental deviation from semiempirical masses. 
The solid curves connect all even-A isotopes for the elements shown. The closed circles belong 
to even-Z and odd-n nuclides. As can be seen from Table V, the uncertainty in the plotted 
points is generally less than 0.5 m.m.u. and, in the case of the Minnesota mass values, is very 
much less than this figure. 


ordinate the difference between the masses computed from the semiempirical 
mass formula by Metropolis and Reitweisner (86) and the observed masses. 
These differences are shown in Fig. 1 as a function of the number of neutrons in 
the nucleus, and in Fig. 2 as a function of the number of protons. Similar plots 
have been employed by Low and Townes (30) and by Wapstra (41). 

(c) Analysis of Figs. 1 and 2 for Shell Structure Effects 

The large sharp peaks occurring at neutron number 50 and proton number 28 
are immediately apparent in Figs. 1 and 2, respectively. These have been previ- 
ously ascribed to the spin-orbit splitting of the gg/2-g7/2. neutron levels and 
frj-fs;2 proton levels. The mass effect associated with the 28-neutron magic 
number is also apparent in Fig. 1, although partially obscured by the 28-proton 
shell, which is a much larger effect. 

This strong spin-orbit coupling model of the nucleus (34, 35, 24) seems to be 
the most useful guide in the study of nuclear shell effects. On this scheme the 
fz). shell (completed at 28 nucleons) is followed by three subshells before the 
g9/2 shell (completed at 50 nucleons) begins. These subshells are f5,2, p3/2, and 
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p1;2 and, if completed in this order, they represent nucleon numbers 34, 38, and 
40, respectively. We shall examine Figs. 1 and 2 for mass effects at these numbers. 


MILLIMASS UNITS 


30 
PROTON NUMBER 


Fic. 2. A plot of proton number versus experimental deviation from semiempirical masses. 
The solid curves connect even-A nuclides for the isotonic numbers shown. The circles belong 
to nuclides having even Z and odd neutron number. The uncertainty in the points is the same 
as in Fig. 1. 


34 Nucleons 

A small effect may be trying to assert itself at 34 neutrons in the nickel isotope 
curve as shown in Fig. 1. However, if instead of the directly determined mass of 
Ni®, we use the value computed from Cu®, taking the Cu®(y,2) threshold to be 
10.85 + 0.20 Mev. (89), and assuming that the 0.56 Mev. y-ray follows the 
2.89 Mev. positron (42) in the decay of Cu®, the point at 34 neutrons drops and 
a smooth curve results. The nickel curve then has a shape typical of that region 
and no shell effect appears.* 

We have previously noted (9) that the most stable nuclide of all, in terms of 
binding energy per nucleon, is Ni®. In the light of Figs. 1 and 2 this stability 
must be attributed first and foremost to its 28-proton configuration, and, second, 
to its position near the center of the normal distribution of the nickel isotopes. 
The stability due to this latter effect must be credited, therefore, to the accident 
of its birth, and not to its possessing 34 neutrons.* 

Fig. 2 gives no indication of a subshell at 34 protons. 

* Note added in proof: We have just learned from Prof. Leon Katz (July. 29, 1953) that he 
has found the Cu®*(y,n) threshold to be 10.62 Mev., in exact agreement with the new value reported 
by Dr. Leo Seren at the 1953 Washington Meeting of the American Physical Society. With this 
datum the Cu®-derived value for the mass of Ni® becomes identical with the directly determined 


one. Thus, there may be evidence in Fig. 1 for a slight extra stability associated with the 34 neutron 
configuration. 
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38 and 40 Nucleons 

The zinc curve in Fig. 1 shows a minimum at 38 and a rise toward 40 neutrons. 
It is unlikely that this is a real effect since, in obtaining the mass of Ga®® from 
Zn® and Zn”, a mass of Zn® which is a little smaller and a mass of Zn7° which is 
a little larger would give values more in keeping with the directly determined 
value. These calculations have been discussed above in Section IV (c). Such 
changes would smooth out the zinc curve and remove the only indication for a 
shell effect in the region of 38 or 40 neutrons. 

An examination of Fig. 2 shows that any effect existing at 38 or 40 protons is 
masked by the strong 50-neutron influence. 


40-50 Neutrons 
In this region the curves of Fig. 1 rise smoothly and no subshells are found or 
expected. This suggests that the 10 gg/2 levels are filling in regularly. 


32 Nucleons 

Since the f7/2-f5;2 levels are widely split, it is quite likely that the 3/2 levels 
fill in before the f5,2. If so, a subshell would exist at 32 nucleons. There is possibly 
in Fig. 2 a vague suggestion of extra stability at 32 protons, but this may dis- 
appear when more accurate mass measurements become available. 

From Fig. 1, 32 neutrons appears to be a particularly unfavored configuration 
as far as stability is concerned. 
(d) Correctness of the 6-Term in the Semiempirical Mass Formula 

In Fig. 1 the even-A isotopes of the even-Z elements are connected by solid 
lines. The position of the odd-A isotopes of these elements (indicated by closed 
circles) with respect to these solid lines is indicative of the correctness of the 
6-term employed by Metropolis and Reitweisner in their calculations (36) based 
on the semiempirical mass formula. One sees, when he examines Fig. 1 with this 
point in mind, that there is no regular error in the 6-term, and that it approxi- 
mates reality fairly well in this region. The deviations which do exist, although 
rather larger than the experimental error, do not appear to be unequivocally 
associated with shell structure. 


(e) Summary 

The only noticeable mass effects which occur in the region 28 < n < 50 and 
28< Z< 40 are those associated with 28 and 50 nucleons. These are properly 
attributed to the wide spin-orbit splitting of nuclear levels with high orbital 
angular momentum. We have found no extra stability which can be associated 
with the filling of the f5,2, Ps/2, and p1/2 subshells. 
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DETERMINATION OF RATES OF CHANGE OF POLARIZABILITY 
FROM RAMAN AND RAYLEIGH INTENSITIES! 


By E. J. Stranspury,? M. F. CRAWrorpD, AND H. L. WELSH 


ABSTRACT 


Mean values of the rate of change of polarizability with respect to internuclear 
distance were determined for several molecules from the ratio of Raman and Ray- 
leigh intensities in the gas. The values obtained are: hydrogen, 1.2 X 107'* cm.?; 
deuterium, 1.1; hydrogen chloride, 1.0; hydrogen bromide, 1.2; nitrogen, 1.6; 
oxygen, 1.4; carbon dioxide (1 vibration), 4.2; methane (» vibration), 4.1. It is 
noteworthy that the values for the partially ionic molecules, hydrogen chloride 
and hydrogen bromide, are nearly the same as for hydrogen and deuterium. 


INTRODUCTION 
In the polarizability theory of the Raman effect, as developed by Placzek 
(11), the intensity of the vibrational Raman band of a diatomic molecule 
depends on invariants of the tensor quantity 


aij) = (da;;/dr)o 


where a;,; is the polarizability of the molecule and 7 the internuclear distance, 
the derivative being taken at the equilibrium position. The mean value, a’, of 
the rate of change of polarizability appears as a molecular constant in other 
semiempirical theories, as, for example, in the recent treatment of induced 


infrared absorption, by van Kranendonk and Bird (13). The evaluation of a’ 
for hydrogen has been the subject of a number of theoretical papers. There- 
fore its experimental determination is of some interest. 

A direct determination of a’ by the measurement of absolute Raman intensities 
is impractical; however, the difficulties are obviated to a large extent by using 
the ratio of the intensities of Raman and Rayleigh scattering in gases. This 
method was used by Bhagavantam (8, 4, 5, 6), who measured the Raman- 
Rayleigh intensity ratio for hydrogen, deuterium, nitrogen, oxygen, and 
carbon dioxide. Since Bhagavantam’s value of a’ for hydrogen is considerably 
lower than the values calculated by wave mechanics (Hirschfelder (8), 
Bell and Long (2)) and the value obtained indirectly from the difference 
in the refractive indices of hydrogen and deuterium (Bell (1)), it seemed worth 
while to undertake a new determination of the Raman—Rayleigh intensity 
ratio for the simpler molecules. In these experiments particular care was used 
to eliminate the effects of stray light which is the main source of error in the 
method. The diatomic molecules investigated were hydrogen, deuterium, 
nitrogen, and oxygen, which have covalent binding, and hydrogen chloride and 
hydrogen bromide, which have partially ionic binding. In addition, the totally 
symmetric (v;) vibrations of carbon dioxide and methane were studied. 

1 Manuscript received June 17, 1958. 
Contribution from the McLennan Laboratory, University of Toronto, Toronto, Ontario. 


2 Holder of scholarships under the National Research Council of Canada, 1950-52. Present 
address: Bell Telephone Laboratories, Inc., Murray Hill, New Jersey, U.S.A. 
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The ratio of Raman and Rayleigh intensities for transverse observation is 
given by 





Tram _ (¥ — 0)! (45a"* + 13y"*) _h 





Tein y" (45a° + 137°) 82'yw,’ 
where y is the frequency of the exciting radiation, w, the vibrational Raman 
shift, a and y? are the mean value and anisotropy of the polarizability tensor, 
a’ and 7’? the corresponding invariants of the rate of change tensor, u is the 
reduced mass of the molecule, and w, its frequency of vibration for infinitesimal 
amplitudes. The mean value of the polarizability can be determined from the 
refractive index, 7, of the gas from the relation 









n— 1 = 2rNa, 






where N is the number of molecules per cubic centimeter. The ratios, a’?/y’? 


and a’/y*?, can be calculated from the depolarization factors, py.ia, and 






by means of the formulae, 






Prayleigh’ 


92 





Oy 6y* 
PRaman — Aber = £. Fy and PRayleigh = 45a? + 77? . 





T ; if Z . ¢ > ae : ™ | 1| ~e > 
rhus, if x, PRaman’ and Prayleigh; and Leet Trayletgh are measured, jae] Can be 






determined. 

In the formula given above the intensities of the rotational wings of both 
the Raman and Rayleigh bands are included. If the rotational lines are widely 
spaced it is advantageous to measure the ratio of the intensities of the Q branches 






alone. The formula then becomes 


Raman ed (p= wo)' (450”* = 13xy"*) is 
Tay tetgh a vy (45a° + 13x7°) Sm uw,’ 





where x, the fraction of the anisotropic intensity in the Q branch, can be cal- 
culated from the probabilities for rotational transitions (Manneback (10), 
Placzek and Teller (12)). 

The Raman intensities corresponding to the totally symmetric vibrations 
of the linear symmetric molecule XY», and the tetrahedral molecule XY, can 
also be expressed in terms of the invariants of a tensor a;,’, where the derivatives 
are taken only with respect to 7, the length of the X—-Y bond. The constant, 
wu, in the above formula is 2m for XY» and 4m for XY4, where m is the mass 


Vy 














of the atom. 





EXPERIMENTAL 






The Raman tube for gases, shown schematically in Fig. 1, reduced the stray 
light to a very low value. It was made from standard wall Pyrex tubing and 
was blackened on the outside except for the window and the illuminated section. 
Diaphragms of black glass, bevelled on the edges, were positioned in the front 
section to prevent light from the lamp from reaching the window directly or 
by reflection from the walls of the tube. The light trap, made by drawing out 
the rear end of the tube into a smooth cone, was well removed from the short 
illuminated section so that no direct light from the lamp could reach it. When 
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Fic. 1. The Raman tube used in the measurement of the ratio of Raman and Rayleigh 
intensities in gases. A, mercury lamp; B, water jacket; C, plane window; D, gas inlet. 


the inside of the light trap was scrupulously clean, the stray light from the 
evacuated Raman tube was about one-tenth as intense as the Rayleigh scatter- 
ing from oxygen at a pressure of seven atmospheres. I]lumination was provided 
by a helical mercury arc with internal water cooling in the electrode pools (14). 
All gases except deuterium were of commercial purity and were introduced 
into the evacuated Raman tube through a tightly packed cotton wool plug 
which served to remove dust particles. Deuterium was prepared by the action 
of calcium metal on heavy water. 

A condensing lens was used to image the front and back diaphragms of the 
Raman tube at the collimator lens and the slit of the spectrograph, respectively. 
Since high dispersion was not required and short exposure times were desirable, 
a glass spectrograph with an F/4 camera lens and a reciprocal linear dispersion 
of 225 cm.—! per mm. at (4358 was used. 

To measure the Raman—Rayleigh intensity ratio, the Raman scattering, 
the Rayleigh scattering, and the stray light from the evacuated Raman tube 
were photographed in separate exposures on the same plate. Stepped slit spectra 
of a Kipp and Zonen standard lamp were used for the intensity calibration 
of the plates. In general five spectra were photographed on each plate. 

(a) The Raman spectrum. The spectral slit width was from 10 to 40 cm.~! 
and the exposure time for Eastman 103aO spectroscopic plates was from one 
to eight hours depending on the gas. 

(6) The Rayleigh spectrum. A group of Wratten filters, whose transmissions 
were measured with a photoelectric spectrograph, was used to reduce the 
intensity so that the exposure time was about one hour. 

(c) The stray light spectrum. The intensity was adjusted by a group of 
Wratten filters to give the same exposure time as for the Rayleigh spectrum. 

(d) The intensity calibration for the Raman spectrum. The exposure time 
was the same as (a). 

(e) The intensity calibration for the Rayleigh spectrum. It was recorded 
with the same exposure time as (6) and (c). 

The same slit width was used for each of the first three exposures, and the 
current in the mercury arcs was maintained at 13 + 0.2 amp. The correct 
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densities were obtained in exposures (d) and (e) by adjusting the current in 
the standard lamp. Since the spectral distribution of intensity and its variation 
with current were known for the standard lamp, these exposures could be used 
to obtain the intensities of the Raman and Rayleigh lines on the same scale. 

Densities were measured with a recording microphotometer. The profiles 
of the Raman, Rayleigh, and stray light lines were plotted on arbitrary intensity 
and frequency scales. The areas under the profiles were then corrected by 
factors determined from the transmissions of the Wratten filters, the intensities 
of the standard lamp, and the dispersion of the spectrograph, to give the relative 
intensities of the lines. Finally, the stray light intensity was subtracted from the 
observed Rayleigh intensity to give Jp, reign. 

The random errors inherent in the photographic measurement of intensities 
were lessened by averaging the results from several plates for each gas. The 
Raman-Rayleigh intensity ratios determined from individual plates showed 
an average deviation from the mean in the range +3% to +10%. Three values 
of the ratio for nitrogen obtained with \4047 as exciting line agreed well with 
five values obtained with \4358, when the frequency dependent factors were 
removed. The value of a’ obtained for deuterium (1.1 &K 10~!® cm.?) was in 
fair agreement with the value for hydrogen (1.2 X 10~'*cm.?), the small dis- 
crepancy probably being caused by the presence of impurities in the deuterium. 
The consistency shown by these results indicates that the methods of hetero- 
chromatic photometry used introduced no systematic error. 

Some error in the measurement of the intensity ratio might be expected to 
arise from the light scattered by dust particles suspended in the gas. However, 
the measured Rayleigh intensities for the various gases, corrected for the 
stray light from the evacuated tube, and reduced to the same conditions of 
slit width and gas pressure, were found to be closely proportional to a’. It 
was therefore concluded that there was no significant contribution to the 
Rayleigh intensities from light scattered by dust particles. 

Since the illumination of the Raman tube was cylindrically symmetrical, 
the light scattered in the direction of observation was unpolarized and the 
measured intensity ratios were independent of the transmission characteristics 
of the spectrograph for polarized light. The conditions for transverse observa- 
tion were not strictly satisfied. since the angle, 8, between the incident rays and 
the direction of observation ha all values from about 30° to 150°. For the angle 
6 the Raman-Rayleigh intensity ratio is : 


_ = we)! (4a!* + 7") Is + Ip c0s"8) + 6y"Iy sin’ _hk 


y' (45a + 7y°)(1, + I, cos’) + 67, sin’ = 8 rua,’ 
where J, and J, are the components of the incident intensity with electric 
vector perpendicular and parallel, respectively, to the plane containing the 
incident ray and the direction of observation. 7, and J, can be calculated by 
applying the Fresnel boundary conditions at the surfaces of the Raman tube 
and the water jacket, taking into account the inverse square law of intensity. 
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The observed intensity ratio corresponds to the value of Zp anan(9)/Lpayietgn 9) 
averaged over the possible values of @. A calculation made for a slightly idealized 
Raman source showed that the effect of the oblique rays was not at all large; 
for the extreme case of completely depolarized Raman scattering and completely 
polarized Rayleigh scattering the Raman—Rayleigh ratio is depressed by only 
15%. The Raman bands and lines (Q branches) used in the determination of 
a’ were fairly highly polarized and the error caused by the oblique rays was not 
greater than 5% in the most unfavorable case. Since this is smaller than the 
experimental error, no correction for the effect of oblique rays was made. 

The data used in the reduction of the experimental results are given in Table 
except for hydrogen and deuterium are those quoted 

TABLE I 

DATA USED IN THE REDUCTION OF THE RESULTS 


I. The values of pra vieien 


7, 24 
Gas We a X 10 PRayleigh PRaman 
(cm.~!) (cm.3)* 


He 1395 0.84 0.012 14 
D2 3118 0.84 0.012 14 
HCl 2990 2.69 0.007 4 

HBr 2650 3.71 0.008 4 

No 2360 1.79 0.036 19 
Oz 1580 1.63 0.065 .26 
CO, 1340 2.69 0.097 20 
CH, 2914 2.67 0.000 .00 


*Calculated from refractive indices for \4358 quoted in the International Critical Tables. 


by Bhagavantam (6). For hydrogen and deuterium p,,,,,,., has not been 
accurately determined; results of direct measurements, ranging from 0.009 
to 0.038, have been given in the literature for hydrogen, and no results are 
available for deuterium. Although this constant has little effect on the derived 
values Of @’, Pravieyn Was determined for these gases, as a matter of interest 
per se, by measuring the relative intensities of the unshifted Rayleigh line and 
the rotational Raman lines. These experiments gave ppg, = 9.014 for 
hydrogen, and 0.010 for deuterium. The mean of these values was used in the 
determination of a’ for both gases. For all the gases except deuterium, hydrogen 
chloride, and hydrogen bromide, px,,,,, has been measured by Cabannes and 
Rousset (7).* The value of this constant for deuterium was assumed to be the 
same as that for hydrogen. The values for hydrogen chloride and hydrogen 
bromide, which have not been measured previously, were determined from the 
intensity ratios of the S branch lines and the Q branch of the Raman bands. 
Since a’ was determined mainly from the isotropic scattering, it is not very 
sensitive to the accuracy Of ppamnan 
RESULTS 

The results are summarized in Table II. The wave length of the mercury 
radiation used to excite the Raman effect is given in the second column; in the 

* Cabannes and Rousset measured the depolarization of the Q branch in hydrogen as 0.05 and 


calculated ppaman Oy assuming that x (the fraction of the anisotropic scattering contained in the Q 
branch) is equal to 0.25. However, the correct value of x is 0.32. 
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TABLE II 
N AND RAYLEIGH INTENSITIES AND DERIVED CONSTANTS 


RATIOS OF RAMA 





| 


Gas | Exciting 4 la’ xX 10'6(cm.2)} 70) e) LD? 51 
| line (A) 7b. ———— x 108 | | ©! iy te 
Rayleigh (¥ — @o) ; 
as Sr eS a eee ee 
Be 4047 18 1.2 16 
D, 4047 /4358 9 1.1 
HCl | 4047/4358 E2 1.0 15 
HBr 4047/4358 1.0 Rea 26 
No | 4047/4358 0.9 1.6 4.0 
Oz 4358 Ro 1.4 4.4 
CO2( 1) 4358 1.3 | 4.2 | 11 
CH4a(v) | J 4358 3.4 | 4.1 24 


| (4047/4358 


| 
| 
| 
| 


cases designated 4047/4358 the Rayleigh line at 4358A was compared with 
the more conveniently placed Raman line excited by 4047A and the intensity 
ratio was corrected by the factor, 1.65, determined from experiments with 
nitrogen for which both Raman lines could be measured easily. The observed 
intensity ratios are given in the form, (Zpaman/Rayteten) (v4/(v — wo)*), which 
is independent of the frequency of the exciting light. For carbon dioxide and 
methane the total intensities of the bands were measured; for the other gases 
the intensities refer to the Q branches only. The ratios are higher than those of 
Bhagavantam (3, 4, 5) by a factor of about two in most cases; however, Bhaga- 
vantam (6) states that his results should be considered as giving only the order 
of magnitude correctly. 

In addition to the derived values of a’, the values of (Ola\ 1)? = (a’? + 13y'?/45) 
X (h/827uw,) are given. The latter quantity, the square of the matrix element of 
the polarizability for the 1—0 vibrational transition averaged over all orienta- 
tions of the molecule, is proportional to the total intensity of the Raman band. 


The values of (0| a! 1)?, which include all of the anisotropic intensity, can be 
considerably in error if pp. ijan 1S NOt accurate; this is particularly true of the 
more highly depolarized bands. 

The results given in Table II show that there is some correlation between 
ae’ and the multiplicity of the chemical bond; for the four compounds with a 
single bond the values of a’ lie between 1.0 and 1.2 X 10-'*cm.? For oxygen 
with a double bond a’ is 1.4 X 10~'* and for nitrogen with a triple bond a’ is 
1.6 X 10-'®. For polyatomic molecules a’ increases with the number of bonds; 
thus, for carbon dioxide with two double bonds a’ is 4.2 X 10-!* and for methane 
with four single bonds it is 4.1 X 10-'®. The degree of covalence of the bond 
appears to have little effect on a’ since the values for the partially ionic molecules, 
hydrogen chloride and hydrogen bromide, are as large as that for hydrogen. 
It has usually been assumed that the Raman scattering from ionic compounds 
should be very weak. 

Although the sign of a’ cannot be determined from experiment, there is no 
doubt that the positive sign should be chosen. Quantum mechanical calculations 
show that a’ is positive for hydrogen, and the isotope effect discussed by Bell 
(1) gives positive values for hydrogen, hydrogen chloride, hydrogen bromide, 
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and methane. The fact that the polarizability of the helium atom, the limiting 
case of the hydrogen molecule as r — 0, is less than the polarizability of the 
hydrogen molecule indicates that a’ is positive for hydrogen. Also, the polariza- 
bilities of neon, argon, and krypton are less than the polarizabilities of methane, 
hydrogen chloride, and hydrogen bromide, respectively. It is therefore fairly 
certain that in general the polarizability increases with internuclear distance. 

For four molecules studied in this investigation the value of a’ has been 
calculated by Bell (1), using an ingenious method based on the difference in the 
refractive indices of isotopic compounds. If a@ and a; are the polarizabilities 
obtained from the refractive indices of the isotopic molecules, a’ is equal to 


(a — a;)/(& — #;), where = (r — r,), the average departure from the equili- 
brium position in the v = 0 state, is calculated quantum mechanically from the 
cubic approximation to the Morse potential curve. Bell’s results are compared 
with those obtained in the present investigation in Table III. The agreement is 


TABLE III 


Value of a’ & 10!* (cm.?) 
From Raman-—Rayleigh From refractive indices 
intensity ratios | of isotopic compounds (1) 


Gas 


H» 1.68 


1.07 
7.71 


HBr 
CH, 


} 
| | 
HCI | ; 1.23 
‘ 
| ° | 


probably within the uncertainties in Bell’s values which are sensitively dependent 
on the accuracy of the potential function used. 

The value of a’ for hydrogen is of particular interest since wave mechanical 
calculations of this constant have been made. The results of these calculations 
are compared in Table IV with the experimental results. The theoretical values 


TABLE IV 
VALUES OF a’ FOR HYDROGEN 


Method a’ X 10'6 (cm.?) 


Raman and Rayleigh intensities 
Bhagavantam (4) 
This investigation 

Refractive indices of hydrogen and deuterium 
Bell (1) 

Wave mechanical calculation 
Hirschfelder (8) 0.89-1.00 
Bell and Long (2) 1.49 
Ishiguro, Arai, Mizushima, and Kotani (9) 1.411 


of a’ have a considerable spread about the present experimental result. The 
most recent paper is by Ishiguro, Arai, Mizushima, Kotani (9) who have taken 
into account the mechanical and electrical anharmonicities of the molecule. 
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A comparison with the experimental results is best made through the Raman- 
Rayleigh intensity ratio; their theoretical value of ¥*T roman! (v — oe) Tg actais 
is 33.6 X 10-*; the experimental value is 18 X 10-*. It is unlikely that the 
experimental result is in error by as much as 90%. 
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INTERPRETATION OF THE FLUCTUATING ECHO 
FROM RANDOMLY DISTRIBUTED SCATTERERS. PART [I 


By J. S. MARSHALL AND WALTER HITSCHFELD 


ABSTRACT 


The radar echo from randomly distributed scatterers fluctuates in intensity 
about a mean equal to the sum of the intensities contributed by each particle. 
If k independent intensity values are averaged, the r.m.s. deviation of the fluc- 
tuation is reduced by a factor yk. Counting the fraction of signals falling above 
one or more thresholds is a practically simpler, though somewhat less rapid, 
alternative to averaging. As a further alternative, one can superimpose many 
independent sweeps of amplitude or signal level on an A-scope, and so observe 
the most probable value. 

In range interval /, a single sweep contains 2//h independent data, where h is 
the pulse length. The intensity averaged under the continuous trace is slightly 
more accurate than the average of the independent data. 

Pulses sent through the same scatterers yield identical returns until the 
scatterers have had time to reshuffle themselves. Successive pulses give indepen- 
dent echoes, however, if their frequencies differ by about 1/7, 7 being the pulse 
duration. Echoes from the same pulse received at separate antennas are similarly 
independent. 

Several (k) independent echo values are helpful in measuring steady echoes 
against a noise background. For fluctuating echoes, k must be large, regardless of 
the relative strengths of echo and background. Large values of & in addition 
materially enhance the detectability of fluctuating echoes weak compared to the 
noise. 


1. INTRODUCTION 


The radar echo from a random array of scatterers is ‘‘incoherent’’ and so the 
average intensity is well known to be the sum of the intensities from the individual 
scatterers. The random fluctuations of the signal about its average value provide 
an important limitation to the accuracy of the information obtainable. 

The purpose of this paper is to investigate these fluctuations theoretically 
(and by analogue computation) with a view to establishing the extent of this 
limitation, and the improvements that may be effected by averaging or by 
modifying the method of pulsing and scanning. 

The present approach is an elementary one. Mathematical proofs of some of 
the points raised are provided by Wallace (5). 

2. PROBABILITY DJSTRIBUTIONS 

Let a signal y = A sinwt be generated at range r = 0, starting at time ¢ = 0, 
and ending at time rt = //c, where c is the velocity of light. The quantities 7 and 
h are the duration of the emitted pulse in time and its length in space; w = 2», 
where » is the radio frequency. The resulting pulse of electromagnetic radiation 
is beamed by an antenna into a region containing scatterers. Some of the energy 
is scattered back towards the transmitting site, where it may be received by the 
transmitting antenna. 

At any time ¢t > h//c, the receiver will be receiving signals from all scatterers 
in the beam between distances 3ct and }(ct — h). These signals will be distributed 

1 Manuscript received June 4, 1958. 
Contribution from the Department of Physics, McGill University, Montreal, Quebec. 
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randomly as to phase. Thus the region from which signals are being received at a 
particular instant is a section of the beam, of length along the beam //2, and of 
cross-sectional area which we will denote by B. Every point in this region has a 
particular phase; i.e. the position of the scattering particle determines the 
phase of the signal received from it. The volume B//2 of the region is uniformly 
distributed over all possible phases. This region will be referred to as the ‘“‘con- 
tributing region”’ for the time of reception ¢, or for the corresponding recorded 
range r = }ic (see Fig. 1). 


oF puL sé 
TION 
posi iME t 





Fic. 1. Positions of pulse and contributing region at time ¢. 


Consider that fraction of the active volume giving signals of phase between 

6 and 6 + Aé. It consists of a series of cross-sectional slices of the beam, having 
- Ad Bh : ; : ; . 

a total volume of 4 which we will denote by V9. Let us treat it as a region 


2e 2 


all-in-phase. 

In the case of randomly distributed particles, the number in volume V at a 
given instant is Vi + G(Vn), where “is the average number in unit volume and 
G(Vn) is a Gaussian deviate whose mean square value is V#. The probability 
distribution of G is then 


[1] P(G) = ( . e Oe Gg? = Via, 


P(G) dG being the probability that G falls into the range of values G to G + dG. 

Since the particles in this volume are scatterers all in the same phase and each 
giving amplitude a, the total amplitude is Via + G(Vfa’). If various sizes of 
scatterers are involved, this becomes 


Via, +. Viiods 2 + G(Viyay? aa Viva? oS Agate 
In either case, let A, denote the amplitude of the signal from the region under 
‘consideration; then 
Ap = Veda + G(Ve >a’), 
where the summations are over unit volume, G as before is a Gaussian deviate 


er . ; . Aé Bh 
whose mean square value is ¢? = Vs >-a?, and V» = ele ee 
oT « 
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For use in summation over the various phases rectangular coordinates may be 
taken. The X-component of A¢ is 


Xo = [Ve Da + G(Ve da?) cos 4]. 


The X-component of the sum of all amplitude elements, and so of the resultant 


phasor, is 


X => Ve(doa) cos 0 + G[ > Vo( ><a”) cos’ 6] = 0 + G(4Bh>_a’), 


obtained on the assumption that }va is independent of 6, i.e. that the mean 
density of the scatterers is the same throughout the contributing region. 

It follows at once that X and Y similarly are Gaussian deviates of mean square 
values 


[2] X? = 1Bhoa’ and Y’ = 1Bh 0a’. 
Since A? = X?+ Y?, therefore 
(3] A? = X?+ Y? = 1Bh> a’. 


Now A? and a? are proportional to the total and individual intensities, respec- 
tively, and so equation 3 states that the average intensity is the sum of the 
intensities from the individual scatterers. 

The foregoing addition of the random phasors being received at a given instant 
is Rayleigh’s (4) problem of the ‘‘random walk’’. While the answer is well 
known, it has been developed here using a form of grouping most convenient to 
later arguments in the present work. 

The probability distributions of X and Y are Gaussian: 


1 —X2/202 


t PX) a 
[4] (X) ate 


and similarly for Y, where the mean square deviation o? of X and Y is given by 
~ 2 2 
[5] o = 4Bh).a’. 


The summation is over unit volume; multiplication by }Bhk makes it a sum- 
mation over the contributing region, i.e. the region heard from at any one instant. 

Using equation 3, 20? may be replaced by A2, which henceferth wil! always be 
done, since A?, the mean square value of the received amplitude and hence a 
number proportional to the intensity, is a most important and easily visualized 
quantity. In subsequent work it will be convenient to disregard the factor of 
proportionality between A? and the mean intensity; so that A? will be the symbol 
both for mean-square amplitude and mean intensity. 

The probability of the vector ending in an element of area dX dY at position 
X, Y in the phasor plane is 


. ; 2 os 1 —~(x2+¥?) /A2 , 
[6] P(X) P(Y) dX dY¥ = — e& @'t’”/™ ax dV. 

TA 
This is a dome of which all cross sections are Gaussian (Fig. 2). The region about 
the origin is the area on which the phasor is most likely to terminate, but zero 








MARSHALL AND HITSCHFELD: INTERPRETATION OF FLUCTUATING ECHO. I 965 


o«. 





Fic. 2. Perspective sketch of the Gaussian dome P(X, Y) (equation 6), with A, X, and Y 
measured in units of o, the r.m.s. value of X and Y. 
is not the most likely value for the amplitude, which is the distance from the 
origin regardless of phase. 
Probability Distributions of Amplitude, Intensity, and Intensity Level 


By rewriting equation 6 in plane-polar coordinates, and integrating over all 
values of 6, the probability distribution for the amplitude A 


[7] P(A) dA = (24/A’) €*°'" dA 
is obtained. 

It is the quantity A?, the intensity of the received signal, that has been 
related in equation 3 to the scattering cross section of the particles }-a? (and so, 
for Rayleigh scattering, to };D*, where D is particle diameter). Thus the prob- 
ability distribution of A®, the instantaneous value of the intensity, is of special 
importance. 


Using P(A?) dA? = P(A) dA, one finds 
[8] P(A’) dA? = (1/4°) €4°* dA’. 


This distribution is a simple exponential; its maximum value is at the origin. 

The range of values of intensity is often very great. For this reason, displays, 
measurements, and analyses in terms of intensity level (in db. = 10 logyA?/A 9?) 
are often preferred. In such an event, the distribution in the logarithm 


19] P(log A*)d log A’ = | exp ) log A” —- ‘ oe ania| 


MA PL OM (C fdloga, 


where 1/M = Inl10 = 2.30258, is required. 
Some properties of the above distribution functions are listed in Table I. 
Sketches of equations 4, 8, 7, and 9 appear on the left of the table. 
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3. PROBABILITY DISTRIBUTIONS FOR SEVERAL INDEPENDENT DATA 






The widths of the above distribution functions suggest that little significance 
can be attached to single echo measurements. Several data are required for the 
reduction of the fluctuations. Only when the fluctuations are small can A?, the 
mean value of the signal intensity, be determined from the measurements with 






any degree of confidence. 






Averaging A* 
It is required to determine P(J,), the distribution of J;, the average of k 





independent intensity values. Thus 






[10] 











The required distribution is the product 


(11] P(A;*) P(As?) ... P(A,?2) dAy? dA;?...dA,;? 


summed over all values of the A®’s which satisfy condition [10]. Thus 








kJ%/A* . . 2 
[12] P(Jp) dJy = ey es eee 

(A*)* 

ay es 





(k — | integration signs). 









where the integration is carried out subject to )> A,’ = kJj. 





The problem can be solved in an elementary way! by the method of induction. 
The integrations can be carried out for k = 2 and k = 3 by changing [12] into 
polar coordinates, and from P(J2) and P(J/;) results for other low values of k can 
be obtained. An inspection of a few such expressions allows the form for any 






value of k to be inferred to be 
: k* ye 
[13] P(J,) dJ, = aia ee eo ae 












This general result is the correct one, because (i) it reduces to the proper forms 
for k = 2 and 8 respectively, and (ii) it can easily be shown that an increase of 






k by 1 leads to the corresponding expression for J;+1. 

Graphs of the distribution P(J,) for k = 1, 2, 10, and 100 are shown in Fig. 3. 
It will be noted that the distributions become increasingly more peaked (at 
A*(1 — 1/k)) and narrower (standard deviation = A?2/+./k) as k increases. 
This is brought out more usefully by the curves of Fig. 5. Here, using Pearson’s 
Tables of the Incomplete ['-Function, the limits below which 2.5, 10, 25, 75, 
90, and 97.5% of the averages of k independent measurements of A? may be 
‘expected to fall have been plotted. The vertical separations of corresponding 
members of this family are the 95%, 80%, and 50% confidence limits. It will be 
seen, for instance, that 10% of the averages of 50 independent measurements fall 
below 0.80A2, while 90% come below 1.175A?. It follows then that the 80% 













‘For a more general method of formulating the problem, see Wallace (4). 
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Vk 


Fic. 3. The probability distribution of Jz, the average of k independent values of signal 
intensity, when the true mean intensity is A? (see equation 13). 


confidence limits are approximately A? + 0.24. When the confidence limits are 
plotted against k~“’”, they appear as almost straight lines converging towards a 
point at kR-GV/) = 0. 

The above solution may be formulated usefully in an inverse manner (as 
suggested by Wallace (5)) by inquiring into the precision to which A? may be 
determined when an observed value of J; is available. The relative probability 
that the true average has a value A?, given a value of J;, is simply (from equation 
13) (A2)~*e-*4*/4*: the distribution in A? follows by normalization: 
eo tysla® (A*)* 


> 


-gA’. 


[14a] Q,(A*) dA? = —— 
J ebay A? (A?)* 


0 


The normalizing integral in the denominator is finite only if k > 2; that is to say, 
from a single intensity measurement no information about A? may be obtained; 
or, to put it more precisely, for a given measurement of J; there is always an 
overwhelming probability that it came from a very large A?. The fact that we 
do not really admit this overwhelming probability is based on considerations 
quite extraneous to the measurement made, e.g. a general knowledge of the 
nature of the disturbance from which the echo comes. 
If k > 2, equation [14a] becomes 


ows c= Ae k-1 se are 
[149] 0, (4?) dA? = a oa etl Gat 
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This distribution is plotted for k = 2, 3, 10, and 100 in Fig. 4 which should be 
compared with Fig. 3. The new curves are in general broader, though as k 
increases, the two distributions become more and more alike. A notable differ- 
ence is that the most probable value of A? equals J, for all values of k. Also, the 


Fic. 4. The probability distribution of A?, the mean signal intensity, in terms of J,, an 
average of k independent values of signal intensity (equation 140). 


% 40 60 80 


Fic. 5. Width of the distribution P(J;), expressed in units of A? (along the ordinate axis) 
The 90% curve, for instance, indicates that 90% of the averages of 50 independent values of 


signal intensity (A?) may be expected to fall below 1.175A?. 
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60 


Fic. 6. Width of the distribution Q;(A2), expressed in units of J, (along the ordinate axis). 
The 90% curve indicates, for instance, that there is a 90% chance that A? (the mean of very 
many measurements) falls below 1.225 Js», Js) being an average of 50 independent intensity 
values. 


nth moment of the Q,(A2) distribution exists only for k > (mn + 2). Limits of 
confidence for Q,(A2), similar to those for P(J;,), are shown in Fig. 6. Some 
properties of both distributions are listed in Table I. 


Averaging A and log A” 

The probability distributions of averages of k independent values of signal 
amplitude A and intensity level log A® cannot be found by the above methods. 
A numerical method was therefore adopted. One thousand random numbers for 
X and another thousand for Y (all taken from Wold (6) ) were combined to form 
sets of one thousand independent values of A and of log A”, obeying, respectively, 
the distributions P(A) and P(log A’). Confidence limits, describing the con- 
vergence of the averages of k independent values of A and of log A?, were then 
drawn, as functions of , ina manner similar to Figs. 4 and 6. These limits, when 
expressed in comparable terms, were found to be very close to those for P(J;,); 
for large values of k, they can in fact be shown to be identical. 

Probability of a Signal Being Above a Given Threshold 

While averaging is the most obvious and theoretically the most satisfactory 
way of dealing with several independent data, the-practical problems involved 
are considerable. Extensive and costly equipment of at best doubtful performance 
is required in the storing and averaging, if a reasonably large number of data are 
required. An alternative method— the classification of the signals and the deter- 
mination of the fraction falling above a given threshold—can be realized by 


much simpler and more reliable equipment, as only on-off type of response is 
required in the circuits. The theory of this method of analyzing the data will 
now be presented, and the resulting reduction of the fluctuations will be com- 


pared with that obtained by straight averaging. 
The probability that A* has a value equal to or greater than C? (a given thresh- 
old) is (using equation [9]) 
4° jar 
[15] —;—dA® = ¢ ©! 4° 
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(The probability for A > C and for log A? > log C? is, of course, the same.) 

If k independent values of A* are examined, let p be the number of them 
falling above the threshold C?. This number p has an a priori probability 
p = ke~©*/4* and follows the binomial distribution: 


= ———— eo e)'( z 7 
[16] P,(p) = (p)! (R — pb)! (? : b 


(17] fof mean value jp 


(and mean square deviation p(1 — j/k). 
Let the calculated value of A?, obtained from the observed value of p/k, be 
~ ~ 
denoted by A?. Examine now the fluctuations in A? about the a priori mean 


intensity A. Equation [15] leads to 





2: el ee Ge 
[18] a" Gas * “eae 


where, for brevity, one may put x = p/k and y = A?/C°. 

Using the binomial distribution [16], one can compute the limits x,y in x, 
within which given fractions F of measurements of x may be expected to fall. 
From equation [18] one can then obtain the corresponding limits yy in y, which 
may conveniently be expressed as multiples of 7 = A?/C?. The results can thus 
be displayed (Fig. 7) as plots of yr/%, as functions of p/k or C?/A?. The necessary 
computations were carried out for all integral values of p when k = 10. In the 
cases k = 50 and k = 100, on account of the unwieldy form of the binomial law, 
a Gaussian approximation for P(p) was used. 

From Fig. 7, the smoothing is seen to depend on the proper choice of C’, 
which in turn varies somewhat with k and with the degree of confidence required: 
when k is large (50 or 100), C? should be between 1.2 A? and 1.6 A?, correspond- 
ing to p/k between 0.2 and 0.3. For smaller values of k, lower thresholds are 
indicated. At k = 10, C? should be about equal to A?, giving j/k of about 0.35. 
The confidence limits of Fig. 7 are fairly flat in the vicinity of the optimum 
points; at k = 50 or 100, for instance, the variation in precision over an 8 db. 
range in C? containing the optimum value is barely 20%. It is worth observing 
that this particular optimum range of C? arises as a compromise from the con- 
tradictory indications of equations [18] (sketched in Fig. 8) and [17]. Equation 
[18] permits the evaluation of A?/C? from j/k with the greatest precision where 
the slope of the curve has its flat minimum (at C? = 2A? or p = 0.135k). The 
precision is extremely poor, on the other hand, near p = 0 and particularly near 
p = k, which are precisely the regions where, according to equation 17, minimum 
scatter in p may be expected. 

In Fig. 9, the best smoothing obtainable by the present technique is com- 
pared with the reduction of the fluctuations when k intensity measurements are 
averaged in the usual way. It may be seen that at k = 10 the 95% confidence 
limits are 70% wider with the present technique than when the data are averaged 
in the ordinary way. At k = 100, the widening of these limits is less than about 
30%, which is probably a loss in precision smaller than would be incurred on 
account of experimental error if ordinary averaging were attempted. 
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Fic. 7. Confidence limits for k= 10, k = 50, and & = 100. It may be seen, for instance, 
that when k = 100, one can place A? with 80% confidence within the limits 0.854? and 1.174? 
when C? = 1.5A?, i.e. when p/k is about 0.22. 


An attempt to check the validity of the theory here presented by numerical 
means was made. Values of p/k for k = 100 were calculated from the set of 
1000 independent values of A assembled for use in the previous subsection. 
These values were plotted in Fig. 10 against the theoretical limits of Fig. 7c. 
If allowance is made for the smallness of the sample studied, the numbers may 
be seen to support the theory. 
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Fic. 8. Plot of equations [18], y = = 


—inx 


2.5 and 97.5% 
. 10 and 90% 
25 and 75% 


40 60 60 100 


Fic. 9. Comparison of confidence limits in classifying and counting technique (heavy points) 
and in averaging (curves, copied from Fig. 5). Note, for instance, that when k = 50, the theo- 


retical 80% confidence band in averaging is 0.35A?, and in classifying and counting 0.45A?. 
3 g g 
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k=100 


+ 90% and 10% 
© 75% and 25% 
« 50% 


Fic. 10. Numerical check of the confidence limits of Fig. 


4. DEPENDENT AND INDEPENDENT DATA 


So far this paper has been concerned with independent data. It is now desir- 


able to discuss just exactly what data are independent, how their number may 
be increased, and whether any information can be derived from nonindependent 
ones. 


Achieving Independence in Time 

Consider first the scatterers in a particular ‘‘contributing region” (cf. Fig. 1). 
The specific value of the phasor A depends on the distribution of the particles 
relative to the phase pattern of the beam. If another pulse is sent out after some 
motion of the scatterers has taken place, a change in A in both magnitude and 
phase is to be expected, unless the motions involve exact half wave lengths only 
(i.e. zero phase change), or are small compared to the half wave length. On the 
other hand, this new value of A will never, in principle, be completely indepen- 
dent of the first, because however much time elapses, the configuration of the 
particles at any time will always depend to some extent on the initial distribu- 
tion. Independence is, however, approached gradually. If the scatterers are 
raindrops, experiment indicates that practically complete independence is 
achieved at 10 cm. wave length after time intervals of the order of 0.01 to 
0.02 sec. (announced by Aaron Fleisher at the Radar Weather Conference, 
University of Illinois, October 1-3, 1951). This means that in one second, 
between 50 and 100 independent data may be obtained from one particular 
contributing region. Looking vertically at cloud or rain by means of a 1.25 cm. 
radar, on the other hand, Bartnoff, Paulsen, and Atlas (1) measured times-to- 
independence of the order of 0.00375 sec. 
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The variation of the returned intensity in time intervals shorter than those 
mentioned is related to radio frequency and beam width, and to the frequency 
and amplitude of the relative motion of the scattering particles. This aspect has 
been studied by Hilst, Fleisher, and others, and is reported in a series of Tech- 
nical Reports by the Weather Radar Research group of the Massachusetts 
Institute of Technology. 


Achieving Independence in Range 


As the pulse travels out from the radar, the signal returned fluctuates as new 
scatterers continually come into the contributing region while others leave it. 
These fluctuations may be studied conveniently by considering the contributing 
region to be divided into m elements of length 34 X 1/n, each containing an 
integral number of wave lengths. The signal phasor contributed from each 
element has rectangular components which are Gaussian deviates of mean 
square value (Bh/4n)>ia*. 

As the contributing region proceeds along the beam with speed half that of 
light, new units are included continually, and old ones dropped out. Sample 
behavior of A as the contributing region progresses outward has been computed 
from tables of randomly sequenced deviates (Wold (6) ) by building up the 
components X and Y¥ as the algebraic sums of 100 (for 2) deviates. The addition 


-Y 


Fic. 11. Typical hodograph (of half-pulse length PL 18), with contributing region divided 
into 100 slices. Ags is a typical instantaneous amplitude vector. 
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of one deviate to X and one to Y, and the removal of the 100th-preceding 
deviates, gives new values for X and Y, and so for A, representing the condition 
existing when the contributing region has moved by (1/100)th its own length. 
After 100 such changes, one has a whole new set of deviates, an independent 
value for A. 

An actual case, developed on a 100-element basis, is shown in Fig. 11. Each 
phasor element (the straight line joining successive points on the hodograph) 
is independent of all others except the hundredth-preceding and the hundredth- 
succeeding, with each of which it shares one of its two vector terms. Thus while 
the first and last points on the hodograph for length 4/2 are the only two inde- 
pendent points (i.e. phasors starting from the origin) on the hodograph, every 
single element (line joining adjacent points) is independent of every one of the 
other 99 elements appearing in this hodograph. When a continuous hodograph 
for several successive lengths 4/2 is prepared, the elements become half-indepen- 
dent, each one being related to the hundredth-preceding and succeeding elements. 

It is rather surprising to find that the probability distribution of each of the 
phasor elements is the same as that for the vector sum of 100 of them. When we 
add them vectorially, the intensity (corresponding to A?) increases proportion- 
ately, on an average, to their number, but (if, as for the moment, background 
noise is neglected) the accuracy of the information conveyed remains constant. 
If the individual terms could be squared and their squares added, the accuracy 
would increase from that for k = 1 to that for k = 100. 

The variation of amplitude and phase with range has been plotted (Fig. 12) 
for 30 half-pulse lengths. The plot of A is effectively the picture resulting on an 
A-scope type display, if the band-width of the amplifiers is approximately 
100¢/h and all circuit elements respond linearly. The variation of A? with range 
is shown in Fig. 13. The pre’ vility distribution in this case reaches a maximum 
at A? = 0 (equation 8) _4 so of course the trace spends most of its time close 
to the axis A? = 0. 

A characteristic of the graphs of Fig. 12 worth special mention is that the 
fluctuations in @ become small when A is large. The explanation is that the 
fluctuations in the increment AA from one amplitude to the next are relatively 
small, as is illustrated by the general appearance of the hodograph of Fig. 11. 
When A is large, the effect of the vector increment on the phase angle is therefore 
small, and vice versa. It may be noted in passing that if the average value of A 
over a length, say h/2, were corrected according to the amount of fluctuation 
in 8, a much improved measure of A would be obtained. 

An assessment was made of the value of the information conveyed by incom- 
pletely independent signals, such as make up the curves of Figs. 12 and 13. 
Mean values of A® (or A) obtained by averaging under the full curves were 
compared with the means of the instantaneous values separated in range by 
h/2, i.e. adjacent independent data. On the whole, the results indicate a slightly 
more rapid convergence to the true average when the dependent values are 
included along with the independent ones. This same conclusion is further borne 
out by observations of the standard deviation. The theoretical value of o (for 
infinite samples), for independent values of A, is A?; for averages of 100 inde- 
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pendent values of A? it is A?/+/100 = 0.10A®. For the 30 half-pulse lengths 
under study the deviations turned out to be 1.008 A? (for the independent values 
of A*) and 0.85 A? (for the 30 averages of 100 values of A*). This indicates again 
that an average of 100 such semidependent values is a somewhat better clue to 
true mean (A?) than any independent value alone, but is, of course, not nearly 
so good as an average of 100 independent values. A more quantitative discussion 
of this problem is given by Wallace (5), who considers a continuous smooth 
transition from one independent point to the next (corresponding to our 100 
intermediate steps). The gain in accuracy of information again appears to be 
disappointingly small. The number of independent data available from a length 
along the beam k X 3h isk and the improvement obtainable by averaging a larger 
number of data is never great. 

It has been seen from the hodograph that the signal variations in a range 
interval h/2 contain the equivalent in information of many independent data. 
It has appeared, however, that averaging under the A? (or the A) plot is not an 
efficient way of bringing it out. It is the successive vector variations in A, which 
we will call AA [= A(t + At) — A(t)], that contain the additional information. 
It may be worth considering briefly the practical problems of obtaining this 
quantity. Quite a direct method, simple in principle, though probably difficult to 
realize, might consist in splitting the signal received by the antenna into two 
parts, to be propagated by two wave guides in parallel. One of the guides is slightly 
longer than the other, so that one half of the signal will be delayed by some time 
interval Aft, equivalent to a range interval $cAt, relative to the other. If the de- 
layed term is subtracted from the other, the resultant is, for a perfectly rect- 
angular pulse, just this AA. The signal is, of course, of low intensity, for the 
method if used at all would be for At < 7. Further, amplifiers would be required 
of band-width 1/ Af. 

It is probably easier to obtain the component of AA in the direction of A(t). 
This component is inferior to AA only in having (because it is just a component) 
an intensity level 3 db. lower. It would be obtained by differentiating the 
“video” or output signal, and full-wave rectifying the derivative to obtain its 
absolute magnitude. Again, amplifiers would be required of band-width 1/ Af, 
but where the AA device was most practicable for At « 7, the dA/dt device 
should be practicable for any At for which the corresponding band-width 1/ Af 
is available. This derivative no longer makes sense, however, when circuits 
preceding the video become saturated. 

Now for any interval At between independent data, a band-width of 1/At is 
sufficient to provide all the independent data. An increase in band-width 
beyond this value adds extra detail which adds little useful information, hardly 
increasing the value of k. It does, however, add background noise of intensity 
proportional to the band-width. The graphs of A and A? in Figs. 12 and 13 
consist of points plotted at intervals of 7/100, corresponding to signals passed 
through amplifiers of band-width 100/7. If we observe or average A? (or A) 
directly, all detail beyond one point per time interval 7 is of little use. Similarly, 
in the case of a receiver, any band-width beyond 1/7 is of little use. If, however, 
we observe AA or its component, or the square of either of these, then added 
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band-width up to 1/At can be utilized to increase k, our measure of useful 
information, efficiently. The increment At is here the time interval over which 
AA is measured. 

In general, suppose that signals from precipitation are required with resolu- 
tion in range to a range interval /. The accuracy of the information can be 
increased by making the half-pulse length 4/2 short compared with /, and aver- 
aging along the range over a travelling interval /. Alternatively, the pulse length 
can be left unchanged, but the increments in the signal magnitude rather than 
the total signal are averaged. In either case, the signal component used is 
weakened by this shortening, and the receiver band-width required is increased. 
One advantage of the short time-constant at the output, not shared by a 
shortened pulse length, is that the former could presumably be made variable 
with range. Thus the sensitivity of the equipment could be made independent 
of range, and the excess power at short range disposed of to provide increased 
precision at short range. 


Achieving Independence by Changing the Frequency 


The above discussion of the fluctuations in time in the return from a given 
contributing region demonstrates that considerable time must be allowed to 
elapse before an independent datum can be obtained. The discussion of the 
fluctuation encountered as the pulse is moving out from the radar shows that an 
appreciable loss in power and a gain in background noise level must be coun- 
tenanced if the number of independent signals along a given bearing is to be 
increased. The number of independent signals from a given contributing region 
can however be multiplied by shifting the frequency of successive pulses, without 
any sacrifice in the signal/noise ratio, and regardless of the time interval between 
the pulses. 

In Section 2 it was noted that any particular range of phase @ to 6+A0 
could be identified with a series of cross-sectional slices of the beam, the volume 
of which was denoted by Ve. In the case of a one microsecond pulse of 3 cm. 
radiation, for example, the series contains 10,000 slices. The deviation from 
average of the magnitude of the phasor element representing that range of phase 
will be proportional to the deviation from average of the scatterers in that series 
of slices. If the wave length of the radiation is changed slightly, so that the 
number of wave lengths in the pulse length is changed by one, then the volume 
Vo is no longer identifiable with a single phase, but instead is distributed over all 
possible phases (Fig. 14). The specific contribution that is made to the resultant 
signal is thus eliminated, even though its chance of contributing, on the average, 
remains unchanged. This applies to each element of volume Vo, and to the corre- 
sponding phasor element Ag. It appears, therefore, that a change of one in the 
number of wave lengths in the pulse length should give an independent value of 
A. 


If the spreading-out of V¢ over 27 (or any integer X 27) makes the new signal 
value completely independent of the old, it would seem reasonable that a 
spreading-out of Vs over any angle should create some degree of independence. 
Wallace (5) gives a complete mathematical analysis of this problem, confirming 
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SLICES ORIGINALLY AT ZERO PHASE 
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SLICES FINALLY AT ZERO PHASE 


0 TT iT 


PHASE CHANGE 


Fic. 14. Showing how the phase corresponding to each slice contributing to A (for 6 = 0°) 
is increased as the number of wave lengths in the pulse length (and so the number of slices) is 
increased by one. For the purpose of the sketch, only eight original slices have been taken to 
make up the volume Vo. 


the above conclusions. He shows in addition that effectively independent returns 
can be obtained by frequency shifting, not only for rectangular, but even for 
somewhat distorted pulses. 

For the purpose of illustrating this effect further, it is sufficient to quote one 
of Wallace’s results (also obtained by Goldstein (3)*) valid for a rectangular 
pulse. The correlation coefficient between the signal intensity A»? received from 
a set of scatterers at frequency v» and the intensity A* obtained from the same 
scatterers at frequency vy + Av is given by: 


sin mrAv\? 
9 == —_— 
[19] p ( ar Av ) ; 


where the correlation coefficient is defined by 
[20] gs Ao A* — (A*)’ 

A* a (A*)? 
Equation [19] is sketched in Fig. 3 of Wallace (5). The correlation between 
A? and A,’ falls off rapidly as the frequency shift approaches 1/7, the reciprocal 
pulse duration. When Av is greater than about 3/47, A* may be regarded as vir- 
tually independent of A,’. 


Achieving Independent Data by Change of Aspect 

Two (or more) independent signals may be obtained simultaneously from the 
scatterers of a given contributing region, if the scattered energy is received at 
two (or more) separate antennas. One of the antennas may also do the trans- 
mitting, so that the method to be discussed involves one ordinary radar set and 
at least one additional receiver. Independence comes about in this procedure 
through the fact that randomly distributed particles present randomly different 
‘ aspects to properly spaced receiving antennas. In Fig. 15, AB represents a thin 
slice of the contributing region from which signals are returned to the radar, all 
at the same phase. Signals arriving from the point A at the separate receiver 
will be changed in phase (delayed) by 2x(r — r’)/X, while those from B will be 


*There is a minor misprint in Goldstein's equation (134). 
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Fic. 15. Two separate antennas receive practically independent signals from the same 
scatterers in the slice AB. 
changed in phase (advanced) by 2r(r — r’’)/X. The phase changes at inter- 
mediate points will be between these limits. At a point C, practically halfway, 
the phase change will be zero. By suitably choosing d, the distance between the 
antennas, it is possible to make the phase changes at A and B equal, respectively, 
to —m and +z, or integral multiples thereof. In such an event, the phases of the 
signals received at the two antennas from each scatterer will differ by effective 
amounts between —z and +7, the exact value of the difference depending on the 
random position of the particle. The difference itself will therefore be random, and 
the two signals will be independent. Using the fact that \ and d are very small 
relative to the range 7, and on the assumption of a uniform and sharply defined 
beam, it is easy to show that for complete independence 
[21] d = phon wy 2 


where @¢ is the beam width. In fact, the boundaries of the beam are not sharp, 
and a more precise theory would be required. 

The above argument follows so closely that given under the heading of fre- 
quency shifting, that it is reasonable to infer a formula similar to [19]: 


> ee apd ay 
p apd /d 
for the correlation coefficient between the two signals. When the beam width is 
defined as the angular distance between half-power points, ¢ is of the order of 
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\/D, where D is the diameter of the antenna. The appropriately refined form of 
equation [21] would then indicate that the smallest value of d for independence 
is of the order of D, or that the two antennas should be contiguous. It would 
follow that while only certain discrete values of d lead to complete independence, 
in practice any separation at all would lead to a useful amount of it. 

An argument very similar to the above would show that if the radar itself 
(i.e. transmitter as well as receiver) were moved sideways, independent signals 
would be obtained from the same scatterers at intervals equal to about one-half 
of those calculated by equation [21]. 

If the beam is not filled with scatterers (as has been assumed above) very much 
the same theory holds, except that @ needs to be redefined as the angle subtended 
by the target at the radar. Two separate receivers and a measured value of the 
correlation coefficient may then be used to determine the size of the target. This 
technique has already been applied in radioastronomy (2). 

When the existence of background noise is taken into account, the apparent 
advantage of the above method of obtaining independent signals is reduced, 
in analogy to what was found above in connection with the suggestion of using 
a shorter pulse length. Clearly, several small antennas, each with its own 
receiver, contribute more independent data than one single large one, but the 
signal intensity carried by each would be small, with an attendant reduction in 


signal/noise ratio. 


5. RADAR DISPLAYS 
The averaging of several independent signals, which has been shown to be an 
essential step in a quantitative evaluation of the return from precipitation, can 
be greatly facilitated by a careful consideration of the display used. In displace- 
ment presentations, or A-scopes, the amplitude A, or some function of A, is 
displayed as a vertical displacement against a horizontal range-sweep. In area- 
displays, the signal function is presented as a brightness modulation of the range- 
sweep, while successive range-sweeps are spread over the area of the display as 
the direction of the beam is rotated to scan an ‘‘area’”’ in space. 


Displacement Modulation: Single Sweep 

A single sweep of an A-scope resembles Fig. 12, or would if the frequency- 
response band-width were about 100/7. The best band-width, however, is the 
reciprocal of the time between independent data. In Fig. 16, therefore, plots 
resembling usual A-scope traces more closely are obtained by drawing straight 
lines through a set of 25 independent points. Such curves have been drawn for 
A?, A, and log A?. 

Alongside the A-scope traces in Fig. 16 have been drawn the probability 
distributions of the signal functions as summarized in Table I, and so of the 
height of the trace. In every case the height of the trace varies widely, and so the 
average height or “height of grass’’ is not the sort of thing that can be determined 
reasonably by a glance at an A-scope presentation. A quantity one can deter- 
mine well, however, from a single sweep, if it is photographed, is the fraction of 
time (over several or many pulse lengths) that the signal spends above any 
specified threshold value. This fraction has been related to the correct value of 
A? in Section 3. The fraction can usually be determined most accurately when 
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P(log A’) log Pllog A*) 


Fic. 16. Top row (left to right): probability distribution of signal intensity (A*), typical 
traces of independent A®-values, and distribution of luminance on ‘‘A?-scope”’ resulting from 
superposition of many independent A?*-traces. Center row, same for signal amplitude, 4; 
bottom row for signal intensity level, log A?. 


its value is of the order of 20%. Any display may be used: A or A? or log A?; 
the height of the optimum threshold will vary most gradually on a logarithmic 
display. Thus with the wide range of values of A? encountered in precipitation 
studies, a logarithmic function should be the most useful one to display. 


Displacement Modulation: Many Sweeps 


The case of several sweeps is more difficult to analyze. When very many sweeps 


are superimposed a new approach suggests itself, however. The superposition 
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may be by long-persistence phosphor, or preferably, by time-exposure photog- 
raphy, where the contribution of each sweep is given equal weight. Such A-scope 
photographs are frequently used and reported; their correct interpretation calls 
for careful consideration. 

With many sweeps superimposed, the detail of the individual lines is lost; 
instead we obtain a continuous gradation of luminance over the display. The 
degree of contrast in the gradation depends on the characteristics of the persis- 
tent phosphor or of the photographic emulsion involved. “‘Saturation”’ is possible 
here, just as it is in electrical circuits, and in most cases it is something to avoid. 
The probability distributions of the signal function are also distributions of 
luminance with height, assuming linear response in the phosphorescent or photo- 
graphic material. Since the processes involved are not usually linear, but approxi- 
mately logarithmic, the logarithm of the probability is more pertinent to the 
description of the appearance of the display or its photograph. Plots of the 
logarithms of the probabilities of the signal functions are also shown in Fig. 16. 

The least helpful display, for many sweeps, is that of A?. Here the log lumin- 
ance decreases linearly with height. The desired quantity A? could be determined 
from the slope of the line, but this-would involve photometry and a knowledge 
of the degree of contrast in the photographic or phosphorescent process. The 
height at which the signal disappears would require a knowledge of the same 
factors, and of threshold sensitivity as well. The A-display is more useful, in 
that the height of maximum luminance and so of most-probable A can be 
determined. A relation between the most probable A(=A,) and A? has been 
given in Table I. Photometry is not needed for this determination: the over-all 
brightness would be reduced until only a narrow region around the maximum 
was perceptible. Unfortunately, the maximum in the curve of log luminance 
against height is rather a broad one. The most useful display of the three would 
be that of log A®. The maximum here is of the same width regardless of its height 
on the screen. Further, as mentioned in connection with single sweeps, the 
motion of the maximum up and down the screen with variations in A? is appro- 
priately more gradual. It is fortuitous and convenient that the most probable 
value of log A® is equal to log A?. The rapid decrease in luminance above the 
height of maximum luminance makes it appear that maximum signal height could 
be read with considerable accuracy. Consideration of the scales involved reveals 
that this is not the case. When A? drops to zero, its logarithm drops to minus 
infinity. The finite lower limit of any logarithmic amplifiers distorts the dis- 
tribution of luminance with height, and for weak signals it must be interpreted 
with care. 

If several independent A’s or A®’s or log A®’s are averaged before presentation 
on a displacement presentation, the resulting distribution in height will be much 
narrower than those shown in Fig. 14, and a firm, narrow trace across the tube 
can be obtained. The only drawback is the need for a memory device of some 
sort before the A-scope. 

Intensity Modulation 

Figs. 12 and 16 show the variations in displacement with range on an ‘“‘A- 

scope’’, or of intensity (luminance) with range along a single sweep on a two- 
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dimensional display. Now suppose that the antenna is rotated slowly, so that 
there are many sweeps on the two-dimensional display for every beam width. 
Then the variation with angle for a given range r will be similar to the variation 
with range for a given direction, except that the former changes discontinuously 
from sweep to sweep. Neglecting motion of the scatterers among themselves, 
each successive signal at the same range but at a slightly changed direction 
incorporates some new volume to the contributing region, discards some old, in 
the same way as an increment in range with a fixed direction (Fig. 17). 


Antenna 

(rotating 4@ 

about this 
axis.) 


Fic. 17. Moving along the beam through range interval 4/2, the half-pulse length, or 
scanning through ¢, the beam width, leads to an independent signal value. 


To picture the random variations in intensity over the surface of a two- 
dimensional display, one must visualize the two-dimensional equivalent of 
Fig. 12. In practice, amplifier band-widths of approximately 1/7 smooth out the 
variations with range until effectively only the essential independent points are 
left. On the other hand, the variations of signal with direction of the beam, at 
fixed range, tend to include many intervening dependent (and so nonessential) 
points. 

Let us consider an example: 

Pulse length 1 usec.(34 = 0.1 mile), 


Beam width @ = 1/2° (=1/2 mile at 56 miles), 


PRF = 1000 sec.~! \ =10 pulses per beam width, 
Rate of sweep = 50° sec.~! § =20 pulses per mile at 56 miles. 


When the degree of independence achieved by relative motion of the scatterers 
as the beam rotates through its own width (0.01 sec.) is neglected, there will be 
nine dependent and relatively useless data at range 7 between successive indepen- 
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dent data. A sample of the random area pattern of A obtained is given in Fig. 18, 
for a region 15 half-pulse lengths by two beam widths, or 1.5 X 1 miles, at 
range 56 miles for the particular case given above. Streaking across the beam is 
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Fic. 18. Typical appearance of a section (15 half-pulse lengths X 2 beam widths) of a 
brightness-modulated display. (Area shown measures 1.5 X 1 miles for the example discussed 
in the text). The numbers represent amplitude values received from element of area resolved 


(h/2 X $). Vertical shading indicates amplitudes exceeding average (A = 12.5); cross shading, 







amplitudes exceeding 1.54 = 18.75. 












apparent, particularly when high threshold values of the signal are applied. 
Similar effects would, of course, be obtained in the case of the A? and the log A? 






pattern. 

If the frequency were varied by one or more megacycles between successive 
pulses, with a saw-toothed pattern over a cycle of 10 pulses, each successive 
value in a column would be independent of its neighbors, just as holds now for 
figures in a row. There would then be 200 independent data per square mile, 
instead of effectively 20 as there are now. Samples are given (for A) in Fig. 19. 

The element of area to which a picture can be resolved is [beam-width X 3h] 
but this is scarcely significant when the randomness of the individual recorded 
data is as great as that displayed in Fig. 19. Without increasing the size of this 
element, it would be useful to average the data within the element. If the 
‘desired elemental area of resolution is greater than beam-width X 3h, the 
averaging might usefully be carried to this larger element, which can have any 
desired shape provided that neither of its dimensions is less than that of the basic 
element. The creation of the larger element could be achieved fairly simply by 
deforming the focus of the cathode ray tube or alternatively by defocusing the 
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Fic. 19. Same as Fig. 18, except that amplitudes in adjoining rows (as well as in the 
columns) are independent. Note how the vertical streaking in the background, so predominant 
in Fig. 18, is largely absent. 
recording camera. The basic element changes its linear dimensions with range, 
of course. Whether the modified element did so would be largely a matter of 
choice. 

If the k-value for the average within the desired element of resolution is too 
low, i.e. leaves too much uncertainty as to the correct signal mean strength, one 
can superimpose successive scans.* Relative motion of the scatterers is almost 
certain to provide effective independence of data for successive scans. Some sort 
of memory device would be required, of course, and the simplest would seem to be 
a photographic record. 

As has been pointed out previously (Section 2) signal intensities or A*’s are 
directly additive, but A2, A, and log A or log A? have fixed relationships among 
them, so that any one of A’, A, or their logarithms or some other function of A? 
can be averaged, so long as the correct factor is introduced. It must be borne in 


mind, of course, that 


A’ # (A)? # log" (log A”). 


Use of A, still more so of log A’, reduces the range of signal strength over which 
the amplifiers and cathode ray tube are required to behave proportionately. 
There is another way, of limited applicability, of reducing the proportional 
response requirement. That is to generate a black-and-white picture, i.e. one of 
extreme contrast, rather than one of continuous gradation. All signals above a 
*A ‘‘sweep”’ is the line traced on the display tube as the pulse moves out from the radar. As the 


antenna rotates once about a vertical axis, or completes one up and down motion ina vertical plane, 
the picture caused by all sweeps is called a ‘‘scan’’. 
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certain threshold are recorded ‘‘white’’, all below that threshold ‘“‘black’’. Then 
the camera or other memory device simply records the amount of white signal 
in a given area in one or any specified number of scans, this being the fraction of 
signal above a given threshold. This fraction has been related to A? in Section 3, 
and its corresponding application to A-scope observations was discussed above. 
The range of signal strength over which this technique works well is only about 
10 db. Amplifiers might possibly be designed, however, in which the signal re- 
corded comes from the first saturated stage of the amplifier, or from the one 
preceding. (Such amplifiers have actually been designed for continuous signals.) 
The fractional technique applied to this type of amplifier should provide equip- 
ment giving reasonable precision over a wider range of signal intensities. 





6. EFFECTS OF BACKGROUND NOISE 





Fluctuating Signals in the Presence of Background Noise 
The echo arriving back at a radar from a random array of scatterers has been 
discussed at length. Its average strength is related to the properties of the distri- 
buted target; but all that is known about its instantaneous value is a probability 
distribution. When the echo enters: the circuits, it combines with a random 
background signal originating in the circuits; this background will have a 
probability distribution of much the same sort, depending on the circuit para- 
meters. On the assumption that the two statistics are, in fact, identical, echo 
plus background has the same distribution as an echo from a stronger target. 

The signal arriving back from a single point target has a specific strength, a 
probability distribution of zero width. When it enters the circuits, and combines 
with the background signal, however, the combination no longer has a specific 
value, but a probability distribution whose width or standard deviation is a 
function of the relative strength of signal and background. The instantaneous 
output signal (the signal leaving the receiver) then can never be predetermined 
specifically, unless the target is a single point (or some equivalent nonrandom 
array) and infinitely stronger than the random background. 

Whatever the probability distribution of the output signal, that for an 
average of k such independent signals will have a standard deviation propor- 
tional to k-“'/, If the random component of the output signal is due to circuit 
fluctuations, those fluctuations can be trusted to achieve independence in the 
interval between consecutive radar pulses. When the major part of it is due to 
the random array of the target elements, say raindrops, the time to independence 
may be great compared with the interval between radar pulses. 

The interplay of the ratio of echo intensity to background noise intensity, 
m?”, and of k, the number of independent signals available, is displayed in Fig. 20 
for both fluctuating and steady echoes. Here lines of constant standard deviation 
g, taken as a fraction of the mean echo intensity, are drawn. Solid lines apply 
when the incoming echo fluctuates about a mean intensity A? before being mixed 
with the background noise in the receiver. The broken lines apply when the 
incoming echo has a fixed intensity S?; here, the fluctuations, of which oa is the 
standard deviation, are due entirely to the noise background. 

In both cases, increasing k consistently decreases the standard deviation 
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or 10" 


Fic. 20. Lines of constant standard deviation (expressed as fractions of mean echo inten- 
sity); kis the number of independent data available, m?is the intensity ratio (echo/background). 
Solid lines for fluctuating, broken lines for steady, echoes. 


(a k-“/*), In the case of a steady signal (broken lines), increasing the echo/ 
background ratio m? consistently decreases the standard deviation, although at a 
gradually diminishing rate. In the case of a fluctuating input signal (solid lines), 
increasing m? decreases o only until m? reaches a value of about ten. Higher 
values of m? do not appear to be any more useful; ¢ cannot, in fact, be reduced 
below A?/V/k. This brings out again the main thesis of the present paper: 
signals from rain and similar random arrays have a measure of uncertainty 
associated with them which can be overcome by averaging independent data, 
but cannot be improved (beyond m? = 10) by increasing the strength of the 
signals. Every such signal carries an uncertainty equivalent to m? = 1 before it 
ever enters the circuits. 

The solid lines of Fig. 20 are based on the theory of Sections 2 and 3 and on 
the results shown in Table I. Thus 


mr aD ee 
‘ai A! VkA?) Vk A? ot 


‘ 


where m? = A?/B? is the ‘“echo/noise ratio’, B? being the mean intensity of the 
background noise; o; is the standard deviation of the fluctuating output signal 
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when one independent datum is available (i.e. for k = 1). This quantity is equal 
(Table I) to the mean intensity of echo plus background. 






For the broken lines, the theory given by Goldstein (3) has been used. Gold- 
stein (his equation (175) ) shows that 
o:/(S° + B’) =v/1 + 2m*/(1 + m’*) 
where o; is the standard deviation of the steady echo mixed with the background 
when k = 1, while S? is the intensity of the steady signal. Using m? = S?/B? 
and ¢ = o1/~/k, the equation 


[24] o/S*= V1 + 2m” /m/k 


for the broken lines is obtained. 













Detectability of Weak Echoes 

Increasing k has just been shown to aid materially in the precise measurement 
of a fluctuating echo in the presence of a background noise; in addition, large 
values of k make possible the detection of fluctuating echoes which are weak 
compared to the background noise. 

In Section 3, the probability distributions P(J,) of the average of k indepen- 
dent signal measurements (equation [13]) and Q,(A?) of the true mean intensity 
in terms of J, (equation [14b]) were derived. Table I lists the corresponding 
standard deviations: 


[25] for P(J,): A’ Vk 
[26] for Q,(A*): Jpk/[(k — 2) Vk — 3}. 


The distributions themselves are sketched in Figs. 3 and 4. Their appearance 
suggests—and this can be proved rigorously—that both functions tend towards 
a Gaussian probability distribution with standard deviations [25] or [26] as k 













increases. 
Since the noise background too follows these distributions, it is evident that 
provided a sufficient number of observations of it are made, an exact value of the 






noise level may be assumed to be known. Let this value be a. 

Consider now that k independent signal measurements, each signal consisting 
of background noise plus possibly an echo, have been made. For simplicity, let 
k be large enough (say greater than about 50) for the Gaussian approximations 
to be usable. If the resulting average of the k intensities, J;, happens to coincide 
with the noise intensity a, the probability distribution of the true signal average 
A? will be a Gaussian centered at a and having the standard deviation of 
equation [26]. It then follows from symmetry that there is a half-chance that the 
true value of the echo + noise should be above a, i.e. that there is an echo. There 
is also a half-chance, of course, that the signal is due only to noise. 

Assume now that a fluctuating echo is being received, for which the true mean 
(echo + noise) power is c (>a). The probability distribution of the received 
signal for any value of k& is then equation [13] or its Gaussian approximation, 
and the probability R of this echo plus noise giving rise to a signal greater-than a 
may be worked out. This R may be interpreted as the chance of getting a return 
greater than a, and hence the chance of getting a return which permits the ob- 
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server to assert with greater than 50% confidence that he is looking at a real 
target. The 50% curves of Fig. 21 show this probability R as a function of 
c — 1 = m’, the echo/background ratio, for three values of k. 


10 ° 2 4 
echo / noise (m®) 
Fic. 21. The probability R that for a given mean echo, the average of k = 25, 100, or 400 
echo + noise intensities will inspire in the observer a confidence greater than 50, 80, or 95% 
that an echo exists. 


For a higher than 50% confidence, it is necessary to analyze the distribution 
Q, (A?) further. It is found, for instance, that, fork = 100, signals (echoes + noise) 
falling at 6 = 1.209a inspire an 80% confidence that an echo exists. The prob- 


ability R is then simply the chance that a signal of mean power ¢ gives rise to an 
average power greater than 1.209a in k = 100 cases. In this way, the 80% and 
95% confidence curves of Fig. 21 may be calculated. These curves are plotted 
against m*; the point 1.0 on the abscissa thus indicates an echo of the same 
power as the noise. This is the echo strength which for low values of & is often 
regarded as ‘‘minimum detectable’’. It has been noted before that in the case of 
random targets severa: independent data must always be examined together. 
For this reason it would appear that for fluctuating signals from random targets 
the threshold of detectability is in all cases well below noise level. 

The results shown in Fig. 21 apply when k independent signals are averaged. 
In Section 3 the alternative and practically more attractive method of analyzing 
k data, by counting the fraction p/k of them which exceed a given threshold, was 
described. The criterion of detectability of echoes utilized in this way may be 
constructed in an analogous manner. The resulting curves for R versus m? are 
very similar to those shown in Fig. 21 for averaging, except that they are 
displaced towards the right by a factor about 1.3, when k = 100; this means that 
in counting, the echo has to be about 1 db. stronger than in averaging, in order 
to be detectable with the same precision. 


SUMMARY 
Randomly distributed scatterers give rise to echoes widely fluctuating in 
intensity about a mean which is generally the only quantity of physical interest. 
Straight averaging of the data (whether of intensity, amplitude, or intensity 
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level is immaterial) is the most obvious way of obtaining the mean. The width of 
the fluctuations reduces by Vk, where k is the number of independent data 
involved. But in practice averaging is both difficult and costly. A more readily 
realized alternative (involving only on-off type of response in the equipment) 
consists of classifying the signals and counting the fraction falling above a given 
threshold. This procedure yields the desired mean slightly less rapidly than ordin- 
ary averaging, provided the threshold is picked judiciously (to within +4 db.), 
or else provided that a sequence of thresholds, spaced not further apart than 
8 db., is available. 

Another simple alternative to averaging consists of superimposing many 
independent traces on an A-scope type display. The details of the individual 
traces are lost, and a continuous gradation of luminance results. The location of 
the brightest portion of the screen corresponds to the peak of the distribution 
curve of the signal function displayed. The mean value of the signal intensity 
may then be deduced easily. This technique fails when the receiver is square-law, 
as the distribution law of signal intensity has its maximum at zero intensity. 
Linear or logarithmic receivers would seem almost equally suitable for the pur- 
pose. 

The above techniques involve independent data. From any emitted pulse, these 
are received at the rate of one per pulse duration. Shortening the pulse dura- 
tion, or alternatively, receiving echoes from each pulse by separate antennas, 
increases the amount of independent information available, though only at a 
cost in signal power. 

The number of independent data may be increased substantially (by a factor 
between 5 and 20) and at no loss of power, by transmitting consecutive pulses 
at slightly different radio frequencies. A difference in frequency of only about 
1/7, where 7 is the unvarying pulse duration, ensures complete independence of 
consecutive pulses. The frequency need only be varied until the natural re- 
shuffling among the particles produces the independence needed. 

In general, receivers with logarithmic response characteristics are most 
suitable with random targets, such as are involved in weather radar. Their 
ability to respond over the wide range of power encountered in such echoes 
reduces appreciably the incidence of saturation. Further, because they are made 
up of several successively saturating stages, logarithmic amplifiers permit 
readily the establishment of the several thresholds which are desirable in the 
classifying and counting method of signal analysis. 

The precision of measuring the intensity of echoes from steady targets observed 
against a fluctuating noise background improves by an increase in either k, 
the number of independent data available, or the echo/noise ratio m*. For 
fluctuating echoes, the echo/noise ratio is not quite so important; its increase 
beyond about 10 adds practically nothing to the precision of measurement, 
while an indefinite improvement in precision may be secured by a sufficient 
increase in k. Large values of this quantity similarly permit the detection of 
fluctuating echoes amounting in strength to only a small fraction of the back- 


ground noise. 
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INTERPRETATION OF THE FLUCTUATING ECHO FROM 
RANDOMLY DISTRIBUTED SCATTERERS. IT! 


By P. R. WALLACE? 


ABSTRACT 


The determination of the probability distribution for the magnitude of a radar 
echo from randomly distributed scatterers is a standard problem of the ‘‘random 
walk"’ type. We are concerned with the information to be gained about the 
scatterers from observation of the echoes, and the methods whereby this informa- 
tion may be extended. The method of averaging independent signals is investi- 
gated, and compared with the method of continuous averaging, which is found 
to be somewhat better. By appropriate variation of frequency of the radiation, 
it is found possible to obtain effectively independent determinations from a given 
region, and hence to derive more precise information about it. Finally, the con- 
tinuous variation of signal with time is studied, and a simple picture of the 
development of the signal proposed. The Markoff method is used throughout. 


1. INTRODUCTION 

If a radar signal is directed toward atmospheric precipitation and the echo 
from that precipitation studied, limited information about the scatterers may be 
obtained. Neither the size of the scattering droplets nor their density may be 
determined directly; but if a drop-size distribution is assumed, information about 
each may be inferred. Since the echo to be expected from a given distribution of 
scatterers is subject to wide statistical fluctuation, no precise information may 
be obtained from an individual echo. 

We note that at any particular time ¢ after emission of the signal, the echo 
being received is a superposition of waves back-scattered from all particles in the 
range of distance between }c(¢t — r+) and 43ct from the station, and within the 
angular range of the outgoing beam, 7 being the duration of the emitted pulse. 
Thus, with the passage of time, we scan a region of space whose near and distant 
frontiers are both receding from us at a rate equal to half the velocity of light. 
Echoes received at intervals of time 7 come from completely different regions 
and the signals received are therefore independent. If it were known that the 
density and size distribution of scatterers were constant over a number of these 
independent regions of length }cr, we might make use of the several independent 
meastrements of signal to gain more precise information; however, this presum- 
ably assumes part of what one would like to know. 

We have undertaken, in Sections 3 and 4, a calculation of the probability 
distribution of the precipitation parameters which can be inferred from a sequence 
of independent measurements, and also from a continuous recording of the echo 
from the same region. The latter is found to be somewhat, though not spectacu- 
larly, superior. Some further insight into this problem may be acquired by study- 
ing the development of the signal as the frontiers of the scattering region recede.* 

1 Manuscript received June 4, 1953. 
2 Contribution from the Mathematical Physics group of the Department of Mathematics, 
McGill University, Montreal, P.Q. 

*The analysis will also be applicable to the variation of signal as the region changes owing to 

lateral sweeps. 
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It has been found possible to make a rather simple picture on the basis of which 
the autocorrelation of signals separated by a given time interval may be studied. 

Finally, a proposal by Marshall (1) permits the determination of independent 
signals from a given region of view. This is accomplished by changing the fre- 
quency of the radiation used. This effectively reshuffles the phases at which the 
radiation is scattered from the various droplets—and does it in a random way, 
provided the phase change is appropriately chosen. We have developed in 
Section 6 the quantitative theory of this method. 

It should be remarked that we are neglecting throughout the random spatial 
reshuffling of the scatterers. This will be justifiable if measurements are made 
within a time appreciably smaller than that required for random radial displace- 
ments of the order of one wave length to take place. Ultimately, of course, this 
effect will itself so redistribute the phases as to create an independent signal. 

2. ECHO FROM RANDOM SCATTERERS BY THE MARKOFF METHOD 


Suppose that a pulse is emitted at wave length \ cm. and constant amplitude, 


and is h cm. long, so that it contains 
[2.1] ny = h/d 


wave lengths. 

Now consider a scatterer at linear distance r from the station. Its contribution 
to the echo will commence at a time 27/c from the beginning of the emission of 
the pulse; this signal will be 


iw(t—2r/c) ¢ ‘ / ‘ 
2.2 z= ae"! for Qr/e <t < (2r/c) +7, 
= 0 fort < 2r/cort > (2r/c) + 7. 
The amplitude a will depend on the size of the scatterer and also on its distance 
from the receiver—the latter owing to solid angle effects. If this complication is 
neglected—a justifiable approximation if relative distances from the source do 
not vary too much—and if the size distribution is assumed everywhere the same, 
the signal is given by 
[2.3] Ma ga ents ene 
n 


where @ is the average scattering amplitude, and K = w/c. The summation is 


over all particles lying in the range of distance 
[2.4] get — th < 1%, < 3et. 


To find the probability of a given value of Z = X + 7Y we must weight the 
probability of this value when the number of scatterers contributing is NV, with 
the probability that these are in fact N scatterers. If the linear density of scat- 
terers is o, this latter probability is given by the Poisson distribution 

ire 
A —N 
py = ie 
where N = }oh is the average value of N. 

To calculate the probability distribution of the signal from N random scat- 

terers, as in all later problems, it is most convenient to use the Markoff method. 
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We shall recall briefly the conventional treatment of the problem by this method. 
We consider an N-dimensional space with coordinates 7,(m = 1, 2, ..., N). 
Each point in this space for which all coordinates are between 0 and h/2 will 
represent a possible disposition of the N particles, all having equal a priori 
probability. Thus the probability that the real and imaginary parts of Z lie in 
the ranges (X, X + dX), (Y, Y + dY) respectively will be the proportion of the 
volume (3h)* lying between the surfaces 


cos 2Kr,= X, X + dX 
and also sin 2Kr,= Y, Y+ dY. 


If a function 6 which is unity in this region, zero outside, is defined, it is clear 
that the required probability is 

ahh ath 9\NV 
[2.6] Py(X, Y) dXdY = J Gy 6 és J arg(2) 6. 

0 0 


Now let us introduce the two-dimensional Fourier transform @¢y(é, 7) of 
Py(X, Y). Taking transforms of both sides of [2.6] we find that 


on = [Jo(p)]” 


where p? = & + 7”, one Jy) coming from the integration over each coordinate 7,. 
Here J designates the conventional Bessel function of order zero. Therefore 
the probability of a signal in (X + dX, Y + dY) is given by 


P(X, Y) = a by Py(X, Y) 


5 , j ‘op [— N(1 — Jo(p))] e “**™ dé dn. 


~ (2m)? 


In the limit of large N it is sufficient to replace 1 — Jo(p) by its approximate 
value jp” for small p; the integral is then easily evaluated to give 


9 
[2.8] P(A) dA = 74 


5 


as 


—A? 42 
fer ae 


where A? = Na®; or alternatively 
[2.9] P,(I) dI = oo 
I 

where J = A®. This is the familiar ‘random walk” solution (2). We note that 
what we have obtained is a distribution in intensity A?. 

Before going on to specific applications, we might mention two advantages 
of the Markoff method in problems of this type. 

1. It yields first the Fourier transform of the distribution required. Thus, even 
if the inversion cannot be carried out, we may obtain the moments of the dis- 
tribution by differentiation; e.g. if g(£) is the transform of the distribution p(x), 


=E= —iq' (0), 
[2.10] (Ax)? = [g’(0)]* — q’’ (0), ete. 


2. In problems concerned with chains of processes, we can make use of the fact 





998 CANADIAN JOURNAL OF PHYSICS. VOL. 31 


that the probability distribution of the result of the chain is the convolution of 
the probability distributions of the individual steps. Therefore, using a familiar 
theorem of Fourier analysis, the transform of the complete probability is the 
product of the transforms of the component probabilities. 


3. DISTRIBUTION OF THE AVERAGE OF INDEPENDENT SIGNALS 

Suppose that, either by lateral sweeps or by measurements at pulse-length 
intervals, we obtain & different independent intensity measurements J;,... , Z;. 
Let us designate their average as J,. Now by the convolution theorem of the 
previous section, the transform of the probability distribution for the sum of the 
determinations is the product of the transforms of the probabilities of the in- 
dividual determinations; since the latter are identical, 
[3.1] m(€) = [ge(é)]* 
where go(€) = transform of P»(J), 


aw(€) = transform of the distribution of RJ, = Y,. 


Since direct calculation gives 


‘ . sft 
[3.2] go(&) = E+a/f’ 


we get, on evaluating the inverse transform of [qo(é)]* by contour integration, 
for the distribution of the sum Y,, 


rk—1 
| k Ye/I 


a0 =k é 
pel (f)” (k— 1)! 
Consequently for the distribution of J, 
k k-1 

ie Rk’ J; —kJk/T 
= —s 

(1)"(k — 1)! 

This result does not quite answer the practical question, namely: to what 
accuracy can we determine the mean intensity 7, which depends on the mean 


(3.4] PLS) 


scattering amplitude from a single particle and the particle density, when our 
information consists of k independent intensity measurements? However, the 
distribution of I for a given J; can be deduced from [3.4]. The relative prob- 
abilities are 
eed l kKJk/ 
[3.5] (Ty e 
and the absolute probability may be obtained by normalization. The normaliza- 
tion integral is only convergent for k at least 2. This arises from the fact that, 
for a single intensity measurement, there is always an overwhelming probability 
that it came from an indefinitely large J! The fact that we do not really admit 
this overwhelming probability is based on considerations quite extraneous to the 
measurement made, that is to say, our conviction on physical grounds of the 
finite character of the disturbance from which the echo comes. 
For k > 2 we find the distribution 

ay I —kJIx/I 

e ; 


[3.6] Q.(1) = (e-2)! F 
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From these we calculate immediately 










ke 2 
(= 2)" = 8) 7 
In general, the 2™ moment is finite only for k > n + 2. In particular, at least 
four independent measurements are needed to give information about J with 
finite probable error, and even then AI = 2J,,a very wide spread. For k = 10, 
AI = 0.47J,o, a result which is considerably greater than the r.m.s. deviation in 
J for fixed J, namely AJ, = 0.31/. In the limit of large 2, the fractional r.m.s. 
errors tend to become equal. 

The distributions P,(J,) and Q,(/) are shown in the accompanying paper of 
Marshall and Hitschfeld (1), Figs. 3 and 4. 


[3.7] 





(Ar)’ = 

















4. DISTRIBUTION OF CONTINUOUSLY AVERAGED SIGNALS 







Suppose now we consider the average signal over m pulse lengths, and see 
whether it offers the possibility of greater accuracy in the determination of J. 
It is clear that it must be at least as good as the averaging of independents, 
since all data contained in them are used in the continuous averaging process. 
Suppose we take the time at which the averaging is begun as t = 0. The 
particles then contributing to the signal lie, according to [2.4], in the coordinate 







range 
(i) —h/2<%m% <0 (Region R;). 







Particles in this region contribute to the average for time < r. 






In the region 


(ii) O< 47, < (m — 1)h/2 (Region R2) 







are particles which do not contribute to the signal either at the beginning or the 
end, but each contributes for a total time +. 







In the region 
(iii) (m — 1)h/2 <7, < mh/2 (Region R3) 










lie particles which contribute to the signal at the end of the interval of averaging, 
but contribute altogether for time < r. 









We are now interested in the intensity 7 = ZZ* (see equation [2.3]). The time 
average of J is immediately calculable and is 


[4.1] ie = Ey\(r,) + 2 Z. cos 2K (tq — rn) E2(Tn; Tn’) 


Tn¥FTn’ 







where 


Ein) = - (1 + 2r,/h) in Ri, 






[4.2] > joie 


m 


9 = a 
«3 ea i ae 






C 
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and 
MiN (fn, Tn’ 
Cc 


E2(tn, Tn’) ) a5 Tr’ in Ri), 


2rn — Tn! 
ae (at least one of rp, 7, in Re and 


Irn ras Tn’ | 2 zh), 


1 _ g MAX (Fn) Mn’) — 3 (m ae ee 


c 


m 


It is evident that it would be difficult to obtain the distribution of J,; however, 
it is quite feasible to get its mean-square deviation 


Al, otis (Te a, (I,)°}}. 


In the calculatiom of 72, the cross-product of the first and second terms in 
[4.1] drops out, as do cross-products between the terms of the second series. Also, 
in the squared terms arising from the second series, cos? 2K (7, — r,) may be 
replaced by 3[1 + cos 4K(r, — 7’)]; the contribution from the oscillating term 
is then 1/m» of that from the constant (”) being the number of waves in the 
pulse), and may therefore be neglected. 

Direct computation of the integrals then yields the result: 


[4.4] (AI,)? = [?(4m — 1)/6m? 


where J is the mean intensity. 

Now the region over which the average has been calculated contains (m + 1) 
independent regions. The r.m.s. deviation of the average of the (m + 1) indepen- 
dent signals is 


[4.5] Ang: = I/\fm + 1. 


The ratio of the mean square deviations, which we shall call the ‘factor of 
improvement”’ f, is 
9 
(Al,)" _ (4m — 1)(m +1) 


6 a j= 5 
- | CAT na) 6m" > 


We list a few values in Table I: We note in particular that there is no improve- 
ment unless the average is taken over more than one pulse duration. At best, 
the reduction of the r.m.s. deviation due to averaging does not exceed 22.5%. 


TABLE I 


2 3 4 5 10 20 50 © 


0.88 0.81 0.78 0.76 0.72 0.69 0.68 2/3 


Put differently, we can give the number of independent signals which give the 
same r.m.s. deviation as the continuous average over the duration of m pulses. 
This is given in Table II. 
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TABLE II 





3 










4.9 6.4 


It should be noted that if we are again interested in obtaining the distribution 
of true mean intensity J from a known signal, a further analysis is necessary. 
Roughly, we may estimate the result in this way: we assume that the continuous 
average is equivalent to the average of (m + 1)/f independent signals, as indi- 
cated in Table IJ. If the average is over 10 or more intervals, to a good approxi- 


mation 


ee ae) er ee 
~ 3(m—3)*(m—1)° 5 3m~‘ 


[4.7] (Al) 
(see [3.7]); so that we get the same factor of improvement as in the converse 


problem. This is to be expected, in view of the similarity of the distributions for 
averages over reasonably large numbers of pulse lengths. 
5. CONTINUOUS VARIATION OF THE SIGNAL INTENSITY* 

We now turn to a somewhat different problem: that of the variation of the 
signal as the region scanned is varied. We ignore the motion of the individual 
scatterers, but consider the evolution of the signal as the region from which the 
signal is received recedes continuously into the distance. Alternatively, the 
analysis applies equally well to the case in which the direction of the outgoing 
beam is altered, and the signal varies owing to a lateral sweeping. Finally, 
although the continuous variation is studied, the results may be used to relate 
signals received from any two regions which partially overlap. 

Let us look for the probability that the complex vector representing the echo 
be Ay at an initial time, and A at a time ¢ later. The signal at the initial time will 
come from a region of length 3/4, which we shall designate 0 < 7 < $h; at the 
later time it will come from the interval ct/2 <r < ct/2 + h/2. Let us assume 
that there are N, particles in (0, ct/2), contributing to the first signal only; N2 
in (ct/2, h/2) contributing to both; and N;3 in (h/2, h/2 + ct/2) contributing to 
the latter one only. The first signal will then have components Xo, Yo given by 


Xo = > cos 2Kr, + >> cos 2Krp, 
Ni Na 


Yo = > sin 2Kr, + )> sin 2Kry. 
Na 


Ni 
The second is given by 


A= z= cos 2K (r, — 3ct) + = cos 2K (r, — $2), 


Ns 


Il 


Na 

Y = > sin 2K(r, — 3ct) + D> sin 2K(r, — ct). 
Na Na 

*Results similar to some of those in this section have been given by Siegert (3) (1943) and 
(4) (1946) in connection with the variation of signal solely due to the motion of the scatterers. We 
include a discussion of this problem in order to point out a simple result (formula |5.3]) and its 
interpretation (Fig. 1), which make possible a graphic picture of the development with time of the 
sort of signal with which this paper ts concerned. 
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Now the probability of a signal intensity represented by a vector in dXq, dYo 
initially and in dX, d¥Y at time ¢ is given by 


a ( 9 )" ey whet Ni nth Na ahhthet N> 
4 eee = dx, dx J dx, . 650 
ct h—ct ct J, I] ‘Viecpe Oh IT ’ 


where the symbol [| signifies a product, and 
59 = 1 if Xo, Y, lies in the required range dXo, d Yo 
= (0 otherwise, and 
= 1 if X, ¥ lies in the required range dX, dY, 
= () otherwise. 

Representing the 6’s in the usual way by Fourier transformations, we can 
determine the probability distribution for given Ni, No, N3. Weighting with the 
usual Poisson distributions for these quantities, we find that the probability 
that the original signal is in (Xe, Xo + dXo; Yo, Yo + dYo) and the second one 
in (X, X + dX; Y, Y+ dVY) is 


4 . ® 
dX od Yod Xd “(3 ) J growl ge de, J ie 
oT 
[5.1] X exp {—o[3ct(1 — Fi) + 3(h — ct)(1 — Fe) + 3ct(1 — Fs)} 
where 


exp {7(& cos 2Kr + mo sin 2Kr)} dr, 


‘ whh 
i ; | exp {i[£) cos 2Kr + no sin 2Kr + ~cos 2K(r — 3ct) 
ie 


het 


+ nsin 2K(r — 3ct)]} dr, 


ehh+}ct 
3 exp {i[f cos 2K(r — 3ct) + n sin 2K(r — }ct)]} dr. 
ry @Pth 
Making use of the fact that oct/2 is very large to expand the F’s to quadratic 
terms as in [2.7], [5.1] comes out in terms of Gaussian error integrals. Evaluating, 
we find that the probability in question is 


P(Xo, Yo; X, Y) dXodYo dX dY 


‘ 


[5.3] 1 i lA Ais aa 2 
= dX) dY)dX av (4ee7)( pagcig ees Le ) 
wl 


where A,? = Xo? + Fe’, 

A = 1 — t/r where 7 is the pulse duration, 

Ay’ is the vector obtained by rotating the vector Ay with components (Xo, Yo) 
through the angle —wt, w being the circular frequency of the radiation, 
and A is the vector with components (X, Y). 

The first bracket in [5.3] represents the probability distribution of the original 
signal; the second, which may be designated P(A» — A), is the probability dis- 
tribution after time ¢, when the initial signal is given by Apo. 

If we ignore variations with the frequency of the emitted radiation—consider, 
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that is, the signal only at intervals of one or several r-f. periods—the distribution 
at time ¢ is a Gaussian (normal error function) centered on the point AAo (a 
point moving uniformly in toward zero from the initial value Ay) and having a 
standard deviation J(1 — A?), which increases from zero at t = 0 to J, the 
value for an independent signal, at ¢ = +r. This is illustrated in Fig. 1. 


Yr =| 


Fic. 1. Straight line represents initial signal amplitude and phase. At later instances, as 
t approaches 7, the phasors are subject to two-dimensional Gaussian probability distributions. 
The surfaces representing these distributions are centered, at the times indicated, at the 
emphasized points, and have standard deviations equal to the radii of the circles. When 
t = 7, the distribution is that of an independent phasor. 


The result [5.3] makes possible the calculation of the frequency distribution 
of the echo, as calculated in Siegert’s papers. Knowledge of the frequency distri- 
bution, however, conveys no useful information about the precipitation which 
returns the echo. The preceding treatment of the development of the signal is the 
basis for autocorrelation studies, and in addition, provides a simple and graphic 
picture which aids in the understanding of the correlation of successive signals. 


6. THE METHOD OF FREQUENCY VARIATION 


Let us consider now the following problem: suppose that a signal of intensity J» 
is received from a given region when the circular frequency is w9; what then will 
be the probability distribution for the intensity J which will be received from the 


same region when the frequency is w? We shall assume, as usual, that there is not 
sufficient time between these measurements to allow the scattering particles to 
undergo any appreciable movement. If, however, there should be an appreciable 
random or turbulent motion, it would clearly have the effect of making the 
successive measurements more independent of each other. 

Let (Xo, Yo) be the real and imaginary parts of the complex vector 


[6.1] Zo = » gn (x. = = 


Tn 
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where the sum is over the N particles in a half-pulse-length range of length 5h. 
Suppose also that (X, Y) are the corresponding parts of the vector 
[6.2] Z=>ye°*" (x ~ “) 
corresponding to an altered frequency w. The probability that (Xo, Yo) lies in a 
given range (dX, dYo) and (X, Y) lies in (dX, dY) is 

\Y ei 
[6.3] ) 606 dx, oie dxy 

h 0 


where 69, 6 are one when (Xo, Yo) and (X, Y) respectively are in the ranges in 
question, zero otherwise. 
Representing 69, 6 as Fourier integrals, transforming to polar coordinates: 
(Xo, Yo) = (Ag, A), 

(X, Y) > (A, 4), 
and for the transformed variables: 

(fo, no) — (po, Xo), 

(€,) — (p, x) 


we obtain, on integration over 69 and @, 


1 . . : 
(Qa)? Aj.dA oAdA J Po dpo dxo | pdpd x . J (po Ao) Jy (pA) b* (po, p) 


(6.4] 
where 
9 vhh 
[6.5]  b(po, p) = “I exp {i[p0 cos (2Kor — x0) + pcos (2Kr — x)]} dr. 
0 
Weighting with the Poisson distribution for the number N of contributing par- 


. ° N ° “eer. 
ticles, and summing, b° (po, p) becomes in this integral 
g-Ne(1—D(p0.0)) 


where No = mean of N = oh/2. 
The next step is to expand 6 to quadratic terms in po, p. If this is done, and the 
integrations over x, xo are performed, we obtain the distribution 


[6.6] A,dAo . AdA j Jo(A po) podpo J so(Ao)ode . exp { = 1No(p 1. po )} 
x Jo(3tNogoppo) 
where Jo is the Bessel function, and 
. _ sin 2hAK 
[6.4 »™ shAK 


Finally, the po, p integrals may be evaluated. Putting J) = Ao’, J = A?, we 
find that the probability of intensity J» at frequency wo and J at frequency w is 





J. e To gIq dl —— exp \- I \ 


[6.8] Mo No(1 — go 


” 7{ Berlin) | 


N,(1 —= 26 ) 





WALLACE: INTERPRETATION OF FLUCTUATING ECHO. II 1005 


The factor in the square bracket, P(J) dJ, is the probability that the intensity 
at frequency w is in J, J + dI, that at frequency w9 being Jo. In Fig. 2, P(J) is 
plotted for several values of hAK/2z. In Fig. 3, go? is plotted as a function of AK 
and Av (solid line). 


x27 Ay! = leldwl _ hak! 
@ eT 


Fic. 2. Probability distributions of the signal intensity J at transmitter frequencies 
vo + Av, given that the intensity is Jp at frequency vo. Each curve is drawn for a fixed value of 
x = 7| Av|, r being the pulse duration, # the pulse length, w = 2rv, K = w/c. mg is the initial 
number of wave lengths per pulse length. The curve for x = 1 is the exponential distribution of 
equation [2.9], indicating that J has become completely independent of Jo. 


The second signal is now completely independent of the first provided gy = 0, 


i.e. 


(n = integer) 


or 


[6.9] 


no being again the number of waves in the outgoing signal. Thus, by virtue of a 
succession of small variations of frequency, independent signals can be obtained, 
and their average used to obtain an estimate of the mean scattering intensity 
for a given half-pulse length. This frequency variation should, of course, remain 
sufficiently small that the scattering cross sections of the particles are not appre- 
ciably affected by it. 

It may be noted that, for Aw/w > 1/mo, go? is never larger than 0.05, so that 
there is a substantial degree of independence of consecutive measurements at any 
such frequency shift. Thus, provided one can afford to make the frequency shifts 
sufficiently large, no very precise control of frequency is necessary. 

If we take account of the fact that the emitted wave will not be precisely a 
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pulse shapes 


-DIOKL = ria vf 
en 


Fic. 3. Variation of p, the correlation coefficient between signal intensities J (at frequency 
vo + Av) and Jy (at frequency vo), with Av. Solid line, for an ideal rectangular pulse; broken 
line, for a trapezoidal pulse. 

It can be proved that p, defined as I & Ip — (1)? / I? — (1)?, equals the factor go? of equa- 


2 


tion [6.7] and y? of equation [6.17]. 


square wave, the result will be somewhat altered. Suppose that at any time ¢ the 
outgoing wave amplitude is modulated by a factor f(¢/r). We may define the 
factor so that its mean value is unity, so that 


al 
| f(u) du = 1. 


t we shall assume to be measured from the initial moment of propagation. We 
shall write F(u) = f(1 — w), so that 


al 
[6.10] J F(u) du = 1 
0 


also. The results of the preceding calculation are now altered in that b(po, p) 
(equation [6.5]) becomes 


9 ahh ) 9 
[6.11] 6'(p0, p) = | exp , ir(™) cos (2Kor — xo) 


0 


+ pcos (2Kr — » | dr. 


If, as will in general be the case, F?(2r/h) = W(2r/h) varies only slightly over a 
wave length of radiation, 
[6.12] 1 — b'(po, p) = 4[a(po? + p?) + 2pp.G] 


where 


[6.13] a= J v(u) du =1+ J (ary du, 
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AF being (F — 1), the deviation of the form factor F from its mean value. Also 
[6.14] G = C(hAK) cos Ax + S(hAK) sin Ax 






where C is the Fourier cosine transform of ¥(u) for argument AK, and S is the 
corresponding sine transform. 
Putting 


[6.15] VC? +S = H(hAK), 








we see that a = /7(0). We may now calculate the distribution exactly as before; 
the probability that Jo lies in dJ,¢ and J in dI is 






be Bs Fe. . \ 
Pa- 7)" " h- ” - se 


[6.16] 
x (ptt a i) dI dl, 


a (I, Io) dI dIo = 







where J is the mean intensity, and 





oe __H(hAK) 
16.17} 1=~~HO) ’ 


H being the magnitude of the Fourier exponential transform of the modulation 






function y. 

We can consider as the simplest case of interest, that in which the intensity 
falls off linearly over the duration of the pulse: 
[6.18] yu) =1+B(u-—3), (O<u <1). 













Introducing u(<7/2) such that 
tan wp = (1 — 38)y/8 







and 








we can easily verify that 





Cy) = - sm ay ® sin (4y — vw) + ® sin y, 









6.19 ; 
(6.19) S(y) = 0 & cos (3 y — u) — - cos: y. 











2 





We then see that C + 7S and therefore H can never be zero, and there is no 

frequency shift for which the measurements are completely independent. 
However, from a practical viewpoint, this result is not unduly perturbing 

since, for the values of AK which give y = 0 for a square wave, namely, 

hAK = 2x X integer, 









y* becomes 





8° ae 


(hAK) ~~ 47° X (integer)” 


Since it is unlikely that 6 will be >1 (ratio of intensities at two ends of the 
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wave 3:1), the correlation is quite small and a high degree of independence 
is assured. Once again, too, the correlation is consistently quite small for 
hAK > 2r. 

A comparison of the results for a rectangular pulse and for one of linearly 
decreasing (or increasing) intensity is given in Fig. 3. 

A second form of modulation which might be considered is one symmetrical 
about the center point of the pulse: 


[6.20} ¥(u) = ¥(1 — u). 


In this case it may be shown immediately that 
j 
A(hAK) = af ¥(u) cos hAK(u — $) du. 
0 


If, on the half-pulse range 


v(u)=1+8'G —w), 


we get 
, 


Tag ~ oo ma | 


1 
eee TAg\ oe neta ae 
H(hAK) = afc + 38’) nAK SO thAK , 
This is zero whenever 


hAK = 4nxn, (n = integer). 


In this case we get complete independence for intervals twice as great as those 
which give it for a rectangular pulse; however, halfway between there is again 
approximate independence in the sense of the previous paragraph. 


7. CONCLUSIONS 


Information obtainable about atmospheric disturbances from radar echoes is 
such that one cannot differentiate between “‘weak”’ scattering from an extremely 
large density of scatterers and “‘strong”’ scattering from a smaller density, so long 
as (number of scatterers) X (intensity per scatterer) is the same. If a drop-size 
distribution is assumed known, however, the number and size individually can 
be inferred from the mean intensity scattered from a given region. 

To establish the mean scattered intensity with reasonable accuracy requires a 
considerable number of individual intensity measurements. These may be made 
over a number of different regions, provided it is known (or assumed) that these 
regions are homogeneous. If the echo from these regions is investigated by 
recording the return from a given emitted signal, it is found to be somewhat 
better to average the echo intensity continuously rather than only at ‘such 
intervals that the signals are independent. 

On the other hand, the mean intensity from the single region which contributes 
to the echo at a given instant can be estimated by using independent signals 
from this one region. This is accomplished by varying the frequency of the 
emitted radiation. For a square wave, or a wave symmetrical about its center 
point, complete independence can be attained by appropriate variations of 
frequency. For an arbitrarily shaped wave, provided its intensity variation is 
not too-drastic or rapid, a high degree of independence can be obtained by 
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varying the frequency by any amount larger than the amount necessary to add 
another complete wave length or two to an emitted signal of constant total 
length. 


* The work reported in this paper was supported by the Geophysics Research 
Division, Air Force Cambridge Research Center, under contract No. AF—19 
(122)-217. 
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A THEORETICAL INVESTIGATION OF THE NUCLEAR RESO- 
NANCE ABSORPTION SPECTRUM OF AI” IN SPODUMENE! 


By G. LAMARCHE? AND G. M. VOLKOFF 


ABSTRACT 


A theoretical investigation is presented of the energy levels of a nucleus of spin 
I = 5/2, a given magnetic moment yu, and electric quadruple moment eQ, placed 
in a uniform magnetic field Ho, and a crystalline electrostatic potential ¢ without 
axial symmetry. The dependence of the energy levels, line frequencies, and tran- 
sition matrix elements on the dimensionless parameter R= 4H /eQ¢zz is cal- 
culated for Al?’ forone orientation of aspodumene crystal in the magnetic field Hp 
over the interesting range of the parameter R linking up the pure quadrupole 
with the Zeeman spectrum. Other orientations of the crystal are discussed briefly. 


In a recent paper (7) a summary was given of the various experimental and 
theoretical studies that have been made to date on the electric quadrupole 
interaction of non spherically symmetric nuclei (with J > 3) with the electric 
field gradients existing at the sites of such nuclei in single crystals. It was 
pointed out that, to date, experiments involving an externally applied static 
uniform magnetic field 17) have been performed in two limiting cases: R < 1 
and R> 1, where R = 4uH)/eQ¢-. is the ratio of the magnetic interaction 
energy »H, to the quadrupole coupling energy eQ¢z2/4. 

A search is being made in this laboratory for a suitable crystal with a field 
gradient tensor lacking axial symmetry (7 # 0) in which the nuclear resonance 
absorption spectrum might be observed over a sufficiently wide range of 
values of R to permit an experimental study of the gradual transition from a 
pure quadrupole spectrum (R = 0) to the Zeeman line split into several 
components by the quadrupole interaction (R > 1). Experimental data are 
available in Reference (5) on the quadrupole coupling constant of Al?’ in 
LiAl(SiO;)2 (spodumene) and on the properties of the electric field gradient 
tensor $;; at the Al*’ sites. The expected dependence on R of the frequencies 
and relative intensities of the nuclear resonance absorption lines of Al?’ in 
spodumene was calculated for one particular crystal orientation to determine 
whether this crystal could be used for the proposed experimental study of the 
transition from the pure quadrupole to the Zeeman spectrum. These calcu- 
lations are briefly summarized below. An attempt to perform the proposed 
experiment with spodumene is now in progress. 

To study the gross line structure in terms of the number and frequency 
separation of its components, disregarding the line width of the individual 
components, we may neglect dipole-dipole interactions, and consider each 
nucleus as being individually acted upon by the given Hy and the given crystal- 
line field gradient ¢;;. The Hamiltonian then consists of two parts: 


[1] 9 = Hot Hz. 


1 Manuscript received May 22, 1953. Contribution from the Department of Physics, Univer- 
sity of British Columbia, Vancouver, B.C. Abstract presented at the Royal. Society of Canada 


meeting, June, 1958. 
2 Holder of a Bursary from the ‘‘Ministére du Bien-Etre Social et de la Jeunesse’’ de la 


Province de Québec. 
1010 








LAMARCHE AND VOLKOFF: SPODUMENE ABSORPTION SPECTRUM 1011 


Here §g, the part due to the quadrupole coupling, is given (Refs. (2) and (3)) 


; = eQ¢.: 2 2 2 ‘9 
2 = — :— - 
where x, y, 2 are the principal axes of ¢;; and 7 is the asymmetry parameter 


defined as 
[3] 1 = (ozz — Oyy)/O2: 
The Zeeman energy §z is given by: 
[4] Dz = — gbl+Ho. 
The problem consists in solving the secular determinant of the operator 
[5] | Gum — Ednne | = 0 


for the energy levels, and in computing the frequencies and relative proba- 
bilities of transition between them. We first of all note that the problem is 
considerably simplified if 7) coincides with one of the principal axes of ¢;;. 

When the quadrupole energy matrix $g is written down in a representation 
in which the x, y, z axes, to which I,, I,, and I, are referred, coincide with the 
principal axes of ¢;;, as is done above, the only nonvanishing matrix elements, 
as shown in Ref. (3), are those for which | m — m’ | = 0 or 2. In a coordinate 
system that does not coincide with the principal axes, the |m — m’| = 1 
matrix elements also do not vanish (Refs. 1 and 6). 

The magnetic matrix § z has only diagonal elements if the axis of quanti- 
zation of spin coincides with the direction of the magnetic field H/o. Otherwise 
it has nonvanishing | m — m’ | = 1 matrix elements as well. 

Hence if, and only if, Hy coincides with one of the principal axes of ¢;; can 
the | m — m’ | = 1 matrix elements of both interactions be made to vanish 
simultaneously. With this choice of the H» direction, the only nonvanishing 
terms of the matrix § are those for which | m — m’| = 0 or 2. For half 
integral J the determinant is then readily factored into two determinants each 
of which leads to a characteristic equation of degree J + }. The eigenstates 
defined by each of these two subdeterminants are linear combinations of one 
of two subsets into which the functions Ym which diagonalize I, may be sub- 
divided. One set includes ¥» with m = 4 + 2n, and the other with m = — 3 
+ 2n, where n is an integer. 

In the most general case when /J, does not coincide with any of the principal 
axes of ¢;;, there will always be the |m — m’| = 1 nonvanishing matrix 
elements present, no matter what representation is used. In this case the 
determinant does not factor easily, the characteristic equation obtained from 
equation [5] is of degree 2 + 1, and each eigenstate may involve all the pm. 

To obtain information useful in planning experiments involving Al?’ in 
spodumene, we choose J = 5/2, 7 = 0.95,* and treat the case when Hy coin- 
cides with the z principal axis of ¢;;. The sixth-order secular determinant then 


*The numerical values used in the present paper were taken from the preliminary communication 
on Al?? in spodumene Ref. (4), and differ slightly from those of Ref. (4). 
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gives rise to two cubic secular equations whose coefficients involve the para- 
meters R and n. These two equations are: 

3R 2 or 2 9 $ ‘ 

30 + 400 (13R° + 24R + 7n° + 21)A 

+ qgoglt BR" + 8R?¥ ISR + BR + 4(1 — 7) = 0 | 
where \ = E/eQ¢,,. We have plotted in Fig. 1 \ as a function of Rfor0 < R 

< 4. The \7 and X¥ (i = 1, 2, 3) are respectively the solutions of equation [6] 
corresponding to all upper and all lower signs.* The corresponding eigenstates 






—-+ 





(6] 


















0.5 





0.4 0.8 1.2 1.6 20 2.4 2.8 .2 3.6 40 


Fic. 1. Energy values for J = 5/2, » = 0.95 in terms of A= E/eQ¢zz as a function of 
R= 4uHo/eQdzz, for 0 < R < 4. Ho is along the z principal axis of ¢;;. The arrows represent 
the possible absorption lines, and are numbered from 1 to 9, the first five being the Zeeman 
transitions for R>> 1. 


are of twoclasses as mentioned above: class L involving ¥» withm = — 3, 3, $, 
and class M with m = — §, —}, 3. 

Transitions produced by an r-f. magnetic field placed at right angles to Ho 
are only possible between levels of different classes (Am = + 1). Their fre- 
quencies are easily calculated by taking the difference between the eigenvalues 
d of Fig. 1 multiplied by the known quadrupole coupling constant of Al?’ in 
spodumene eQ¢.-/h = 2.960 + 0.010 Mc./sec. They are plotted (in Mc./sec.) 
in Fig. 2 as a function of R for 0 < R < 4 (or of Hy = ReQdz2/4u). 

i a values of R these solutions agree with those obtained by the perturbation theory of 
e].\4/). 
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Fic. 2. Transition frequencies v (in Mc./sec.) as a function of R (or of Ho = ReQdzz/4y). 
Lines have been labelled intense or weak in accordance with relative values of the squares of 
matrix elements being greater or smaller than unity on the arbitrary scale used in Fig. 3. 


Finally, the relative transition probabilities are determined largely by the 
squares of the absolute value of the matrix elements of the perturbing operator 
which is proportional to (I, + zI,). Since the elements of this operator are 
well known in the ¥ representation, and since the eigenstates of our problem 
can be expressed on solving the secular determinant as known linear combina- 
tions of ¥m, the required matrix elements may be easily computed. The squares 
of these matrix elements are plotted in Fig. 3 in arbitrary units as a function 
of R. In Fig. 3a we have the squares of the matrix elements corresponding to 
the Zeeman lines, i.e. transitions between levels which are adjacent for R > 1, 
while the others, the so-called ‘‘forbidden lines’, are given in Fig. 30. 

From these graphs we may note the following. At Hy = 0 the pure quadru- 
pole spectrum will consist of two strong lines of almost the same frequency 
(0.789 and 0.758 Mc./sec.) plus a much weaker line at the sum of these two 
frequencies (this line would be absent if 7 = 0). As Ho increases, the spectrum 
gradually changes as shown in Fig. 2 until eventually the Zeeman line, split 
into five components by the quadrupole interaction, is established. At inter- 
mediate fields (around 260 gauss) there is apparently a chance of finding 
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Fic. 3. Squares of the matrix elements of the perturbing operator (in arbitrary units) asa 
function of R = 4yHo/eQ¢zz. 


“‘extra lines’”’ as shown in Fig. 2 (vs for 200 < Hy < 400 gauss and »; for 100 
< Ho < 750 gauss). 

If Ho coincides with one of the other two principal axes of ¢;;, the solution 
of the problem is in principle the same, and differs only in numerical details. 
The axes can be relabelled, and y and ¢,, conveniently redefined. 

In the more general case where H, does not coincide with any of the principal 
axes of ¢;;, transitions between any pair of the six levels become possible in 
principle, so that 15 lines may be expected, although some may again be of 
too low intensity or too low a frequency to be observable. 


One of the authors (G.L.) wishes to acknowledge financial assistance from 
the Ministére du Bien-Etre Social et de la Jeunesse de la Province de Québec 
in the form of a bursary, and the other (G.M.V.) the continuing financial aid 
of the National Research Council of Canada in the form of grants-in-aid of 
the experimental program which gave rise to this work. 
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THE DELAYED CATHODE GLOW PRODUCED IN A NEON DISCHARGE 
TUBE AT LOW PRESSURE 


By DONALD J. ECKL 


When the light output of the absorption tube used in the measurements 
of the decay of the populations of neon atoms in the 2p°3s levels, described 
elsewhere in this journal,* was «xamined by observing on an oscilloscope the 
signal voltage produced by a photomultiplier tube, a satellite pulse was 
observed. This pulse occurred some time after the original light output pro- 
duced by the excitation voltage across the discharge tube had dropped to zero. 
When a baffle was placed so that light from the cathode region of the discharge 
tube could not reach the photomultiplier, the satellite pulse was no longer 
observed. 

Photographs showing the photomultiplier tube output for the excitation 
pulse and the accompanying satellite are shown in Fig. 1. A curve of the delay, 





Fic. 1. Light intensity output from absorption tube at various neon pressures obtained 
with a 931A photomultiplier tube. Duration of pulse is 500 usec. 
a. Pressure is 2.6 mm. Amplitude is one-half scale of other three figures. 
b. Pressure is 1.9 mm. 
c. Pressure is 1.0 mm. 
d. Pressure is 0.6 mm. 
measured as the time between the end of the excitation pulse and the maximum 
amplitude of the satellite, as a function of pressure is given in Fig. 2 and has 
the form 
k k = 125 mm-ussec. 
T=-+4+1, 
p L = 28 usec. 
where T is the delay in microseconds and is the pressure in millimeters. 


* Can. J. Phys. 31: 804. 1958. 
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Fic. 2. Cathode glow delay against pressure. 


A spectroscopic check of the light from the satellite shows the presence of 
spectral lines originating from the 2p°3>p levels. Since all electronic levels above 
the 2p°3s levels have extremely short lives, the source of such levels associated 
with the satellite pulse is probably positive ions. The delay then should be 
associated with the recombination time for such ions in the neighborhood of 
the cathode. The light in the satellite pulse would represent a source of light 
from atoms in an essentially field-free region. 

RECEIVED APRIL 24, 1953. 

McLENNAN LABORATORY, 

UNIVERSITY OF TORONTO, . 
TORONTO, ONT. 





INFLUENCE OF STRIATIONS ON THE PLASTIC DEFORMATION OF 
SINGLE CRYSTALS OF TIN 


By K. A. JACKSON AND B. CHALMERS 


In recent years, much work has been done on substructures formed during 
plastic deformation. The purpose of this note is to set forth some observations 
on the effect of an existing substructure on a substructure subsequently formed 
by deformation. 

Single crystals (approximately } in. X 1 in. X 6 in.) of Chempur tin 
(99.986%) were grown from the melt by a technique previously described (2). 
The crystals were grown so that the (001) planes were parallel to the top 
surface of the crystals, and the [110] direction was parallel to the direction of 
growth. This is the preferred direction of growth for tin, and produces stria- 
tions (3) running parallel to the axis of the crystal. Fig. 1 is a Laue back- 
reflection photograph of such a crystal. The orientation difference between 
striations corresponds to a rotation of the lattice about the specimen axis. The 
crystals were deformed in tension on a Hounsfield tensometer. 

The two planes of maximum resolved shear stress in these crystals approxi- 
mated to the (100) and (010) planés, which slip lines indicated to have been 
the active slip planes. The crystals had therefore deformed by conjugate slip. 
Fig. 2 shows the appearance of the top surface of an etched crystal that had 
been deformed in this way. The banded structure is typical. Laue back- 
reflection photographs showed that the bands were regions of different orien- 
tation, the orientation difference being a rotation about the C-axis of the 
lattice in one band with respect to the lattice of an adjacent band. Fig. 3 is a 
Laue back-reflection photograph of a crystal after deformation, taken at the 
boundary between two bands. The splitting of the X-ray spots due to striations 
is a vertical displacement on the photograph corresponding to a rotation about 
the [110] axis. The additional splitting of the X-ray spots corresponding to a 
rotation of the lattice about the [001] direction is due to the deformation bands. 
This orientation difference between the deformation bands increased with 
increasing deformation. For elongations up to 5% each band deformed homo- 
geneously, shown by the absence of asterism in Laue photographs taken on 
single bands. For elongations greater than 10%, asterisms in the X-ray photo- 
graphs showed that the lattice in the bands had become distorted. 

The boundaries of the deformation bands followed the striation boundaries 
wherever possible. This is shown in Fig. 2, taken of a crystal extended 37% to 
show up the bands. The band boundaries did not change appreciably from their 
initial positions with increased deformation, but the greater orientation differ- 
ence between the bands makes them more readily seen and photographed. 
Fig. 4 is a photomicrograph showing the coincidence of a deformation band 
boundary and a striation boundary. In many cases the deformation bands 
were larger than the striations, and more than one striation was contained in a 
single deformation band. The lateral boundaries of the bands were quite 
jagged, since there were no striation boundaries to follow; and in a crystal 
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grown so that there was an area free of striations, the deformation bands 
formed in this region had quite irregular boundaries. 

Barrett and Levenson (1) attribute the formation of deformation bands to 
the operation of different sets of slip systems in different band-shaped regions. 
The observed orientation difference between bands can be explained in these 
terms. The two predominant slip planes are the (100) and the (010) planes. 
Slip on the (100) planes predominates in some bands, and on the (010) planes 
in the other bands. The local lattice rotation thus produced is such that there 
is no change in the mean orientation of the specimen as a whole, a condition 
which is imposed by the grips on the tensile machine. 

The striations running parallel to the crystal axis are attributed to arrays of 
dislocations formed during solidification from the melt (3). These low angle 
boundaries will impede the motion of dislocations through the crystal. The 
dislocation can pass through these boundaries in some cases, and so more than 
one striation can be contained in a single deformation band. Since the defor- 
mation bands will form independently of the existence of available striations, 
the presence of striation boundaries will constitute a natural place for the 
deformation band boundaries to form. Thus, in deformed striated crystals, 
most of the deformation band boundaries coincide with the striation bound- 
aries. 

1. BARRETT, C. S. and LEvenson, L. H. Trans. Am. Inst. Mining Met. Engrs. 135: 327. 
1939; 137: 112. 1940. 


2. CHALMERS, B. Can. J. Phys. 31: 132. 1953. 
3. TEGHTSOONIAN, E. and CHALMERS, B. Can. J. Phys. 29: 30. 1951; 30: 388. 1952. 


RECEIVED May 19, 1953. 

DEPARTMENT OF METALLURGICAL ENGINEERING, 
UNIVERSITY OF TORONTO, 

TORONTO, ONT. 
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Fic. 1. Laue back-reflection X-ray photograph of a striated crystal. 

Fic. 2. A deformed and subsequently etched crystal photographed under oblique illumina- 
tion, X 2. 

Fic. 3. Laue back-reflection photograph of a deformed crystal showing orientation differ- 
ence between deformation bands. 

Fic. 4. Photomicrograph of an etched crystal showing the coincidence of a striation 
boundary and a deformation band boundary, X 500. 








THE CALCULATION OF VIBRATIONAL TRANSITION PROBABILITIES 
OF DIATOMIC MOLECULES 


By R. W. NICHOLLS, W. R. JARMAIN, AND P. A. FRASER 


The intensity of a band of a diatomic molecular band system is controlled 
in part by the square of an integral of the form (2) 


fur R. vs" dr = (v’, R.v"). 
¥,° and ye.) are vibrational wave functions and R.(r) is the electronic 
transition integral. The form of R,(r) is unknown in general and as a first 
approximation is assumed to be slowly varying in the important range of r 
so that it may be replaced by an average value, R,, for the whole band system 
and may thus be taken outside the integrals; (v’, R, v’’) ~ R,(v’, v’’). Relative 
approximate vibrational transition probabilities, p(v’, v’"), may then be defined 
by 

p(v’, v”’) = (v’, v”")?, 

Accordingly, an extended numerical and analytic study of integrals of the 
form (v’, v’’) has recently been made and. applied to molecular transitions of 
astrophysical interest. The integrals (v’, r v’’) and (v’, r? v’’) have also been 
studied for eventual use with experimental intensity measurements in the 
determination of a linear or quadratic form for R,(r). 

Since the Morse potential combines the merits of relative accuracy and 
relative simplicity, the principal work up to the present has been with vibra- 
tional wave functions in Morse potentials. In the Morse potential U;(r) 
= D,{1 — exp—a;,(r — r.;)]?, appropriate to the electronic state i, where the 
symbols have their conventional meaning (2); the v*® normalized vibrational 
wave function is given by (4) 


(9) _ (0) as 1 (Ki—28e-—1) /2 7 Ki—20-1 
i = WN; (exp —}2,) af Lx;—0-1 (2:) 


ignoring the phase factor, and the vibrational energy by 

E{? = (we)i(v + 4) — (were) i(v + 4)’. 
z; = K; exp —a;,(r — re;), LEi—?°>*(z;) is the generalized Laguerre polynomial 
of degree v in 2;, and K; = (we) :/(weXe) i. 

Numerical evaluation of the appropriate exact wave functions and numerical 
integration of (v’, v’’), etc., has been performed for the First and Second 
Positive band systems of nitrogen, and extensive tables obtained. For higher 
quantum numbers the Laguerre polynomials exhibit much cancellation, and 
are very difficult to evaluate. A type of recursion relation has been found that 
relates any one of the polynomials to all those of lower degree. This was used 
since the excessive cancellation was thus avoided. The relation is as follows, 


where L,(z) = LE-°*(z): 
L,(z) = [2 — (K — 20)] Ly-1(2) 
ie! ~ (2k —2)!(v—2 
— (v— 1)(K — 22) & R(k — 1) (: ue _ Ly-2(2) 
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The numerical work, even with this form of recursion relation, is long and 
tedious. 

Consequently an approximation was sought and one has been found that 
allows analytic evaluation of the overlap and other integrals, and which does 
not break down for higher quantum numbers, as would previous approxi- 
mations (5, 1). Essentially, the approximation involves a slight change in 
description of the original vibrational states that allows analytic integration. 
The a; and ae of the Morse potentials appropriate to the electronic states 
between which the transition takes place are each replaced by a mean a, say 
a = (a, + a)/2, and compensating adjustments made in all constants de- 
pending on a and a. Thus normalization and consistency of approximation is 
maintained throughout. This change in description involves a small distortion 
of the potentials and a second order change in the vibrational energy, and 
means that the new potentials and wave functions are probably as adequate 
as the originals for calculation of approximate transition probabilities, bearing 
in mind that the Morse potential is not the ‘true’ potential. 

A working rule, derived partly from theory and partly from experience, 
indicates that the maximum absolute difference between the values of an 
overlap integral calculated using this approximation and by numerical methods 
is roughly given by $| da/a| , where da = (a — a2)/2, when the overlap of 
the potentials is considerable, and by an even smaller fraction of | 6a/a | when 
the overlap is small. This indicates clearly the applicability of the approxi- 
mation to a particular transition, and inspection of tables of molecular data (2) 
shows that the method may be widely used. When not strictly applicable, the 
approximation will at least lead to correct trends. An improvement designed 
for such cases, i.e., | da/a| > 5%, has given results in excellent agreement 
with results of numerical integration. This improvement involves a horizontal 
shift of the ‘new’ potentials to bring the ‘new’ wave functions approximately 
into line with the ‘old’ wave functions, a different shift being necessary for 
each band. This technique approximately doubles the work, but the method 
is still much more rapid than the numerical method. 

The only limiting factor on the extent of a table of (v’, v’’)? calculated using 
this approximation is the digit capacity of the desk calculator available, 
because of the wholesale cancellation that arises in evaluation of the formulae, 
the cancellation increasing with (v’ + v’’). This cancellation is intrinsic in the 
problem, and can be traced back to the Laguerre polynomials. For example, 
using a 10 digit capacity calculator, one may calculate (v’, v’’) for values of 
(v’ + v’’) < 6 for the First Positive system of nitrogen, while for the Vegard- 
Kaplan bands, the table may be extended to (v’ + v’’) < 9, the cancellation 
being less severe owing to the smallness of the (0, 0) overlap. Formulae have 
been derived for bands which have one of the quantum numbers equal to 
0, 1, 2, or 3, in which the cancellation is carried out algebraically, the remainder 
exhibiting no cancellation. Work is continuing on these ‘remainder’ formulae 
for other quantum numbers since their use allows great extension of the tables. 

A sample formula, resulting from the approximation made, before the 
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cancellation has been algebraically carried out, is shown below mainly to 
demonstrate the type of numerical work involved in final evaluation. For 


example, 
(0, 2) Net 2 p, (K2— 3) 4,2 (Ks — 3)(Ka — 4) 


(0,0) ~ WU! — 2 2) + Ry — 8) J’ 





with 


_ (1) ay bt ee 
(0,0) = (E) te V(T(Ki — 1) T(Ke — 1)] * 


Taking the ratio of (v’, v’’) to (0, 0) provides a great simplification in compu- 
tation; the only gamma functions that have to be evaluated are then those 
in the (0, 0) integral. Even ‘this can be avoided. Similar formulae give 
(v’, r v’’)/(0, r 0) and (v’, r? v’’)/(0, r? 0). : 


TABLE I 
No, B'ly— A*Z{, First Positive 
P(r’, »””) = (v', 0”)? 








0 1 | 2 





v ye vi 
0 33; 32, 19, 
(.34p) (.325) (19>) 
, Pe 005 10, 
| (405) | (.00) (.105) 
| 
ee ae Ea 
| (200) | (%) | ¢MD 
| 





p-Values (i.e., (v’, v’’)?) have been computed using this approximation 
for 11 transitions some of which are in process of publication (3). In the 
table are shown a few such p-values for the First Positive System of Nao, 
B*II, — A*Z}. This table contains but a portion of the triangle of 28 bands, 
(v’' + v’’) < 6, that have been obtained by this method, as mentioned above. 
Values obtained by numerical integration of ‘exact’ wave functions are also 
shown, for comparison, in parentheses. Agreement is not always so completely 
satisfactory, but is most often completely adequate for the purposes to which 
these numbers are put. A program of computing p-values by this approximate 
and rapid method for many transitions of astrophysical interest is now under 
way (3), the justification being the validity of the approximation as shown in 
those transitions for which comparison with numerical integrations has been 
possible. Molecules so far studied are: Nz, Not, NO, O2+, CN, Oc, and C2; and 
those under course of study include: OH, CH, CO, and CO?. 


The present work is part of a continuing experimental and theoretical re- 
search program into molecular excitation processes supported by Contract 
No. AF 19(122)-470 of the Air Force Cambridge Research Centre. 
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LETTERS TO THE EDITOR 





Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of issue. No proof will be sent to the authors. 


“‘Phosphorescent”’ Effect in Lead Selenide 


During recent work on the photoconductive properties of lead selenide, it was observed 
that this substance exhibited a very strong ‘‘phosphorescent”’ effect when cooled to liquid air 
temperatures. The lead selenide was obtained in the form of a thin layer by evaporating it 
in vacuo onto the end of the inner wall of a small Dewar. Aquadag electrodes were previously 
painted on this surface so that only a small area 8.X 1 mm.? was used for the experiments. 
Connections to the electrodes were made by means of 1 mm. diameter tungsten leads which 
were brought out at the top of the Dewar. The outer wall consisted of a periclase window 
opposite the electrodes, and a graded seal 1 in. in diameter. Periclase is transparent to infrared 
radiation out to 8 microns and may be sealed directly to soft glass (3). The graded seal was 
used for joining this soft glass to the pyrex glass at the top of the Dewar. 

Fig. 1 gives the results of measurements on the conductivity of the layer plotted on a log 
scale against the inverse of the temperature. To obtain the various temperatures a small 
quantity of mercury was placed in the Dewar. The mercury was then frozen using liquid air 
to 90° K., and then gradually allowed to warm up by means of a small air jet. The temperature 
was measured using a copper-constantan thermocouple immersed in the mercury. 
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Fic. 1. Conductivity of lead selenide. 


The curve ABC is that normally obtained for a semiconductor such as lead selenide, the 
point C corresponding to 90° K. The first part of the curve AB is linear, conforming to the 
usual equation ¢ = oo exp —e/2KT. The curved part BC may be ascribed to impurities and 
photoconductive effects due to room temperature radiation (2). On illuminating the layer with 
a small tungsten lamp the conductivity changes to the point E£. If this illumination is now 
removed the conductivity returns to the point D, and not to C as would be expected. Allowing 
the layer to warm up gradually, the conductivity curve follows DFB, joining the original curve 
at B, and on recooling the conductivity retraces the curve BC. Along FD, lead selenide there- 
fore behaves as a metal. Keeping the temperature constant at 90° K. the conductivity is 
observed to remain at D for at least one hour without any noticeable change. This effect does 
not appear to depend on the amount of impurities present since it has been observed in layers 
of very different sensitivity. A similar effect has been observed for lead sulphide by Chasmar 
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and Gibson (1) who ascribed it as due to transitions to metastable levels by the photoelectron s. 


The curve DFB shows the way in which these trapped electrons are thermally released, if this 
is assumed to be the correct interpretation of this effect. 

This note is published by permission of the Defence Research Board. The author wishes to 
thank Dr. A. M. Crooker for the interest he has shown in this problem. 


1, CHASMAR, R. P. and Gipson, A. F. Proc. Phys. Soc. (London), B, 64: 595. 1951. 
2. Watts, B. N. Proc. Phys. Soc. (London), A, 62: 456. 1949. 
3. Younc, A.S. J. Sci. Instr. 28: 207. 1951. 
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